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ABSTRACT 


NONNEGATIVE  MATRICES  WITH  PRESCRIBED  SPECTRA 

Patricia  D.  Egleston 

Nonnegative  matrices  are  rectangular  arrays  of  nonnegative  real  numbers.  Non¬ 
negative  matrices  are  prevalent  in  many  areas  of  study.  Statistics,  Economics,  and 
Chemistry  are  examples  of  disciplines  that  use  nonnegative  matrices.  In  this  disserta¬ 
tion  we  examined  and  produced  new  results  on  nonnegative  matrices  with  prescribed 
spectra.  This  problem  is  called  the  nonnegative  inverse  eigenvalue  problem  (NIEP). 
The  NIEP  asks  which  lists  of  n  complex  numbers  occur  as  the  spectrum  of  an  n  x  n, 
entry-wise  nonnegative  matrix.  A  recent  paper  by  Leal-Duarte  and  Johnson  stated, 
“this  problem  has  attracted  considerable  attention  over  50  years  and,  despite  many 
exciting  partial  results,  remains  quite  unresolved.”  In  this  dissertation  we  solved  the 

NIEP  for  certain  lists  of  four  complex  numbers. 

Since  the  general  NIEP  is  difficult  to  resolve,  two  natural  variations  were  consid¬ 
ered.  The  real  nonnegative  inverse  eigenvalue  problem  (RNIEP)  asks  which  lists  of  n 
real  numbers  occur  as  the  spectrum  of  a  nonnegative  n  x  n  matrix.  The  symmetric 
nonnegative  inverse  eigenvalue  problem  (SNIEP)  asks  which  lists  of  n  real  numbers 


IV 


occur  as  the  spectrum  of  a  nonnegative  n  x  n  matrix.  The  symmetric  nonnegative 
inverse  eigenvalue  problem  (SNIEP)  asks  which  lists  of  n  real  numbers  occur  as  the 
spectrum  of  a  symmetric  nonnegative  n  x  n  matrix.  We  resolved  the  SNIEP  in  most 
of  the  cases  for  lists  of  five  real  numbers  and  lists  of  six  real  numbers.  In  a  paper 
by  Johnson,  Laffey,  and  Loewy,  it  was  shown  that  the  RNIEP  and  the  SNIEP  are 
different  problems.  We  have  shown  that  the  smallest  list  size  in  which  the  RNIEP 
and  the  SNIEP  are  different  is  five. 

A  matrix  A  =  [oy]  is  said  to  be  subordinate  to  a  given  graph  G  if  whenever  vertex 
Vi  and  Vj  are  not  adjacent  then  entry  =  0.  We  resolved  some  cases  of  the  NIEP 
where  the  realizing  nonnegative  matrix  was  subordinate  to  a  bipartite  graph. 

In  conclusion,  the  results  in  this  dissertation  have  taken  us  several  steps  closer  to 
the  resolution  of  one  of  the  outstanding  problems  in  matrix  theory. 
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CHAPTER  1 


Introduction 

We  will  be  studying  the  nonnegative  inverse  eigenvalue  problem  (NIEP).  The 
NIEP  is  to  determine  necessary  and  sufficient  conditions  for  a  list  of  n  complex 
numbers  to  be  the  spectrum  of  a  nonnegative  n  x  n  matrix.  While  this  problem  is 
simple  and  very  easy  to  understand,  it  remains  unsolved.  This  problem  has  been  the 
subject  of  much  study  over  the  last  50  years  giving  various  partial  solutions.  In  this 
dissertation  we  resolve  some  of  the  special  cases  of  the  NIEP . 

The  next  seven  chapters  of  this  dissertation  present  a  potpourri  of  results.  Because 
of  the  diversity  of  material  presented,  we  dedicate  Chapters  2,  3,  and  5  to  presenting 
background  information  regarding  the  subproblems  for  which  we  have  new  results. 
We  present  our  new  results  in  Chapters  4,  6,  7,  and  8. 

We  begin  with  a  discussion  of  nonnegative  matrices  with  various  known  pertinent 
results  and  theorems  in  Chapter  2.  In  Chapter  3,  we  examine  the  development  of  the 
NIEP.  In  this  chapter  we  give  the  results  of  Suleimanova,  Perfect,  Loewy,  and  London 
regarding  their  solutions  to  subproblems  of  the  NIEP.  Here  we  see  the  solutions  to  the 
2  X  2  and  3  X  3  NIEP  as  well  as  the  trace  zero  solutions  to  a  list  of  four  or  five  complex 
numbers.  We  then  examine  the  history  of  the  real  and  symmetric  nonnegative  inverse 
eigenvalue  problems  in  Chapter  5.  Here  we  have  the  results  of  Johnson  et.al.  regarding 
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the  difference  between  the  two  problems.  We  present  the  sufficient  conditions  of 
Salzmann,  Kellogg,  Fiedler,  Borobia,  and  Radwan  regarding  the  real  nonnegative 
inverse  eigenvalue  problem  (RNIEP)  or  the  symmetric  nonnegative  inverse  eigenvalue 
problem  (SNIEP),  and  conclude  with  a  brief  discussion  on  the  Soules  construction  of 
a  symmetric  nonnegative  matrix. 

We  present  new  results  in  Chapters  4,  6,  7,  and  8.  In  Chapter  4,  we  present  the 
partial  solution  to  the  NEIP  in  regard  to  a  list  of  four  complex  numbers  with  nonzero 
trace  with  our  main  result  in  Theorem  14.  Chapter  6  contains  the  analysis  of  a  list  of 
five  real  numbers  and  shows  that  this  is  the  smallest  list  in  which  the  two  problems, 
the  RNIEP  and  SNIEP,  are  different.  This  is  presented  in  Corollary  4.  In  Chapter  7, 
we  analyze  the  list  of  six  real  numbers  relative  to  the  RNIEP  and  SNIEP  and  resolve 
many  cases.  Finally,  the  last  chapter  of  results,  Chapter  8,  contains  a  partial  solution 
to  the  question  of  when  a  list  of  complex  numbers  that  sums  to  zero  is  the  spectrum 
of  a  nonnegative  matrix  subordinate  to  a  bipartite  graph. 
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CHAPTER  2 


Nonnegative  Matrices 

Nonnegative  matrices  are  used  in  many  areas  of  pure  and  applied  mathematics. 
Numerical  analysis,  statistics,  molecular  chemistry,  and  economics  are  examples  of 
areas  where  nonnegative  matrices  appear.  For  instance,  the  Leontief  input-output 
matrix  in  economics  is  a  nonnegative  matrix.  In  combinatorics,  the  adjacency  matrix 
of  a  graph  is  nonnegative.  A  matrix  where  each  entry  represents  the  probability  of 
occurrence  is  nonnegative. 

Throughout  this  dissertation  we  will  restrict  our  attention  to  the  matrices  that  are 
square  with  real  or  complex  entries.  We  will  define  a  matrix  A  to  be  a  nonnegative 
matrix  if  every  entry  of  A  is  real  and  nonnegative,  and  a  matrix  A  to  be  a  positive 
matrix  if  every  entry  of  A  is  real  and  positive.  We  will  denote  a  nonnegative  matrix 
A  by  A  ^  0  and  a  positive  matrix  A  by  A  >  0.  An  order  relation  between  two  real 
matrices,  A  and  B,  will  be  defined  as  follows.  We  write  A^HifA  —  and 

A>  jBifA  —  jB>0.  Furthermore,  given  A  =  [ajj],  we  define  |A1  to  be  the  matrix 
containing  the  absolute  values  of  the  entries  of  A.  That  is,  jA]  =  [[ciijl].  Note  that 
for  any  A,  \A\  ^  0. 

The  notation  Mn  will  denote  the  set  of  all  n  x  n  matrices  with  complex  entries. 
Given  a  matrix  A  e  Mn  and  x  €  C”,  we  consider  solutions  to  the  equation  Ax  =  Ax, 
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where  x  ^  0  and  A  is  a  complex  number.  If  a  number  A  and  a  nonzero  vector  x  satisfy 
this  equation,  then  A  is  called  an  eigenvalue,  of  A  and  x  is  called  an  eigenvector  of  A 
associated  with  A.  Notice  that  the  collection  of  eigenvectors  associated  with  A  and  the 
zero  vector  form  a  subspace  of  C”  called  the  eigenspace  of  A.  Associated  with  the 
matrix  A  is  an  degree  polynomial  p(A)  =  det(A/  —  A).  The  eigenvalues  of  A  may 
also  be  defined  as  the  n  roots  of  the  polynomial  p  (A).  We  call  this  polynomial  the 
characteristic  polynomial  of  A.  The  set  of  eigenvalues  of  A  is  called  the  spectrum  of 
A  and  is  denoted  cr(A).  The  spectral  radius  of  A,  denoted  p{A),  is  the  nonnegative 
real  number  defined  by  p  (A)  =  max{\X\  :  A  6  cr  (A)}. 

In  this  dissertation  we  will  use  the  following  vector  norms  on  C”. 

1.  the  sum  norm  or  l\  norm  on  C”  given  by  =  [xil  +  •  •  •  +  la^nli 

2.  the  max  norm  or  l^o  norm  on  C”  given  by  Ija^Hoo  =  max{\x\\ , . .  • ,  |a^n|}> 

3.  the  Euclidean  norm  or  I2  norm  on  C"  given  by  ||a:||2  =  (lCi=i  l^il  )  • 

The  unit  ball  {x  :  |la;||  <  1}  of  C"  is  compact  in  any  of  the  vector  norms. 

Associated  with  each  vector  norm  on  C"  is  a  natural  matrix  norm  on  Mn  that  is 
induced  by  the  vector  norm  as  follows.  Define  1|  •  |1  on  Mn  by  ||All  =  moj;iix||=il|Ax|l. 
One  such  norm  is  the  spectral  norm  or  the  operator  norm  induced  by  the  Euclidean 
norm  on  C".  This  spectral  norm  is  denoted  |1A||2  and  is  equivalently  defined  as 
IIAII2  =  max{y/X  :  X  e  a  (A* A)}.  It  is  known  that  for  any  matrix  norm  1|A1|  ^  p  (A) 
and  the  spectral  radius  formula  is  given  by  p  (A)  =  Ito  HA^'ll  * . 
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When  investigating  the  spectrum  of  a  given  matrix,  it  is  often  useful  to  find 
special  properties  of  the  spectrum  without  finding  all  the  eigenvalues.  Note  that  for 
the  following  real  matrix 


0  1 

-1  0 


the  spectrum  c  (A)  =  {i,  —i}  is  not  real.  If  we  consider  another  example 


then  the  spectrum  a{B)  =  {-2,  -3}  is  real  and  negative.  In  neither  case  does  the 
spectral  radius  p  (A)  =  1  and  o  {B)  =  3  belong  to  the  respective  spectrum. 

Considering  the  fact  that  real  matrices  may  have  only  complex  eigenvalues  or  only 
nonpositive  eigenvalues,  it  is  remarkable  to  note  that  a  nonnegative  matrix  must  have 
at  least  one  positive  eigenvalue.  More  precisely,  the  spectral  radius  of  a  nonnegative 
matrix  must  be  an  eigenvalue. 


We  begin  the  discussion  of  nonnegative  matrices  with  the  result  of  0.  Perron  [26]. 


In  1907,  Perron  studied  positive  square  matrices.  He  found  that  a  positive  square 
matrix  always  has  its  spectral  radius  as  an  element  of  the  spectrum.  Therefore  we 
call  the  spectral  radius  of  a  positive  or  nonnegative  matrix  A  the  Perrm  root.  An 
eigenvector  associated  with  the  Perron  root  is  called  the  Perron  vector.  The  following 
due  to  Perron  [26]  states  the  conditions  on  both  the  Perron  root  and  the  Perron  vector. 
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Theorem  1  (Perron’s  Theorem).  [11]  IfAeMn  and  A  >  0,  then 

1. piA)>0; 

2.  p  (i4)  is  an  eigenvalue  of  A; 

3.  there  is  an  x  E  C”  with  x  >  0  and  Ax  —  p  (A)  x; 

I  p{A)  is  an  algebraically  (and  hence  geometrically)  simple  eigenvalue  of  A; 

5.  |A|  <  p  (A)  for  every  eigenvalue  p  (A),  that  is,  p  (A)  is  the  unique  eigenvalue 
of  maximum  modulus. 


Notice  that  in  addition  to  the  Perron  root  being  an  element  of  the  spectrum,  both 
the  Perron  root  and  the  Perron  vector  of  a  positive  matrix  must  be  positive. 

G.  Frobenius  extended  many  of  Perron’s  results  to  irreducible,  nonnegative  ma¬ 
trices.  Recall  that  a  matrix.  A,  is  said  to  be  reducible  if  there  exists  a  permutation 


matrix  P  e  Mn  such  that 


P'^AP 


B  C 
0  D 


where  B  and  D  are  square  matrices  and  0  is  a  zero  matrix.  An  irreducible  matrix  is 
therefore  a  matrix  which  is  not  reducible.  Frobenius  [9]  proved  that  for  an  irreducible, 
nonnegative  matrix.  A,  p  (A)  is  positive  and  belongs  to  the  spectrum  of  A.  Frobenius 
also  proved  that  the  Perron  vector  remains  positive.  All  of  these  results  are  given  in 
the  following  theorem  called  the  Perron-Frobenius  Theorem. 
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Theorem  2  (Perron-Frobenius  Theorem).  [11]  Let  AeMn  and  suppose  that  A 
is  irreducible  and  nonnegative.  Then 

1.  P  (A)  >  0; 

2.  p  {A)  is  an  eigenvalue  of  A; 

3.  there  is  an  x  €.  C"  with  x  >  0  such  that  Ax  =  p  (A)  x;  and 

I  p  (yl)  is  an  algebraically  (and  hence  geometrically)  simple  eigenvalue  of  A. 

The  result  (5)  of  Perron’s  Theorem  is  no  longer  true  when  A  ^  0.  When  A  ^  0 
and  irreducible  we  have  the  following  result. 

Theorem  3.  [11]  Let  A  G  Mn-  Suppose  A  is  nonnegative  and  irreducible  and  suppose 
the  set  S  =  {Xi  =  p{A),X2,...,  Afc}  of  eigenvalues  of  maximum  modulus  p  (A)  has 
exactly  k  distinct  elements.  Then  each  eigenvalue  Xi  E  S  has  algebraic  multiplicity 
one,  andS  =  {e"-^p{A):t  =  0,l,...,k-  1}.  Moreover  if  X  is  any  eigenvalue  of  A, 
then  A  is  an  eigenvalue  for  allt  =  0,  1, . . . ,  fe  —  1. 

Notice  that  the  necessary  requirements,  (1)  the  Perron  root  being  positive  and 
belonging  to  the  spectrum  and  (2)  the  Perron  vector  being  positive,  are  maintained 
as  in  Perron’s  Theorem  but  now  the  matrix  no  longer  must  be  strictly  positive.  This 
raises  a  natural  question.  Can  we  remove  the  irreducibility  condition  from  the  Perron- 
Frobenius  Theorem?  In  order  to  investigate  this  question,  we  will  need  Theorem  4. 


7 


Theorem  4.  [11]  Let  A,  B  £  M„.  If  \A\  ^  B,  then 

p{A)  ^  p(|A|)  ^  p{B).  In  particular  ifO^A-^B,  then  p{A)  ^  p{B). 

Proof  For  every  m  =  1,2,...  we  have  \A^\  <  \A\^  <  B^.  Hence  for  the  spectral 
norm,  ||A”^1|2  ^  1|  I^P  lb  <  ll^'^lb  for  all  m  =  1, 2, . . ..  By  using  the  spectral  radius 
formula  we  deduce  that  p{A)  ^  p(|^l)  ^  p{^)-  If  0  ^  A  ^  B,  then  p{A)  ^ 
p{S).  □ 

In  the  Perron-Frobenius  Theorem,  the  requirement  of  the  matrix  A  being  irre¬ 
ducible  as  well  as  nonnegative  led  to  the  necessary  condition  that  the  Perron  vector 
is  positive.  When  the  matrix  A  is  just  nonnegative,  the  condition  becomes  that  the 
Perron  vector  is  nonnegative.  The  condition  that  p  (A)  be  an  element  of  the  spectrum 
still  holds  true.  These  facts  are  proved  in  Theorem  5. 

Theorem  5.  [11]  //  A  €  and  A  ^  0,  then  p  (A)  is  an  eigenvalue  of  A  and  there 
is  a  nonnegative  vector  x  ^  0,  x  ^  0,  such  that  Ax  =  p  (A)  x. 

Proof.  Let  e  >  0.  Define  A{e)  =  [aij  -f  e]  >  0.  Let  x  (e)  represent  the  Perron  vector 
of  A  (e)  such  that  a;  (e)  >  0  and  x  (e)^  =  1.  As  the  set  of  vectors  {x  (e) :  e  >  0}  is 

izzl 

contained  in  the  compact  set  {x  :  x  £  ||ic||i  ^  1},  there  is  a  monotone  decreasing 

sequence  ei,  £2)  •  •  •  with  limSfe  =  0  such  that  limx  {sk)  =  x  exists.  Since  x  (sk)  >  0 

^  k-^oo  k-^oo 

for  all  A:  =  1,  2,  . . . ,  it  must  be  that  x  =  lima:  (e/t)  ^  0.  Note,  a;  =  0  is  impossible 

fe-^oo 

since  £  a:*  =  lim  rc  {sk)i  =  1.  By  Theorem  4,  we  have  p  (A  (e*))  >  p  (A  (£*+0)  ^ 

i=:i  i=l 

■  •  •  ^  p  (A)  for  all  A:  =  1,  2,  ... .  So  the  sequence  of  real  numbers  {p  (A  (£fc))}fc=i,2,...  is 
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a  monotone  decreasing  sequence.  Therefore  p  =  ^li^p  (-4  (efe))  exists  and  p  ^  p(j4). 
Now  because  Ax  =  lim  A  (sk)  x  (sk)  -  lim  p  {A  {sk))  x  {sk)  =  px  and  x  ^  0,  it  follows 

k^oo  fc->oo 

that  p  is  an  eigenvalue  of  A  so  p  ^  p  (A).  Thus  p  =  p  (A).  ^ 


The  following  result  will  be  used  in  Chapters  4  and  8.  We  state  it  here  for  reference. 
Define  a  matrix  B  as  follows 

0  Bi,2  0  •  •  •  0 

0  0  ^2,3  •  •  •  0 

B  = 

0  0  •  •  •  0  Bk-i,k 

[bm  0  •••  0  0 

where  the  block  Bjj+i  is  rij  x  n^+i  for  j  =  1, . . . ,  A:  -  1  and  Bk,i  is  %  x  nj  to  be  m 
super  diagonal  farm  [22].  This  form  has  a  special  relationship  with  its  eigenvalues 

as  Theorem  6  shows. 


Theorem  6.  [22]  Let  A  be  an  n  x  n  complex  matrix  in  super  diagonal  form  and 
suppose  that  wi, . . . ,  Wm  are  the  nonzero  eigenvalues  of  the  product  Ai,2A2,3  •  *  •  Ak,i- 
The  spectrum  of  A  consists  of  the  n-km  zeros  and  the  km  k^^  roots  of  the  numbers 


j  •  *  *  5 

The  discoveries  of  Perron  and  Probenius  have  given  some  powerful  necessary  con¬ 
ditions  in  the  study  of  nonnegative  matrices.  The  spectrum  of  a  nonnegative  matrix 
must  contain  the  Perron  root  which  is  positive  in  all  cases  except  the  zero  matrix 
where  it  is  zero,  and  the  Perron  vector  must  be  nonnegative.  For  example  we  can 
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deduce  that  a  list  of  negative  numbers  will  never  be  the  spectrum  of  any  nonnegative 
matrix.  This  is  because  the  Perron  root  is  not  part  of  the  set.  So  addressing  the 
special  properties  of  nonnegative  matrices  has  become  a  source  of  interesting  mathe¬ 
matical  problems,  we  will  explore  some  of  these  ideas  in  the  remaining  chapters. 
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CHAPTER  3 


Nonnegative  Inverse  Eigenvalue  Problem 

3.1  Beginning  of  the  Nonnegative  Inverse  Eigenvalue  Problem 

In  1938  Kolmogorov  asked  when  is  a  given  complex  number  an  eigenvalue  of  a 
nonnegative  matrix.  When  there  is  no  restriction  on  the  order  of  the  matrix,  then 
any  complex  number  is  an  eigenvalue  of  some  nonnegative  matrix.  A  simple  solution 
comes  from  Mine  [22]:  if  a +  ^i,  where  a  and  P  are  any  real  numbers,  is  the  complex 
number  in  question,  then  the  circulant  matrix 

|al  1^1  0  0 

0  bl  1^1  0 

A  {a,  P)  = 

0  0  [al  \I3 

m  0  0  |« 

will  have  jal  ±  |/5|  and  1q;|  ±  1^1  i  as  eigenvalues.  Now  let  B  (a,  be  the  8  x  8  matrix 

0  A{a,p) 

B{aJ)  = 

A{a,P)  0 

Theorem  6  shows  that  the  eigenvalues  of  B  {a,  13)  are,  ±  \a\  ±  j^]  and  ±  |q:|  ±  \P\  i. 
Note  that  one  of  the  last  four  eigenvalues  is  the  given  a  +  /3^. 

It  can  be  easily  shown  that  a  single  complex  number  with  imaginary  part  not 
equal  to  zero  cannot  be  an  eigenvalue  of  nonnegative  matrices  of  order  1  or  order  2. 
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In  the  case  of  an  order  1  nonnegative  matrix,  there  is  exactly  one  eigenvalue,  namely 
the  Perron  root.  Since  any  complex  root  of  a  characteristic  polynomial  with  real 
coefficients  must  occur  with  its  conjugate  pair,  in  the  case  of  an  order  2  nonnegative 
matrix  we  know  that  if  one  root  is  complex  then  so  is  the  second.  But  one  root 
must  be  the  real  Perron  root  and  so  the  second  must  be  real  also,  thereby  eliminating 
the  possibility  of  a  complex  root.  Therefore,  the  smallest  dimension  for  which  a 
complex  number  with  nonzero  imaginary  part  can  be  an  element  of  the  spectrum  of 
a  nonnegative  matrix  is  n  =  3.  The  general  problem  as  set  above  by  Kolmogorov 
without  any  restriction  on  the  matrix  thus  has  a  complete  solution.  To  extend  this 
problem  it  becomes  necessary  to  restrict  the  question  to  certain  classes  of  nonnegative 
matrices;  e.g.  stochastic  matrices,  doubly  stochastic  matrices,  nonnegative  circulants. 
Recall  that  a  stochastic  matrix  is  a  nonnegative  matrix  with  all  row  sums  equal  to 
one,  and  a  doubly  stochastic  matrix  is  a  stochastic  matrix  with  the  added  restriction 
that  all  column  sums  equal  one.  History  credits  Dmitriev  and  Dynkin  [5,  6]  for  partial 
results  in  the  case  of  stochastic  matrices  and  Karpelevich  [14]  for  completely  solving 

it. 

3.2  Growth  of  the  Nonnegative  Inverse  Eigenvalue  Problem 

In  1949  Suleimanova  [32]  extended  Kolmogorov’s  question  to  what  is  now  called 
the  nonnegative  inverse  eigenvalue  problem  (NIEP).  The  NIEP  is  to  determine  neces¬ 
sary  and  sufficient  conditions  for  a  list  ofn  complex  numbers  to  be  the  spectrum  of  a 
nonnegatwe  nxn  matrix.  Suleimanova  answered  this  question  when  the  list  of  length 
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n  consists  of  only  real  numbers  where  exactly  one  number  is  positive.  Suleimanova’s 
sufficient  condition  for  a  =  {Ai,  A2,  ■  •  ■ ,  An)  where  Ai  >  0  ^  A2  ^  ^  A„  to  be  the 

spectrum  of  a  nonnegative  matrix  is  that  ^  Aj  ^  0.  We  will  prove  Suleimanova  s 
result  in  Chapter  5.  But  note  for  now  that  this  result  can  be  extended  to  any  list 
of  real  numbers  a  =  {Ai,  A2,  A3, ,  A„}  such  that  Ai  +  XI  ^  Perron 

\j<0 

root  dominates  the  sum  of  all  the  negative  numbers.  This  extension  is  due  to  Perfect 
and  is  discussed  further  in  Chapter  5. 

3.3  Necessary  Conditions  for  the  NIEP 

If  a  matrix  A  is  nonnegative  then  its  trace  must  be  nonnegative.  Since  powers 
of  A  are  nonnegative  so  are  their  traces.  This  necessary  condition  is  described  as 
follows;  Let  A  be  an  n  x  n  nonnegative  matrix.  Suppose  cr  =  {Ai, ... ,  An}  is  the  list 

n 

of  eigenvalues  (repeats  included)  of  A.  Then  define  the  moment  Sk  =  Aj  — 
trace  (A*^)  for  k  —  I,  2,  . . ..  The  necessary  condition  can  be  stated  as  Sfc  ^  0  for 

fc  =  1,2,.... 

Since  the  characteristic  polynomial  of  a  nonnegative  matrix  has  only  real  coeffi¬ 
cients,  the  complex  eigenvalues  must  occur  in  conjugate  pairs.  This  is  indicated  by 

the  necessary  condition  a  =  a  where  a  =  {Ai, . . . ,  An}- 

The  third  condition  due  to  Johnson,  Loewy  and  London  is  described  in  the  fol¬ 
lowing  theorem.  The  condition  appeared  first  in  Loewy  and  London  [20].  Johnson 
independently  found  the  condition  in  [12].  We  will  henceforth  call  this  necessary 
condition  the  JLL  condition. 
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Theorem  7.  [20]  Let  a  =  {Ai,  A2,  •  •  • ,  A„}  be  the  set  of  eigenvalues  of  a  nonnegative 
n  X  n  matrix  A  —  [oy]-  Then  ^  n'^  ^Skm  fo^  k,m  —  1,  2,  . . .. 

Proof.  Consider  the  diagonal  matrix  D  =  diag  (an,  022,  ^33)  •  •  •  >  o-nn)  and  let  A  = 

C  +  D.  Since  A  is  nonnegative  so  are  C  and  D.  It  follows  that  A^  -  C'^  -  for 
m  =  1,  2,  ...  is  nonnegative  too.  This  gives  trace  {A^)  >  trace  (C™)  +  trace  (D^)  = 
trace  {C^)  +  ^a^. 

/  71  \  ”*  n 

Then  by  Holder’s  inequality  we  have,  (  Eoiij  <  Therefore, 

trace  {A^)  =  Sm 

n 

^  trace  {C'^)  +  d!A 

i=i 

n 

^  n^~Hrace  {C'^)  + 

i=l 

(n  \  TTi 

j 

=  n^~^trace  (C^)  +  (si)*” 

Which  results  in  n^'^Sm  -  >  n^-Hrace{C^).  But  since  trace  {C^  ^  0 

this  implies  n'^-hm  -  sf  ^  Q  which  then  gives  n^-^m  ^  sf  for  m=  1,  2,  .... 
Applying  this  inequality  (which  is  true  for  fc  =  1)  to  the  nonnegative  matrix  A^  we 

get  <  n'^~'^Skm  for  all  k,m  =  l,2,  .... 
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Given  a  list  a  =  {Ai,  A2, .  •  • ,  A„},  we  summarize  the  four  necessary  conditions  for 
the  NIEP  below. 

1.  If  />  =  moxi^i<n{|Ail  ;  A^  6  a},  then  pea  (Perron  root  condition), 

2.  a  =  cr^  i.e.,  cr  must  be  closed  under  conjugation, 

3.  ^  Af  ^  0  for  A:  =  1, 2, . . i.e.,  all  moments  are  nonnegative,  and 

i-l 

4.  the  JLL  condition;  ^  for  all  k,m  — 1,2,  — 

We  say  that  a  list  cr  =  {Ai, . . . ,  A„}  of  length  n  is  realizable  by  an  n  x  n  nonnegative 
(resp.  symmetric  nonnegative)  matrix  if  there  exists  a  matrix  A  ^  0  (resp.  symmetric 

matrix  A  ^  0)  such  that  a  [A)  is  the  given  list  a. 

It  is  easy  to  give  examples  that  conditions  (1),  (2),  and  (3)  are  not  sufficient. 

^  By  the  Perron-Frobenius  Theorem,  since  the  maximal 

element  is  not  simple,  the  realizing  matrix  must  be  reducible.  However,  the  given  list  a 
cannot  be  divided  into  two  sublists  each  satisfying  all  the  necessary  conditions  (1),  (2), 
and  (3). 

Now  condition  (4)  does  not  follow  from  conditions  (1),  (2),  and  (3).  We  see  this 
by  an  example.  Suppose  a  =  {^^2,^,-i},  then  (1)  -  (3)  are  satisfied  but  (4)  fails  for 

k  =  l  and  m  =  2. 

Forn  ^  4  condition  (4)  is  not  sufficient.  To  see  this  consider  a  =  {\/2,  a/2,  i,  “O- 
It  is  easy  to  check  that  a  satisfies  conditions  (1)  -  (4).  If  a  is  the  set  of  eigenvalues 
of  a  nonnegative  matrix  A,  then  according  to  the  Perron-Frobenius  Theorem  A  is 


15 


reducible.  Hence  {\/2,  i,  -i}  must  be  realizable  by  some  nonnegative  matrix.  But  we 
have  already  seen  that  is  not  possible. 

3.4  Stochastic  Inverse  Eigenvalue  Problem 

Many  subproblems  have  emerged  from  the  NIEP  because  of  its  complexity.  The 
stochastic  inverse  eigenvalue  problem  is  one  of  these  which  is  to  determine  necessary 
and  sufficient  conditions  for  a  list  of  n  complex  numbers  to  be  the  spectrum  of  a 
stochastic  n  x  n  matrix.  The  next  result  shows  that  when  there  exists  a  nonnegative 
matrix  A  with  spectrum  a  =  {Ai,  A2,  As, . . . ,  A^}  and  positive  Perron  vector,  then 
A^^A  is  diagonally  similar  to  a  stochastic  matrix. 

Theorem  8,  [22]  If  A  is  a  nonnegative  matrix  with  Perron  root  Xi  and  a  positive 
Perron  vector,  then  Xf^A  is  diagonally  similar  to  a  stochastic  matrix. 

Proof.  Let  a:  =  {xi,  X2,  be  the  positive  Perron  vector  of  A.  If  D  = 

diag{xi,  X2,  Xn),  then  x  =  Du  where  u  =  (1,  1,  1)^-  As  1  is  the  Per¬ 

ron  root  of  (Ar^A)  D  =  X[^D-^AD,  we  have  X[^D-^ADu  =  X^^D-'^Ax  = 
Xf^D'^Xix  =  D~^x  =  u.  Therefore  u  is  the  Perron  vector  of 

Xf^D~^AD  =  D~^  (^r^A)  D.  And  since  a  square  nonnegative  matrix  A  is  stochastic 
if  and  only  if  u  =  (1,  1,  . . . ,  1)^  is  an  eigenvector  corresponding  to  the  Perron  root 
1  of  A,  we  have  the  desired  result.  ^ 

Johnson  discusses  the  stochastic  inverse  eigenvalue  problem  in  [12].  He  calls  it  the 
row  stochastic  eigenvalue  problem.  The  Perron  root  of  a  row  stochastic  matrix  is  one; 
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and,  if  the  Perron  root  is  normalized  to  be  one  in  the  NIEP,  then  the  row  stochastic 
eigenvalue  problem  is  equivalent  to  the  NIEP.  (The  list  {0,0,...,  0}  is  omitted  as  we 
have  the  zero  matrix  as  a  solution.)  If  A  is  an  arbitrary  n  x  n  nonnegative  matrix, 
then  A  is  permutationally  similar  to 

0 

P'^AP  = 

*  Am 

where  A*  is  rij  x  n,,  i  =  1, . . . ,  m,  and  irreducible.  The  Ai  are  called  the  irreducible 
components  of  A  and,  of  course,  if  A  itself  is  irreducible,  then  m  =  l.  Now  a  (A)  = 
Dill  is  the  same  as  cr  (Ai  ©  •  •  •  ©  A,n).  Since  each  A*  is  irreducible,  it  has 

a  positive  Perron  vector  Xi  corresponding  to  its  Perron  root  Aj.  If  Di  is  the  positive 
diagonal  matrix  obtained  by  placing  the  entries  of  Xi  in  order,  down  the  diagonal, 
then  D~^AiDi  has  all  its  row  sums  equal  to  Xi  using  Theorem  8.  We  may  assume 
that  Xi  >0.  If  A  is  initially  normalized  so  that  Ai  =  1  =  maxi^i^mXi,  and  An  is 
arbitrarily  chosen  to  be  x  m,  i  =  2, . . . ,  m,  with  nonnegative  entries  and  row  sums 
1  —  At,  then 

DfUiDi  0  ...  0 

Aai  D^^A2D2  ...  0 

A=  i  0  ; 

Ami  0  •  •  •  D^AmDm 

is  row  stochastic  and  satisfies  a  (i)  =  {A).  Thus  we  have  the  equivalence  of  the 
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row  stochastic  inverse  eigenvalue  problem  and  the  (normalized)  NIEP . 

While  the  row  stochastic  inverse  eigenvalue  problem  is  equivalent  to  the  NIEP ,  the 
doubly  stochastic  inverse  eigenvalue  problem,  is  not.  This  question  is  addressed  by 
Johnson.  Johnson  [12]  gave  necessary  and  sufficient  conditions  on  a  list  of  n  complex 
numbers  so  that  it  is  the  spectrum  of  a  doubly  stochastic  matrix.  Recall  a  doubly 
stochastic  matrix  is  a  square  matrix  that  has  row  sum  one  and  column  sum  one. 
The  following  two  observations  show  that  the  doubly  stochastic  inverse  eigenvalue 
problem  is  equivalent  to  the  row  stochastic  inverse  eigenvalue  problem  for  the  2x2 
case,  and  not  equivalent  for  n  ^  3. 

1.  Every  2x2  row  stochastic  matrix  is  similar  to  a  doubly  stochastic  matrix. 

Proof.  A  2  x  2  row  stochastic  matrix  has  the  general  form 

a  1  —  a 

P  l-P 

for  0  ^  a,  /3  <  1,  which  is  similar  to 

l-S+a  1-a+P 
2  2 

) 

l-Q+j3  1-^+a 
2  2 

a  doubly  stochastic  matrix.  ^ 

2.  There  exist  row  stochastic  matrices  that  are  not  similar  to  any  doubly  stochastic 

matrix  for  order  3  or  more. 
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Proof.  Consider  the  following  row  stochastic  matrix, 

0  1  0 

>1=  1  0  0  • 

1  0  0 

A  is  not  similar  to  any  doubly  stochastic  matrix.  For  if  one  did  exist,  it  would 
have  the  form 

a  b  1-a-b 

B  =  c  d  I- c- d 

1  —  o  —  c  1  —  b  —  d  a  +  b  +  c  +  d— 1 

So  for  A  and  B  to  be  similar,  both  the  traces  and  the  2x2  principal  minor 

sums  would  have  to  be  equal.  Thus, 

(a)  2{a  +  d)  +  {b  +  c)-l  =  0  and 

(b)  3{ad-bc)A{a  +  d)  +  2{b  +  c)-2  =  -l. 

But  since  the  (3,3)  entry  of  B  must  be  nonnegative,  this  means  a-hd  ^  1  — (fe+c). 
By  (a),  we  have  2{a  +  d)  =  l-{b  +  c)  resulting  in  the  conclusion  that  a  +  d  =  0 
implying  o  =  0  and  d  =  0.  This  means  that  6  +  c  =  1.  This  makes  (b)  become 
3  (—6c)  =  —36(1  —  6)  =  — 1  or  equivalently  6(1  —  6)  =  |.  This  has  no  solution 
for  0  <  6  <  1  showing  that  such  a  B  cannot  exist.  ^ 

For  more  information  regarding  the  doubly  stochastic  matrix  see  [12].  We  will  now 
return  to  the  NIEP  and  discuss  some  known  solutions. 
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3.5  2x2  and  3x3  NIEP  Solutions 

Two  cases  of  the  NIEP  are  completely  solved.  The  two  cases  are  when  o  —  {Ai,  A2} 
and  when  cr  =  {Ai,  A2,  A3}-.  We  will  begin  when  the  set  o  contains  two  numbers.  We 
observed  in  this  case  Ai  and  A2  must  be  real. 

The  list  cr  =  {Ai,  A2}  must  satisfy  Ai  +  A2  ^  0  where  Ai  ^  |A2l-  The  following 
(symmetric)  nonnegative  matrix  realizes  a: 

Ai  +  A2  Ai  —  A2 
Ai  —  A2  Ai  +  A2 

Now  in  the  case  when  a  =  (Ai,  A2,  A3},  Loewy  and  London  [20]  found  the  necessary 
and  sufficient  conditions.  They  are  given  in  Theorem  9. 

Theorem  9.  [20]  Let  a  =  {Ai,  A2,  A3}  he  a  list  of  three  complex  numbers,  and  assume 
that  a  satisfies  the  following  conditions 

1.  maxIlAd  :  \i  e  cr}  e  a, 

2.  a  =  a, 

3.  Si  =  Ai  +  A2  +  A3  ^  0,  and 

4.  si  ^  3S2. 

Then  a  is  realized  by  a  nonnegative  matrix  A. 

Proof.  If  a  is  real,  the  result  will  follow  from  Suleimanova’s  result  or  Perfect’s  result 
mentioned  earlier. 
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If  <7  is  complex,  then  a  =  {pe*^,pe~*®,r},  where  0  <  0  <  tt  and  0  <  p  ^  r.  As 
we  know  the  hypotheses  (l)-(4)  are  homogeneous,  we  may  assume  p  =  1-  This  gives 
a  =  {e*^,  r},  where  0  <  6  <tt  and  r  ^  1. 

Now,  the  trace  condition  Si  =  Ai  +  A2  +  A3  ^  0  and  <7  =  {e  ,  e  ,  r}  imply 
2  cos 0  +  7  ^  0;  and  substituting  into  sf  ^  3s2  the  information  a  —  {e*®,e 
implies  3  (2  cos  20  +  r^)  -  (2  cos  0  +  rf  ^  0. 

Thus,  (r  —  2  cos  (|7r  +  0))  (r  —  2  cos  (|7r  —  0))  ^0  follows. 

But  for  0  <  0  <  TT,  r  ^  1,  we  have  r  —  2  cos  (jTt  +  0)  ^  0  which  means  r  — 

2  cos  (|7r  —  0)  ^  0. 

So,  if  a  is  of  the  form  a  =  e"*®,  r},  then  our  last  two  hypotheses  become 


1.  2cos0  +  r  ^  0  and 


2.  r  -  2  cos  (Itt  -  0)  >  0  respectively. 

Now  construct  a  nonnegative  matrix  A  having  <7  =  {e*^,  r}  as  its  spectrum.  Let 
U  be  the  orthogonal  matrix 

I 

x/2  ^/3  -1 

\/2  0  2 
x/2  -73  -1 


U 


76 


and  let 


1  1  1 
1  1  1 
1  1  1 
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3  0  0 


Then  we  have, 

U^JU=  lo  0  0 
^0  0  0 

So  cr  =  {e*®,  r}  is  the  spectrum  of  the  real  matrix 

0  0 

A'  =  0  cos  6  sin  6 

^0  —  sin  9  cos  9 

Therefore,  by  the  matrix  multiplication 

A  =  UA'U'^  = 


=  U 


- 

“ 

r 

0 

0 

0 

0 

0 

0 

0 

0 

u'^  +  u 

0 

cos  9 

sin  9 

0 

0 

0 

0 

—  sin0 

cos  9 

COS 


—  cos  9  cos  (^TT  +  9)  cos  (ItT  —  9) 

(iyr-^)  -cos  9  cos(|7r  +  6') 
j^COS  (^TT  +  9)  cos  (Itt  —  9)  —  cos  9 

It  follows  that  A  is  nonnegative. 

The  following  geometrical  interpretation  of  Theorem  9  comes  from  the  work  of 
Dmitriev  and  Dynkin  [5,  6]. 


□ 
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Corollary  1.  [5,  6]  a  =  {r^  ^  pe  where  0  ^  p  ^  r  is  the  set  of  eigenvalues  of 

a  nonnegative  matrix  if  and  only  if  pe^^  belongs  to  the  closed  triangle  whose  vertices 

2Tri  _  2-7ri 

are  r^re  ^  ,re  3  . 

This  ends  the  solution  of  the  NIEP  for  any  lists  of  length  two  or  three.  The 
next  section  provides  the  solution  to  lists  of  length  four  or  five  when  the  sum  of  the 
numbers  in  the  list  is  zero.  This  is  called  the  trace  zero  NIEP. 

3.6  Trace  Zero  NIEP 

A  special  set  of  matrices  are  those  with  trace  zero.  If  A  ^  0  and  trace  (A)  =  0,  then 
it  follows  that  the  main  diagonal  entries  must  be  zero.  Reams  [29]  completely  solved 
the  trace  zero  NIEP  for  lists  of  length  four  as  we  see  in  the  following  Theorem  10. 

Theorem  10,  [29]  Let  a  ==  {Ai,  A2,  A3,  A4}  he  a  list  of  four  complex  numbers.  If 
Si  =  0,  S2  ^  0,  S3  ^  0,  and  As  a  ^  si,  then  there  exists  a  nonnegative  4x4  matrix 
with  spectrum  a. 

Proof.  Notice  that  the  last  inequality  among  the  hypotheses  of  the  theorem  is  the 
JLL  condition  with  n  =  4,  fe  =  2,  and  m  =  2. 

Let 

Pi  =  Ai  +  A2  +  A3  +  A4, 

P2  =  A1A2  +  A1A3  +  A1A4  +  A2A3  +  A2A4  +  A3A4, 

Pz  —  A1A2A3  +  A1A2A4  +  A1A3A4  +  A2A3A4,  and 
Pa  =  A1A2A3A4. 
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i.e.  Pi,P2,P3,P4  are  Newton’s  elementary  symmetric  polynomials. 
Newton’s  identities  for  symmetric  functions  [33]  state  that 

Si  -  Pi  =  0, 

S2  “■  PlSl  +  2p2  =  0, 

S3  -  P1S2  +  P2S1  -  3p3  =  0,  and 
S4  —  P1S3  +  P2S2  “  P3S1  +  4p4  =  0. 

We  can  write  these  equations  in  matrix  form  with  a  companion  matrix  as 


- 

- 

- 

" 

0 

1 

0 

0 

-4 

0 

0 

0 

-Si 

-3 

0 

0 

0 

0 

1 

0 

-Si 

-3 

0 

0 

-S2 

-Si 

-2 

0 

0 

0 

0 

1 

-S2 

-Si 

-2 

0 

-S3 

-S2 

-Si 

-1 

-P4 

P3 

-P2 

Pi 

-S3 

-S2 

-Si 

-1 

-S4 

-S3 

-S2 

-Si 

Letting  pi  =  Si  =  0  and  multiplying  both  sides  on  the  right  by 
diag  -l)  we  get 


- 

_ 

0  10  0 

10  0  0 

0  10  0 

0  0  10 

0  10  0 

^010 

4 

0  0  0  1 

f  0  1  0 

4 

r— 1 

0 

-P4  P3  -P2  0 

I— 1 

0 

1 

a  £2  0 

4  3  2 
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It  is  easily  verified  that 


- 

-*-1 

- 

“ 

1 

0 

0 

0 

1 

0 

0 

0 

0 

1 

0 

0 

0 

1 

0 

0 

4 

0 

1 

0 

4 

0 

1 

0 

4 

£2. 

3 

0 

1 

1 

1 

_£2. 

3 

0 

1 

then  multiplying  on  the  left  of  (1)  with  this  matrix  we  get 


1 

- 

- 

" 

1 

0 

0 

0 

0 

1 

0 

0 

1 

0 

0 

0 

0 

1 

0 

0 

0 

1 

0 

0 

0 

0 

1 

0 

0 

1 

0 

0 

S2 

4 

0 

1 

0 

S2 

4 

0 

1 

0 

0 

0 

0 

1 

4 

0 

1 

0 

4 

12 

0 

1 

4 

—  £2. 

3 

0 

1 

-P4 

P3 

-P2 

0 

4 

£1 

3 

0 

1 

354-5^ 

12 

12 

£2 

6 

0 

Finally,  performing  a  similarity  on  this  matrix  we  have 


p 

- 

- 

- 

" 

1 

0 

0 

0 

0 

1 

0 

0 

1 

0 

0 

0 

0 

1 

0  0 

0 

1 

0 

0 

§2. 

4 

0 

1 

0 

0 

1 

0 

0 

S2 

4 

0 

1  0 

0 

0 

1 

0 

£3. 

4 

S2 

12 

0 

1 

! 

0 

0 

1 

0 

£3. 

4 

0 

0  1 

0 

£2 

12 

0 

1 

354  — ^2 
12 

£3. 

12 

£2 

6 

0 

0 

_£2 

12 

0 

1 

4S4“S2 

16 

£2 

12 

- 1 

0 

which  is  nonnegative  and  similar  to  the  companion  matrix  with  eigenvalues  Ai,  A2,  A3, 
A4  as  required.  ^ 

Reams  gave  sufficient  conditions  for  the  trace  zero  NIEP  for  the  lists  of  five  num¬ 
bers.  The  5x5  trace  zero  NIEP  was  completely  solved  by  Laffey  and  Meehan  [18]. 
Since  the  proof  involves  results  from  combinatorics,  we  merely  state  the  results. 


25 


Theorem  11.  [18]  Let  Ai,  A2,  A3,  A4,  A5  be  complex  numbers  and  let  Sk  =  Ai*4-A2*^4- 
Aa*^  +  A4*  +  As*^  for  k  =  1,  2,  ...  .  Assume  Si  =  0.  Then  a  —  {Ai,  A2,  A3,  A4,  A5} 
is  the  spectrum  of  a  nonnegative  5x5  matrix  if  and  only  if  the  following  conditions 
hold: 


1-  Sfc  ^  0  for  k  =  1,2,  3,  4,  5, 

2.  4s4  ^  s\,  and 

3.  12s5  -  5S2S3  +  5s3-y/4s4  -  si  ^  0. 

A  refinement  of  the  JLL  condition  in  the  case  of  trace  zero  lists  of  odd  length 
comes  from  Laffey  and  Meehan  [17].  Laffey  and  Meehan’s  result  above  uses  this 
refinement  in  condition  (2). 

Theorem  12.  [17]  Let  A  be  a  nonnegative  n  x  n  matrix  with  trace  {A)  =  0.  Then, 
if  n  is  odd,  (n  —  1)  S4  ^  si,  that  is,  {n  —  1)  trace  (A^)  ^  {trace  {A?yf . 


This  inequality  is  the  best  possible.  Consider  the  matrix  for  odd  n, 


A  = 


0  1 


©•••© 


1  0 


0  1 
1  0 


©  [0] , 


then  equality  occurs.  For  n  even,  the  matrix 


- 

0  1 

0  1 

A  = 

©•••© 

1  0 

1  0 

gives  an  example  that  shows  Theorem  12  fails  when  n  is  even. 
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This  concludes  the  known  results  regarding  the  NIEP  for  lists  of  length  n  ^  5. 
In  the  next  chapter  will  return  to  the  list  of  four  complex  numbers  in  an  attempt  to 
solve  the  general  NIEP. 


27 


CHAPTER  4 


Lists  of  Length  Four 

In  this  chapter,  we  will  examine  a  list  of  four  complex  numbers.  We  first  review 
what  has  been  solved  in  the  smaller  dimensions. 

The  1x1  case  is  trivial.  Clearly  the  eigenvalue  must  be  nonnegative  to  satisfy 
the  Perron  root  condition  and  so  1  x  1  nonnegative  matrix  is  the  Perron  root.  Now, 
when  the  list  contains  two  elements,  the  eigenvalues  must  be  real.  This  is  because 
the  characteristic  polynomial  of  a  nonnegative  matrix  has  real  coefficients,  so  the 
eigenvalues  must  occur  in  conjugate  pairs.  As  one  must  be  the  Perron  root,  the  other 
cannot  be  a  complex  number  with  nonzero  imaginary  part.  We  have  seen  the  solution 
to  this  case  in  Chapter  3. 

Now  in  the  case  where  the  list  has  three  elements,  we  have  seen  in  Chapter  3  that 
Theorem  9  by  Loewy  and  London  solves  the  NIEP.  The  conditions  from  Theorem  9 
are  not  sufficient  for  the  list  of  four  complex  numbers.  To  see  this,  consider  o  = 
{\/2,  \/2,2,  -i).  As  this  set  satisfies  all  the  necessary  conditions,  it  appears  to  be  a 
good  candidate  for  the  spectrum  of  a  4  x  4  nonnegative  matrix.  But  according  to 
the  Perron-Frobenius  Theorem,  if  a  is  realizable  by  a  nonnegative  matrix  then  it  is 
reducible.  The  only  possible  partition  requires  the  grouping  a  —  {z,  — i,\/2}.  But 
this  d  fails  the  necessary  requirements  to  be  a  spectrum  of  a  3  x  3  nonnegative  matrix 
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(as  seen  in  Chapter  3),  so  a  cannot  be  the  spectrum  of  a  4  x  4  nonnegative  matrix. 

Thus  in  search  of  sufficient  conditions  for  a  solution  to  the  NIEP  for  lists  of  four 
numbers,  we  divide  the  problem  into  two  parts.  One  is  when  a  is  real,  and  the  other 
is  when  o  has  two  complex  numbers  which  are  complex  conjugates  of  each  other.  The 
first  one  has  a  complete  solution  given  by  Loewy  and  London  which  will  be  presented 
in  Chapter  6.  But  for  the  second  case,  when  o  —  {Ai,  A2,  a  +  bi,  a  —  where  Ai  is 
the  Perron  root,  A2  is  real,  and  6  0,  no  solution  is  found  in  the  literature.  We  will 

begin  with  the  case  where  we  have  normalized  the  spectrum  by  dividing  by  Ai  and 
setting  a  =  0. 

4.1  NIEP  for  (7  =  {1,  A,  bi,  —bi} 

We  will  analyze  the  list  <7  =  {l,A,6i,— 6*}.  If  6^0,  we  know  that  nonnegative 
matrices  with  a  as  spectrum  cannot  be  symmetric  since  symmetric  matrices  have  real 
spectra. 

Analysis  on  the  A  and  b  helps  us  limit  the  region  under  consideration.  Clearly,  as 
1  is  the  Perron  root,  then  both  —1  ^  A  ^  1  and  —1  ^6^1  must  be  true.  Moreover, 
as  we  have  both  bi  and  —bi,  we  may  assume  0  <  6  <  1.  So  when  we  graph  the  region, 
we  have  a  rectangular  region  in  terms  of  A  and  b  where  —  1  ^  A  ^  1  and  0  <  ^  1. 

Since  <t  must  satisfy  the  the  JLL  condition,  the  region  for  possible  A  and  b  becomes 
smaller.  The  necessary  condition  4s2  ^  s?  gives  4  (1  +  A^  —  26^)  ^  (1  +  ^)  •  This 
in  turn  restricts  the  region  to  0  <  6  ^  and  —  1  ^  A  ^  1.  No  other  JLL 

inequality  restricts  the  region  further. 
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Figure  1:  Regions  for  {1,  A,  bi,  —bi} 

We  will  denote  the  regions  possible  by  the  following: 

1.  Region  I  will  be  for  the  area  defined  when  —  1  ^  A  ^  0  and  ^  ^  6  ^ 

o 

\/s  —  2A  +  3), 

2.  Region  II  will  cover  the  area  for  -1  <  A  ^  0  and  0  <  6  ^ 

3.  Region  III  has  area  0  ^  A  ^  1  and  0  <  6  <  and 

4.  Region  IV  covers  the  area  defined  by  0  ^  A  ^  1  and  ^  <  6  ^  (SA^  -  2A  +  3). 

We  prove  that  for  (A,  b)  in  regions  I,  II,  or  III  there  exist  nonnegative  matrices 
whose  spectrum  is  {1,  X,bi  —  bi}.  However,  region  IV  is  still  unknown  at  the  time  of 
this  submission. 

4.1.1  Solution  of  <7  =  {1,  A,  bi,  —bi}  when  A  =  — 1  and  0  <  6  ^  1 

A  solution  to  this  case  which  realizes  regions  I  and  II  is  relatively  simple.  Consider 
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the  following  matrix: 


0  0  1  +  b  1-b 

0  0  1-6  1  +  6 
1-6  1+6  0  0 
1+6  1-6  0  0 
Here  cr(A)  =  {1,  — 1, 6i,  — 6i}  for  0  <  6  ^  1.  We  prove  this  construction  by  using 
Theorem  6. 

Theorem  13.  Let  a  =  {1,  —l,ib,  —ib}  be  such  that  0  ^  6  ^  1.  Then  a  is  the  spectrum 
of  the  nonnegative  matrix  A  where 

0  0  1  +  6  1-6 
0  0  1-6  1  +  6 
1-6  1+6  0  0 
1+6  1-6  0  0 


Proof.  Let 

and 


^1.2  =  ;; 


^2,1 


1  +  6  1-6 
1-6  1+6 

1-6  1  +  6 
1  +  6  1-6 
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Then 


1  +  6  1  —  6  1  —  6  1  +  6 

—  T 

1-6  1  +  6  1  +  6  1-6 
2(l  +  6)(l-6)  (l  +  6)'  +  (l-6)' 

{1  +  bf  +  il-bf  2(l  +  6)(l-6) 

2  (1  -  6^)  2  +  26^ 

2  +  262  2  (1  -  62) 

1-62  1  +  ^,2 
1  +  62 

The  characteristic  polynomial  of  Ai,2-A2,i  is 

5-  1  (1  -  62)  -1  (1  +  62) 

-i(l+62)  5-1(1-62) 

=  52  -(1  -  62)  5  -  62 

=  (•5-1)  (S  +  b^). 

Now,  setting  the  characteristic  polynomial  equal  to  zero  gives  5  =  1  or  5  =  —62.  But 
by  Theorem  6,  we  know  the  eigenvalues  of  the  matrix  A  are  just  the  square  roots  of 
5.  This  shows  that  the  spectrum  of  A  is  cr  =  {1,  —1, 6z,  — 6i}.  D 
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4.1.2  Solution  of  cr  =  {1,  A,  bi,  —bi}  when  0  <  b  ^  ^  and  0  ^  A  ^  1 

In  this  case  (A,  b)  belongs  to  region  III.  An  example  of  a  nonnegative  matrix  in 
this  region  is: 


where  a  (A)  =  {1,  ^i,  -^0- 

Note  that  in  general,  the  list  {l,bi,—bi}  for  0  <  6  ^  ^  satisfies  Corollary  1. 
Hence  there  exists  Ai  ^  0  with  spectrum  {l,bi,  -bi}.  Since  A  ^  0  the  nonnegative 
matrix  A  =  [A]  ©  Ai  is 

AO  0  0 

0  I 

A  ^ 

0  1(1  -  \/36)  I  3 

0  I  (l  +  V^b)  I  (l  —  V^b)  I 

which  realizes  a. 

Therefore,  a  solution  exists  to  the  NIEP  for  a  =  (1,  A,  bi,  —bi}  with  Perron  root  1 
when  Q  <  b  ^  ^  and  0  ^  A  ^  1. 

4.1.3  Solution  of  cr  =  {1,  A,  bi,  —bi}  when  0  <  6  ^  ^  and  — 1  <  A  ^  0 

Here  we  examine  region  II.  We  saw  in  Chapter  3  that  Reams  [29]  completely  solved 
the  NIEP  for  a  list  of  four  numbers  that  have  trace  zero.  The  sufficient  conditions 


33 


are  Si  =  0,  S2  ^  0,  S3  ^  0,  and  4s4  ^  sj.  The  solution  due  to  Reams  is 

0  10  0 

T  0  1  0 

"  .  (2) 

^001 
4 

4S4  — 52  SS  ^  f\ 

16  12  4 

We  will  use  this  result  to  find  the  solution  for  region  II.  The  solution  is  given  in 
Theorem  14. 

4.1.4  Solution  of  a  =  {1,  X,bi,—bi}  when  ^  ^  ^  ^  \J  I  —  2A  +  3)  and 
-1  <  A  <  0 

Here  we  examine  region  I.  Recall  the  upper  limit  for  b  in  region  I  comes  from  the 
JLL  inequality 

^  4s2  (1  +  A)^  ^  4  (1  +  A^  —  2b^) 

•4=^  1  4-  2A  4-  A^  ^4  4"  4A^  —  86^ 

44  0  ^  3A^  —  2A  +  (3  —  86^)  , 

which  is  a  quadratic  inequality  in  b  and  A  where  the  roots  of  0  =  3A  —  2A  +  (3  —  86  ) 
in  terms  of  6  are  A  =  ^  (2  ±  ^9662- 32).  Now  since  A  is  a  real  number  it  follows 
that  966^  -  32  ^  0  44  966^  >  32  44  62^|446^^.  Moreover,  A  =  ^  is  a  double 

■v/T 

root  to  the  quadratic  equation  when  b  =  —. 

Now  b  =  yjl  (3A2  -  2A  4-  3)  is  a  root  of  the  equation  6^  =  |  (SA^  -  2A  +  3)  from 
the  JLL  inequality.  But  b  is  assumed  to  be  a  real  number;  so  to  verify  this,  we  observe 
that  3A2  -2A+3  >  0  for  all  real  A  as  the  roots  of  this  quadratic  are  complex  numbers. 
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For  -1  <  A  ^  0  and  f  ^  ^  yjl  (SA^  -  2A  +  3)  there  exists  a  nonnegative 

matrix  realizing  tr  -  {1,  A,  bi,  -bi}  as  its  spectrum  and  we  prove  this  in  Theorem  14. 

4.1.5  Solution  of  a  =  {l,X,bi,-bi}  when  ^  ^  b  ^  (SA^  -2A  +  3)  and 
0<  A^  1 

This  is  region  IV.  Recall  the  upper  limit  for  b  in  region  IV  comes  from  the  JLL 
inequality.  For  0  ^  A  ^  |  and  ^  ^  5  ^  (3A2-2A  +  3)  Theorem  14  provides  a 

nonnegative  matrix  that  realizes  a.  The  remaining  part  of  region  IV  when  i  <  A  ^  1 
is  still  unresolved  at  the  time  of  this  submission. 

4.1.6  Results 

The  following  theorem  gives  a  nonnegative  matrix  realizing  a  =  {l,X,bi,-bi} 
when  (A,  6)  belongs  to  regions  I,  II  and  parts  of  IV.  Recall  we  completely  solved 
region  III  in  section  (4.1.2).  A  part  of  this  result  is  found  in  Reams  [28]. 

Theorem  14.  Let  a  =  {1,  A,  bi,  -bi}  where  -1  <  A  <  0  and 

0<b^  y/l  (3A2  -  2A  +  3)  orO^X^I  and^^b^  ^l(3X^-2X  +  3).  Suppose 
a  satisfies  all  the  necessary  conditions  of  the  NIEP.  Then  a  nonnegative  matrix  B 
which  realizes  o  as  its  spectrum  is  A where  A  is  Reams  ’  solution  (2)  for  the 
list  {1  -  X  -  ^,bi  -  -bi  -  and  I  is  the  identity  matrix. 

Proof.  Suppose  {1,  A,  bi,  -bi}  can  be  realized  by  a  nonnegative  matrix  B  with  con¬ 
stant  diagonal  entries  then  B  -  has  trace  zero.  A  matrix  with  eigenvalues  of 
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B  -  namely  -bi  -  ^},  is  given  by  Reams  as 

0  10  0 

i  0  10 

4 

? 

4  0  0  1 

4 

4.  £i  n 

16  12  4  ^ 

where  s-  are  the  moments  corresponding  to  the  eigenvalues  of  R  -  If  we  show 

A  ^  0  then  it  follows  A  +  ^  0  and  the  spectrum  of  A  +  is  the  list 

{1,  X,bi,—bi}.  We  will  analyze  the  entries  of  A. 

Throughout,  Sj  will  denote  the  moment  of  the  eigenvalues  of  B,  namely  cr  = 
{l,X,bi,-bi}.  Note  that  we  replace  1  +  A  with  Si  in  the  following  calculations.  To 
show  the  (2, 1)  and  (4, 3)  entries  are  nonnegative  we  need  only  show  that  Sj  ^  0. 
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By  assumption  a  satisfies  the  JLL  conditions,  so  we  have  4s2  ^  Thus  4s2  s?  ^  0 
from  which  follows  the  desired  conclusion  oi  $2  ^  0  for  —  1  <  A  ^  |  and  0  <  6  ^ 
y'l(3A2-2A  +  3). 

To  show  the  (3, 1)  and  (4, 2)  entries  are  nonnegative  we  need  to  show  that  Sg  ^  0. 
Now, 

s'g  =  trace 

=  trace  ^  ~  ^ 

=  trace  -  3  +  3  B  -  /  j 

8S3  —  6SiS2  +  sf 
^  8  ■ 

Reams  [29]  shows  that  if  cr  =  {r,  A2,  A3,  A4}  where  r  is  the  Perron  root,  then 
8S3  —  6siS2  +  5^  =  3  (r  +  A2  —  A3  —  A4)  (r  —  A2  +  A3  —  A4)  (r  —  A2  —  A3  +  A4).  So  in 
our  case  we  have  r  =  1,  A2  =  A,  A3  =  6i,  and  A4  =  —hi.  With  substitution  and 
recalling  1  +  A  ^  0  from  the  Perron  root  condition,  we  have 

8S3  -  6S1S2  +  =  3  (1  +  A)  (1  -  A  +  2hi)  (1  -  A  -  2hi) 

=  3  (1  +  A)  [(1  -  A)^  +  46^] 

^  0. 

Thus  S3  ^  0.  Now  to  show  the  (4, 1)  entry  is  nonnegative  we  will  need  the  equivalent 
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of  S4  in  terms  of  the  moments  of  B. 


S4  =  trace 
=  trace 


trace  \  ^ 


'l  +  A' 


5^+6 


B^ 


1  +  A' 


5  + 


'1  +  A' 


3  2 

=  S4  -  S,S3  +  jS,S2  -  jg  +  54 

3  3 

=  S4  —  S1S3 + ^5452  ~ 


Therefore 


44  -  (4)  =^{si-  S1S3  +  ^Si52  -  ^4)  -  (^4  (4^2  4) 

16S4  —  16Si53  +  8S4S2  ~  4^2  ~ 

4 

Rewriting  the  numerator  using  Si,  S2,  S3,  and  S4  in  terms  of  1,  6,  and  A,  we  have 


I6S4  —  I6S1S3  +  8S4S2  ~  4s2  ■“  s| 

=  16  (1  +  A^  +  2b^)  -  16  (1  +  A)  (1  +  A^)  +  8  (1  +  A)^  (1  +  A^  -  26^)  -  (1  +  A)^ 

-  4  (1  +  A2  -  2l^f 

=  3  +  I6fe4  _  4A  -  3262a  +  2A2  -  +  3A^ 

=  2A2  +  3A^  -  a  (4  +  4A2  +  3262)  4.  (3  +  156") 

which  has  positive  terms  with  the  only  possible  exception  of  -A  (4  +  4A2  +  3262). 
But  this  term  is  nonnegative  when  -1  <  A  ^  0.  Therefore,  for  -1  <  A  ^  0  and 
0  <  6  ^  (3A2  -  2A  +  3)  we  have  4^4  -  (52)^  ^  0.  Thus  B  is  nonnegative. 

Assume  0  ^  A  ^  |  and  ^  ^  6  ^  {3A2-2A  +  3).  Consider  from  above,  the 

polynomial  P  (A,  6)  =  3A^  -  4X^  +  2X^-4  (862  +  1)  a  +  (166^  +  3)  over  0  ^  A  ^  | 
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and  ^  <  6  ^  a  rectangle  which  contains  the  region  0  ^  A  ^  3  and  3  ^  ^ 

y^l  (3A2  —  2A  +  3).  We  will  show  P  (A,  6)  ^  0  in  this  rectangle. 

Now  ■ 

—  -  12A^  -  12A^  +  4A  -  4  (86^  +  l) 
dX 

=  4[3A^-3A2  +  A-  (862  +  1)] 

^  4  ^3A^  —  3A2  +  a  — 
using  6  =  ^  and  gf  =  36A2  -  24A  +  4  =  4  (QA^  -  6A  +  1)  =  4  (3A  -  l)^  so  f  <  0 
on[0,|]x[f,y^andgj^0. 

Hence  P  (A,  6)  is  decreasing  and  concave  up  for  (A,  6)  in  the  rectangle.  ^  —  0 
when  A  =  I  for  6  in  •  Note  that  (A,  6)  =  ~  (^’  V^) 

points  on  the  parabola  coming  from  the  JLL  condition. 

Now,  P{lb)  =  166‘  -  +  ?  with  P(l,f  )  =  0  and  P(l.-/i)  =  i- 

Consider  f  (1.6)  =646(i,-^)  (i+f)-  This  shows  f  (i, 6)  >  0  if  f  «  I- < 

Clearly  P  (0, 6)  ^  0  for  ^  <  b  <  ^1- 

Since  P  (A,  6)  is  decreasing  and  concave  up  with  respect  to  A  on  [0,  jlx 
and  P(0,6)  >  0,  P(i,6)  ^  0  for  f  <  6  <  \/|,  it  follows  ‘hat  P(A.i.)  ?  0  tor 
(A.b)  in  [0,flx[^,\/|]-  Since  this  contains  the  region  0  <  A  «  A  and  ^  «  6  « 
(3A2  -  2A  +  3)  it  follows  that  4s'4  ^  (s'2)^. 

Thus  A  is  nonnegative  as  claimed.  ^ 
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As  an  example,  the  following  nonnegative  matrix  realizes  {1, 6>  y  32  y  32 

1  0  0 

s  1  » 

0  H  1 

_  1043  ri  J_ 

|_  110592  13824  24_ 

A  final  observation  occurs  when  we  attempt  to  use  the  following  result  in  the 
unknown  case.  Here  we  can  construct  various  real  matrices  having  the  spectrum 
{1,  A,  M,  -bi}.  This  method  did  not  always  provide  nonnegative  matrices. 

Theorem  15.  [2]  {1,  A2,  •  ■  • ,  A„}  is  realizable  by  a  nonnegative  matrix  if  and  only  if 
there  exists  some  real  matrix  B  of  order  n-l  with  spectrum  {A2,  •  •  ■ ,  A„}  and  some 
[n  —  l)-simplex  S  C  with  6  S  such  that  BS  C  S. 

In  our  case  {1,  A,  bi,  -bi}  is  realizable  by  a  nonnegative  matrix  if  and  only  if  there 
exists  a  3  X  3  real  matrix  B  with  spectrum  {A,  bi,  -bi]  and  a  tetrahedron  S  containing 

fjp 

(0, 0, 0)  such  that  point  p  =  (pi,P2,P3)  e  S  B  [pi,P2,P3]  e  S. 

There  are  many  choices  for  B  with  spectrum  {A,  bi,  -bi}.  One  such  matrix  is 

0-6  0 
6  0  0- 

0  0  A 

Now  consider  the  tetrahedron  with  vertices  Vi  -  {xi,yi,Zi)  for  i  =  1,2, 3, 4.  The 


which  is  in  the  left  portion  of  region  IV 

7_ 

24 

0 

1043 

4608 

4199 
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centroid  of  the  tetrahedron  with  vertices  vi,V2,V5,  V4  is 

centroid  =  7  (fxi  +  0:2  +  +  X4) ,  (t/i  +  2/2  +  ys  +  Va)  ,  (2:1  +  22  +  23  +  ^4))  • 

4 

This  centroid  divides  the  line  segment  joining  a  vertex  to  the  centroid  of  the  opposite 
face  in  the  ratio  3:1.  The  centroid  of  the  face  triangle  is  the  intersection  of  all  its 
medians. 

Select  va  =  (0, 0, 3).  We  desire  the  tetrahedron  to  be  oriented  so  that  the  centroid 
is  located  at  the  origin.  Then  the  other  end  of  the  line  segment  joining  Va  to  the 
centroid  is  (0,0,  -1).  Note  this  point  is  the  centroid  of  the  triangle  formed  by  Vi,  V2, 
and  H3. 

Let  the  vertices  be  v\  =  (a:i,0,2;i),  V2  —  {x2,y2,  ^2)^  ^3  ==  {^3,yz->^3)i  S'lid  ^4  — 
(0,0,3).  Thus  the  centroid  of  this  tetrahedron  is 

centroid  =  7  {(xi  +  X2  +  ^z) ,  (^2  +  yz)  >  i^i  +  22  +  23  +  3)) . 

4 

This  gives  V3  =  (-Xi  —  X2,  -yz,  -  22  -  3). 

To  use  Theorem  15  we  consider  the  real  matrix 

0-5  0 
B=  6  0  0 

0  0  A 

where  we  are  rotating  a  tetrahedron  by  90  degrees  and  shrinking  it  by  6  relative  to 
the  ay-plane  and  then  compressing  it  by  A  along  the  z-axis.  We  want  to  know  if  every 
point  in  the  tetrahedron  when  acted  on  by  B  will  remain  in  the  original  tetrahedron. 
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As  a  tetrahedron  is  a  convex  set,  any  point  p  =  {pi,P2^P3)  the  tetrahedron  can 
be  described  by  the  convex  sum  p  =  Yli-iUVi  where  ^  0  for  i  =  1,2, 3, 4  and 

We  can  then  test  to  see  if  there  exists  a  nonnegative  matrix  associated  with  a  = 
{1,  A,  —bi,  bi}  by  solving  the  following  system  of  linear  equations  for  the  coefficients 


r  1 

1 

CO 

h 

1 

to 

Vl  y2  2/3  Di 

h 

bpi 

Zi  Z2  23  Za 

h 

AP3 

1111 

U 

1 

Here  p  =  {pi,P2,P3)  belonged  to  the  original  tetrahedron. 
Now  the  matrix 

Xi  X2  X3  Xi 

Vl  1)2  1)3  Di 

A  = 

Zi  Z2  Zs  Z4 
1111 


is  nonsingular  by  construction  as  the  vertices  are  noncoplanar;  therefore,  we  may  find 
its  inverse  to  solve  for  the  coefficients. 

Since  the  extreme  points  of  the  tetrahedron  are  its  vertices,  if  the  vertices  when 
acted  upon  by  B  remain  part  of  the  original  tetrahedron,  then  every  point  of  the  tetra¬ 
hedron  when  acted  on  by  B  will  remain  part  of  the  original  tetrahedron.  Therefore 
we  need  to  test  those  four  vertices,  Vi  for  i  =  l,2, 3, 4,  by  solving  for  their  correspond- 
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ing  coefficients.  If  these  coefficients  remain  nonnegative  then  the  matrix  C  whose 
columns  are  the  coefficients  corresponding  to  for  i  =  1, 2, 3, 4  will  be  a  nonnegative 
matrix  realizing  a  using  Theorem  15. 

Consider  the  following  example.  Let  vi  =  (4, 0,-4)  and  V2  =  (— 1,5,  — 1).  Then 
yj  =  (_3^_5j2)  while  V4  =  (0,0,3).  Solving  this  problem  with  Mathematica  when 
A  =  I  and  b  =  0.6,  we  find  the  matrix  C  is  approximately 


0.3622 

0.0206 

0.5339 

0.0833 

0.7222 

0.6806 

-0.4861 

0.0833 

0.2422 

0.8001 

-0.1261 

0.0833 

-0.3267 

-0.5017 

-0.1261 

0.75 

Notice  this  attempt  fails  to  take  vi,V2,  V3  into  the  original  tetrahedron  as  each  convex 
sum  coefficient  fails  to  be  nonnegative,  but  the  matrix  C  does  indeed  have  spectrum 
{l,|,0.6i,-0.6z}. 

4.2  NIEP  for  a  =  {1,  A,  a  4-  bi,  a  —  bi} 

Let  us  now  consider  the  following  four  cases  relative  to  the  complete  list  of  four 
numbers  where  a  ^  0. 

1.  Consider  A  =  1  and  a?  +  b'^  <  1.  Since  the  Perron  root  is  repeated  any  non¬ 
negative  matrix  realizing  <j  must  be  reducible.  By  the  sufficient  conditions  for 
{1,0  -t-  bi,a  -  bi}  to  be  realizable  by  a  nonnegative  matrix  from  Theorem  9. 
They  are 

(a)  1  -h  2a  ^  0,  the  trace  is  nonnegative,  and 
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(b)  (1  +  2af  <  3  (1  +  2a^  -  26^),  from  the  JLL  condition. 


If  (a)  and  (b)  are  satisfied,  then  by  Theorem  9  there  exists  a  nonnegative  matrix, 
a' ,  realizing  {l,a  +  bi,a  -  bi}.  Using  this  matrix,  we  obtain  the  following 
reducible  nonnegative  matrix  having  spectrum  a . 


A  = 


A'  0 


0  1 


2.  Suppose  A  =  — 1  and  +  b"^  <  1.  Then  a  =  {1,  — l,a  +  bi,a  —  bi},  and  for 
the  trace  condition  to  hold  we  need  o  ^  0.  Should  a  nonnegative  matrix  exist 
having  such  a  spectrum,  it  could  not  be  reducible  as  no  partition  would  work 
with  the  A  =  -1.  Therefore,  the  only  possible  nonnegative  matrix  must  be 
irreducible. 


Since  |-1|  =  1  we  have  S  =  {1,  -1}  in  Theorem  3.  This  means  a  +  bi  and  a-bi 
are  on  the  imaginary  axis  giving  o  =  0.  Since  <  1  we  get  —  1  <  6  <  1. 

But  we  have  seen  solutions  to  {1,  —1,  bi,  —bi}  for  0  <  6  <  1  in  subsection  (4.1.1). 


3.  A  =  -1  and  +  b‘^  =  1.  Then  we  have  four  eigenvalues  of  the  same  modulus 
and  they  must  occur  at  1,  —1,  i,  —i-  Therefore,  the  only  possible  spectrum 
with  these  conditions  that  have  a  nonnegative  matrix  associated  with  it  is  cr 
{1,  -l,i,  -i}.  We  have  seen  such  a  matrix  in  subsection  (4.1.1). 

4.  Suppose  A  =  1  and  +  62  ^  1.  By  the  Perron-FVoebenius  Theorem  we  know  a 
repeating  Perron  root  requires  the  nonnegative  matrix  with  spectrum  cr  to  be 
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reducible.  Therefore  the  only  possibility  is  when  a±bi  is  in  the  closed  triangular 
region  with  vertices  1, 6^,6"^  and  thus  {1, a  +  bi,a-  bi}  is  the  spectrum  of 
some  nonnegative  3x3  matrix  from  Corollary  1,  call  it  B.  Then 

B  0 

y4  = 

0  1 

realizes  a  =  {1, 1,  a  +  bi,  a  —  bi}. 

5.  Suppose  1A|  <  1  and  a'^  +  b^  =  1.  If  we  consider  an  irreducible  nonnegative 
matrix  then  we  must  abide  by  Theorem  3  which  says  that  if  we  have  three 
values  with  modulus  one,  the  remaining  value  A  must  have  two  other  distinct 
numbers  in  a  with  the  same  modulus  as  A.  That  is  not  possible.  When  we 
consider  the  reducible  nonnegative  matrix,  then  a  ±  bi  must  be  in  the  closed 
triangular  region  with  vertices  1,  e^,  and  A  ^  0.  In  this  case.  Corollary  1 
guarantees  B  ^  0  with  spectrum  {l,(i  +  bi,(i  —  bi}.  Then 


has  spectrum  a. 

6.  Suppose  |A1  <  1  and  a"  +  6'  <  1.  This  case  is  still  unresolved. 

As  a  final  result  of  this  chapter  we  make  the  following  observation. 

The  example  below  shows  that  it  is  not  possible  to  solve  the  4  x  4  NIEP  by 
considering  only  matrices  whose  main  diagonal  elements  are  equal  as  m  Theorem  14. 
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Suppose  a  =  {r,  A,  a  +  bi,  a  —  bi}  has  such  a  solution,  then  Reams’  expression  for 
^3  =  3  (^  +  -^2  ~  -^3  ■“  A4)  ~  A2  +  A3  —  A4)  (r  —  A2  —  A3  +  A4)  gives 

S3  =  3  (r  +  A  -  2a)  (r  -  A  +  2bi)  (r  -  A  -  2bi)  ^  0.  This  implies  (r  +  A  -  2a)  >  0. 
But  for  the  following  matrix  from  [28] 

0  6  0  0 

1  ii  2  ,  _i_  2 _ L 

^  3  3  3  v/3 

1  2 _ L  n  2  ,  J_ 

■‘■3  y/3  3  3  ~ 

1  2  ,  J_  2  _  _1_  U 

^  3  v/3  3  v/3  3 

the  spectrum  is  a  (A)  =  {6,  — 1, 3  +  3  —  i}.  Here  (r  +  A  —  2a)  =  6  —  1  —  6  =  — 1  <0. 
Hence  {6,  —1, 3  +  3  —  i}  cannot  be  realized  by  a  matrix  B  ^  0  with  equal  entries 

along  the  main  diagonal. 

4.3  Conclusion 

In  this  chapter  we  have  identified  nonnegative  matrices  realizing  a  list  of  four 
complex  numbers.  In  the  case  where  a  =  {1,  A,  bi,  —bi},  the  main  result  is  Theorem  14 
which  summarizes  the  known  regions.  A  part  of  region  IV  is  still  unresolved.  In  the 
general  case  we  addressed  the  boundary  point  possibilities  leaving  open  the  case  when 
<7  =  {1,  A,  a  4-  bi,  a  —  bi}  where  [A]  <  1  and  a^  +  6^  <  1. 

We  also  have  a  construction  technique  to  find  matrices  having  certain  lists  as  spec- 
trums  using  Theorem  15  and  Mathematica.  Although  this  process  was  not  fruitful 
in  providing  a  nonnegative  matrix,  it  gave  the  opportunity  to  test  many  different  real 
matrices. 
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CHAPTER  5 


Real  versus  Symmetric  Nonnegative  Inverse  Eigenvalue  Problem 

The  general  NIEP  is  open.  Some  progress  has  been  made  in  the  problem  of  the 
existence  of  a  nonnegative  matrix  with  a  prescribed  real  spectrum.  Chapter  5  focuses 
on  lists  of  real  numbers. 

5.1  Real  Nonnegative  Inverse  Eigenvalue  Problem  (RNIEP ) 

The  RNIEP  is  a  subproblem  of  the  NIEP  with  one  additional  constraint.  The 
RNIEP  restricts  the  elements  of  cr  =  {Ai,  A2,  •  •  • ,  A„}  to  be  from  the  set  of  real 
numbers.  The  RNIEP  is  to  determine  the  necessary  and  sufficient  conditions  on 
a  =  {Ai,  A2,  . . . ,  An}  C  R  so  that  there  exists  a  nonnegative  matrix  having  a  as  its 

spectrum. 

As  we  have  seen  in  Chapter  3,  Suleimanova  [32]  gave  a  sufficient  condition  for  the 
existence  of  a  nonnegative  matrix  when  the  list  contains  real  numbers.  The  simple 
proof  given  here  for  Suleimanova’s  Theorem  is  due  to  Perfect  [24]. 

Theorem  16.  [22]  Let  a  =  {Xi,  \2,  ■■■ ,  An}  be  a  list  of  n  real  numbers  satisfying 

Ai  "t“  A2  "h  ■  ■  ’  4"  An  ^  0 
Xj  <  0,  for  j  =  2,  3, . . . ,  n. 
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Then,  there  exists  a  nonnegative  n  x  n  matrix  with  spectrum  a . 


Proof.  Let 

/(A)  =  n(A-Ai) 

j=l  j  =  l 

We  assert  that  all  02,03,... c„  must  be  positive.  For  the  polynomial  /  (—A)  = 

A"  +  ^  (-1)^”^  Cj>P~^  has  exactly  (n  -  1)  positive  roots.  Therefore,  by 
i=i 

Descartes’  rule  of  signs,  the  number  of  variations  of  sign  in  the  sequence 

1,  Oi,  — 02,  C3,  — C4,  .  .  .  ,  (  — 1)  Cji 

is  exactly  (n  —  1).  However,  Oi  =  Ai  +  A2  +  •  •  •  +  A^  is  nonnegative,  and  it  follows 
that  all  the  other  Cj  must  be  positive.  Hence  the  companion  matrix  of  /  (A), 

0  10  0  0 

0  0  1  0  0 

0  0  1 

Ofi  Oji—\  Oji—2  •  •  •  ^2  Cl 

is  nonnegative  and  has  eigenvalues  Ai,  A2, . .  • ,  A^. 


The  following  corollaries  are  due  to  Perfect  [24]. 

Corollary  2.  [22]  Let  a  =  {Ai,  Aa, . . . ,  A^}  be  a  list  of  n  real  numbers  where 

Ai  ^  Aa  ^  Ap  ^  0  >  Ap+i  ^  An  and 

Ai  +  Ap+i  +  •  •  •  +  An  ^  0. 

Proof.  Use  the  proof  of  Theorem  16  to  construct  a  nonnegative 

(n  —  p  +  1)  X  {n  —  p  +  1)  matrix  B  with  eigenvalues  Ai,  Ap+i,  Ap+a, . . . ,  An-  Then  the 

nonnegative  matrix  B  0  diag  (Aa,  •  -  • ,  Ap)  has  spectrum  a.  n 

Corollary  3.  [22]  Let  a  be  a  list  of  n  real  numbers,  and  suppose  that  it  is  possible 
to  partition  o  into  sublists  in  such  a  way  that  (i)  each  sublist  contains  one  or  more 
nonnegative  numbers  and  (ii)  the  sum  of  the  negative  numbers  (if  any)  in  each  sublist 
does  not  exceed  the  largest  positive  number  in  the  sublist.  Then  a  is  the  spectrum  of 
a  nonnegative  matrix. 

The  proof  of  this  corollary  follows  from  applying  Corollary  2  to  each  of  the  sublists 

and  taking  a  direct  sum  of  nonnegative  matrices. 

Suleimanova  also  presented  the  following  conjecture  in  1965  regarding  the  RNIEP. 
Conjecture:  If  the  numbers  1,  Aa, . .  • ,  An  satisfy  1  +  '£(^=2  >  0,  for  /c  =  1, 2, 3, . . . , 

and  lAj]  <  1  for  i  =  2, 3, . . . ,  n,  then  the  numbers  are  eigenvalues  of  some  stochastic 

matrix. 
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This  conjecture  is  not  satisfied  for  n  ^  5  as  the  list 
peated  Perron  root,  satisfies  the  hypotheses  but  cannot  be  the  spectrum  of  a  non¬ 
negative  matrix  as  seen  in  Chapter  3.  Although  this  conjecture  is  not  true,  it  was 
Suleimanova’s  work  which  inspired  many  new  results  in  the  RNIEP.  One  improvement 
to  Suleimanova’s  condition  comes  from  Salzmann  [30]. 

Theorem  17  (Salzmann’s  Theorem).  [30]  iet  Ai  >  A2  ^  ^  be  real  numbers 

such  that 

1.  I  (Aj  -t-  A„_i+i)  <  i  E  Aj  when  i  =  2,  3,  . . . ,  [^] ,  and 

2.  EAi^O. 

7=1 

Then  there  exists  an  n  x  n  nonnegative  diagonalizable  matrix  having 
a  =  {Ai,  A2,  A„}  as  its  spectrum.  Furthermore  if  the  inequalities  (1)  and  (2) 

are  strict,  then  there  exists  a  positive  diagonalizable  matrix  with  spectrum  a. 

Salzmann’s  condition  clearly  extends  Suleimanova’s  condition.  Moreover,  if  we 
consider  a  =  {7, 2, -2, -3, -3},  by  Corollary  3,  there  is  a  nonnegative  reducible 
matrix  A  with  spectrum  a.  However  Salzmann’s  theorem  shows  that  A  must  be 
similar  to  a  positive  and  thus  an  irreducible  matrix. 

Kellogg  [15]  produced  a  set  of  sufficient  conditions  to  solve  the  RNIEP.  He  noticed 
that  if  there  exists  a  real  matrix  A  having  spectrum  cr  =  {Ai,  A2, .  •  • ,  An}  such  that 
Ai  >  [Ajl  for  2  <  j  ^  n,  then  A  is  similar  to  a  positive  matrix  provided  that  as 
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elements  meet  certain  conditions.  Theorem  18  contains  Kellogg’s  sufficient  conditions 
for  a  real  matrix  having  spectrum  a  to  be  similar  to  a  positive  matrix. 

Theorem  18.  [15]  Let  A  be  a  real  matrix  having  real  eigenvalues  Ai,  A2, .  •  • ,  A„  such 
that  Ai  >  A2  ^  ^  A„  and  Ai  >  1A„|.  Let  r  be  the  greatest  integer  with  K  ^  0. 

Define  Si  =  An+2-i  for  i  =  2, . . . ,  s  =  n  -  r  +  1  and 

K  =  {ie  {2,...,  min  {r,  s}}  :  A,  +  <  0}.  If  M  +  Si)  +  Sh  >  0  for  all 

heK  andifXi+Y,  >  0  ^  “  positive  matrix. 

ieK  j=r+l 

Kellogg  proved  Theorem  18  by  using  the  following  ideas.  Since  A  is  a  real  matrix 
having  real  eigenvalues  Ai,  A2, . .  • ,  An  such  that  Ai  >  A2  ^  ^  An  and  Ai  >  |Anl, 

then  A  is  similar  to  a  strictly  positive  matrix  if  and  only  if  there  is  an  (n  -  l)-simplex 
S  in  the  set  of  vectors  orthogonal  to  A^’s  eigenspace  for  Ai  such  that  AS  C  Ai5o. 
Here  5o  is  defined  to  be  the  interior  of  S  where  Sq  excludes  any  points  of  S  that 
require  a  zero  coefficient  on  any  of  the  n  —  1  vertices  of  S  in  writing  the  convex  sum 
to  obtain  the  points.  Kellogg  produced  a  real  (n  -  1)  x  (n  -  1)  matrix  A  having  the 
normalized  eigenvalues  1,  •  •  •>  ^  showed  that  A  takes  S  into  Sq. 

It  was  Fiedler  [7]  who  observed  that  changing  the  strict  inequalities  in  Kellogg  s 
conditions  allowed  Theorem  18  to  apply  to  nonnegative  matrices.  In  Theorem  19, 
Fiedler  gives  the  extenstion  of  Theorem  18  for  the  nonnegative  case.  We  call  this  as 
Kellogg’s  Theorem. 
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Theorem  19  (Kellogg’s  Theorem).  [7]  Let  Ai  ^  A2  ^  ^  A^,  Ai  ^  |A„|,  and 

let  r  be  the  greatest  integer  with  A^  ^  0.  Let  5i  =  Xn+2-i  for  i  =  2,...,s  =  n  —  r  +  1. 
Define  K  =  {i  €  {2, . . . ,  min  {r,  s}}  :  Xi  +  5i  <  0}.  If  Xi  +  +  <^i)  +  ^  0 

ieK,i<h 

S 

for  all  h  £  K  and  if  Xi  +  +  <^1)  +  ^  ^  then  a  —  {Ai,  A2,  . . . ,  A„}  is 

ieK  j=r+l 

the  spectrum  of  some  n  x  n  nonnegative  matrix. 

In  the  following,  Fiedler  [7]  also  gives  sufficient  conditions  for  a  set  of  real  numbers 
to  be  the  spectrum  of  a  nonnegative  matrix. 

Theorem  20  (Fiedler’s  Theorem).  [7]  // Ai  ^  A2  ^  ^  An  and  Ai  +  A„  + 

E  >  IE  l^i  +  An+i-il,  then  a  =  {Ai,  A2,  . . . ,  A„}  is  the  spectrum  of  some  n  x 

j=l  i=2 

n  nonnegative  matrix. 

Fiedler  proved  that  if  the  conditions  of  Salzmann’s  Theorem  are  satisfied,  then  the 
conditions  of  Fiedler’s  Theorem  are  satisfied;  moreover,  if  the  conditions  of  Fiedler’s 
Theorem  are  satisfied  then  the  conditions  of  Kellogg’s  Theorem  are  satisfied.  That 
is  Salzmann  Fiedler  Kellogg.  None  of  the  converses  are  true.  To  see  this 
consider  ui  =  {5,4,  -2,  -4}  and  (72  =  {3, 2,  -1,  -1,  -3}.  Here  notice  that  ai  satisfies 
Fiedler’s  condition  but  not  Salzmann’s  and  that  (72  satisfies  Kellogg’s  condition  but 
not  Fiedler’s.  Therefore  the  most  general  sufficient  condition  from  the  above  theorems 
comes  from  Kellogg. 

But  even  Kellogg’s  Theorem  has  its  limitations.  Notice  that  the  set 
a  =  {1,0.6,  -0.3,  -0.35,  -0.45,  -0.5}  does  not  satisfy  Kellogg’s  Theorem,  and  conse¬ 
quently  none  of  the  prior  theorems  can  identify  (7  is  a  spectrum  of  some  nonnegative 
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matrix.  As  it  turns  out,  a  is  the  spectrum  of  a  6  x  6  nonnegative  matrix.  Borobia 
proved  this.  Borobia  [2]  was  able  to  extended  Kellogg’s  conditions  in  Theorem  21. 

Theorem  21  (Borobia’s  Theorem).  [2]  Let  o  =  {Ai,  A2,  . .  ■ ,  A„}  he  a  set  of  real 
numbers  with  Ai  ^  Aa  ^  ^  A„.  Let  r  be  the  greatest  integer  with  Ar  >  0.  If  some 

partition  Ji  U  •  •  •  U  Jj  of  J  =  {Ar+i,  . . . ,  A^}  exists  such  that  Ai  ^  A2  ^  ^  A^  ^ 

E  >  T,  ^  satisfies  Kellogg’s  Theorem,  then  a  is  the  spectrum  of 

X^Ja  AC.7s  — 1  X^Ji 

some  n  X  n  nonnegative  matrix. 

To  see  an  application  of  Borobia’s  Theorem  we  return  to  the  example 
a  =  {1,0.6, -0.3, -0.35, -0.45, -0.5}.  Kellogg’s  Theorem  fails  for  this  a  because 
here  we  have  AT  =  0  as  0.6  +  (-0.5)  ^  0  but  1  +  (-0.3)  +  (-0.35)  +  (-0.45)  ^  0. 
Therefore  we  partition  a  in  the  following  way;  (1,0.6,  (-0.3)  +  (-0.35),  (-0.45)  + 
(—0.5)}.  Now  applying  Kellogg’s  Theorem  to  (1,0.6,  —0.65,  —0.95},  we  have  K  —  (2} 
ag  0.6  +  (-0.95)  <  0.  So  checking  the  remaining  conditions  we  have  1  +  (-0.95)  ^  0 
and  1  +  0.6  +  (-0.65)  +  (-0.95)  ^  0;  therefore,  there  exists  a  nonnegative  matrix 
realizing  a  by  Borobia’s  Theorem. 

5.2  Symmetric  Nonnegative  Inverse  Eigenvalue  Problem  (SNIEP) 

One  additional  restriction  to  the  RNIEP  gives  the  SNIEP.  The  SNIEP  requires 
finding  necessary  and  sufficient  conditions  on  a  =  (Ai,  Aa,  •  •  • ,  A„}  C  R  so  that 
there  exists  an  n  x  n  symmetric  nonnegative  matrix  having  a  as  its  spectrum. 

Fiedler  [7]  proved  that  the  conditions  of  Kellogg’s  Theorem  are  also  sufficient  for 
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the  existence  of  a  symmetric  nonnegative  matrix.  This  result  is  in  Theorem  22. 

Theorem  22,  [7]  Let  Ai  ^  A2  ^  ^  A„  satisfy  the  conditions  of  Kellogg’s  The¬ 

orem.  Then  there  exists  an  n  x  n  symmetric  nonnegative  matrix  with  eigenvalues 
Ai,  A2,  . . . ,  A„. 


Soules  [31]  gives  a  sufficient  condition  using  the  all  one’s  vector  as  the  Perron  vector 
for  the  existence  of  a  symmetric  nonnegative  matrix  having  real  a  as  its  spectrum. 
His  conditions  are  in  Theorem  23.  We  will  discuss  Soules  construction  further  in 
Chapter  6. 


Theorem  23.  [31]  Let  [1, 1, . . . ,  1]^  be  the  Perron  vector  corresponding  to  the  Perron 
root  Ai.  Let  Ai  ^  A2  ^  ^  A„  and  A2  +  X^fc=i  (k+iVt  ^  ®  where 


^fc=l  {k+l)k 


=  ^^] ,  then  there  exists  an  n  x  n  symmetric,  doubly  stochastic  matrix  A  such 


that  {Ai,  A2, . . . ,  An}  is  the  spectrum  of  XiA. 


Radwan  [27]  shows  that  the  conditions  of  Theorem  22  and  the  conditions  of  The¬ 
orem  23  are  not  comparable.  For  example,  ci  =  {5, 3,  —2,  —2,  —2,  —2}  satisfies 
Theorem  23  but  not  Theorem  22,  and  (T2  =  {5, 3,  —2,  —2,  —4}  satisfies  Theorem  22 
but  not  Theorem  23.  Radwan  also  points  out  that  Theorem  22  does  not  work  for 
a  =  {9,6,  —2,  —4,  —4,  —4};  however,  the  following  Theorem  24  does. 


Theorem  24  (Loewy’s  Theorem).  [27]  Let  Ai  ^  A2  >  0  ^  A3  ^  ^  An 

real  numbers.  If  for  some  partition  Ji  U  J2  0/  J  =  {—A3,  — A4, . . . ,  —An),  we  have 
Ai  ^  Y,\€Ji  ^  ^  Sj=i  ^  {Ai,  A2, .  • . ,  An}  is  the  spectrum  of 

some  n  X  n  symmetric  nonnegative  matrix. 
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As  Loewy’s  Theorem  is  applicable  to  a  cr  containing  exactly  two  positive  elements, 
the  question  arises  as  to  whether  or  not  a  general  partitioning  can  be  used  for  the 
SNIEP.  Here  is  where  a  result  similar  to  Borobia’s  Theorem  for  the  RNIEP  would 
be  useful  in  the  SNIEP.  Radwan  answers  this  question.  Radwan  [27]  proved  that  the 
conditions  of  Borobia’s  Theorem  are  also  sufficient  for  the  existence  of  a  symmetric 
nonnegative  matrix.  This  result  is  in  the  following  Theorem  25. 

Theorem  25  (Radwan’s  Theorem).  [27]  Le^  Ai  ^  A2  ^  ^  Am  >  0  nonneg¬ 
ative  numbers  and  i^i  ^  [j,2  ^  >  0  be  positive  numbers.  Let  S  be  an  integer 

such  that  l^S  ^N.  Lef  Ji  U  J2  U  •  •  •  U  6e  an  S-partition  of  J  =  ^^2,  ■  •  • ,  I^n} 

where  E  M  ^  E  ^  ^  E  A*-  Define  Tj  =  E  M  for  j  =  1,  2,  ...,  S.  If 

iieJi  tieJi 

crs  =  {Ai,  A2,  . . . ,  Am,  -Ts,  -Ts-u  ■  •  • ,  satisfies  the  hypotheses  of  Kellogg’s 

Theorem,  then  a  =  {Ai,  A2,  ...,  Am,  -Tn,  •••,  “Mi}  ^^e  spectrum  of 

some  (M  +  A  +  1)  X  (M  +  A  +  1)  symmetric  nonnegative  matrix. 

A  nice  result  due  to  Fiedler  helps  in  the  construction  of  symmetric  nonnegative 
matrices  when  spectra  of  two  symmetric  nonnegative  matrices  are  known.  In  [7],  we 
find  the  following  two  theorems  to  be  very  useful  in  this  construction.  We  use  5„ 
to  denote  the  set  of  lists  of  length  n  realizable  by  an  n  x  n  symmetric  nonnegative 

matrix. 

Theorem  26  (Fiedler’s  Symmetric  Theorem).  [7]  If  {ai,  02, ... ,  ^m}  G  Sm  with 
as  the  Perron  root,  {fiup2, . . . ,  e  with  fix  as  the  Perron  root,  and  Oi  ^  fix, 
then  for  any  cr  ^  0,  {oi  +  cr,  —  cr,  02, . . . ,  ^2,  •  ■  ■  >  Pn}  ^  Sm+n- 
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Theorem  27.  [7]  Let  A  be  a  symmetric  m  x  m  matrix  with  eigenvalues 
and  let  u,  l|ull2  =  I,  be  a  unit  eigenvector  corresponding  to  ai;  let  B  be  a  symmetric 
n  TL  n  matrix  with  eigenvalues  o.'^d  let  v,  ||'y||2  =  1,  o  unit  eigenvector 

corresponding  to  /Si .  Then  for  any  p,  the  matrix 

A 


C  = 


pvviT 


has  eigenvalues  Oi, . . . ,  am,  -  ■  ■  ,Pn,  7i)  72?  where  71, 72  are  eigenvalues  of  the  ma¬ 


trix 


P 

Pi 


5.3  The  RNIEP  and  the  SNIEP  are  Different 

The  RNIEP  and  the  SNIEP  are  not  equivalent  problems.  In  the  case  when  a  = 
{Ai,  A2, . . . ,  An}  is  real  we  consider  both  the  RNIEP  and  the  SNIEP.  We  will  restate 
the  question: 

1.  the  real  nonnegative  inverse  eigenvalue  problem  (RNIEP)  asks  which  lists  of  n 
real  numbers  occur  as  the  spectrum  of  a  nonnegative  matrix  A  G  Mn,  3-iid 

2.  the  symmetric  nonnegative  inverse  eigenvalue  problem  (SNIEP)  asks  which  lists 
of  n  real  numbers  occur  as  the  spectrum  of  a  symmetric  nonnegative  matrix 

A  e  Mn- 

Obviously,  if  the  SNIEP  has  a  solution  for  a  given  a  then  the  RNIEP  also  has  a 
solution.  In  particular,  in  the  case  Ai  ^  A2  ^  ^  ^  0  then  both  problems 
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trivially  have  a  solution,  namely  the  diagonal  matrix  D  =  diag  (Ai,  A2,  •  •  • ,  A„).  When 
at  least  one  of  the  A’s  is  negative,  are  these  problems  equivalent?  In  other  words,  if 
there  is  a  solution  for  the  RNIEP  does  it  imply  there  is  a  solution  for  the  SNIEP? 
This  question  was  raised  by  Hershkowitz  [10]  in  1978. 

We  have  seen  several  sufficient  conditions  for  the  RNIEP  in  the  previous  sections. 
We  also  have  seen  that  all  those  conditions  actually  solve  the  SNIEP.  But  in  1996, 
Johnson,  Laffey,  and  Loewy  showed  that  the  two  problems  are  different. 

The  following  is  a  brief  description  of  the  result  of  Johnson,  Laffey,  and  Loewy 
[13].  According  to  Boyle  and  Handelman  [3]  if  Ai, . . . ,  A„  meet  certain  conditions, 
it  is  possible  to  append  sufficiently  many  O’s;  A^^-i  =  )^n+2  —  •  •  •  —  Am  —  0  so 
that  Ai, . . . ,  An,  A„+1,  An+2,  •  •  •  ,  Am  are  the  eigenvalues  of  a  positive  matrix  A  G  M„. 
However,  m  and  rank  A  may  have  to  be  large  relative  to  n.  It  is  possible  to  give 
examples  for  which  m  must  be  very  large.  If  A  were  a  symmetric  nonnegative  matrix, 
then  rank  A  would  be  very  low  relative  to  m.  In  [13],  it  is  shown  that  if  A  G  Mm  is 
a  sufficiently  low  rank  nonnegative  matrix,  then  there  is  another  nonnegative  matrix 
a'  G  M^'  with  m  <m  such  that  the  nonzero  part  of  the  spectrum  of  A  agrees  with 
that  of  A.  This  rules  out  the  possibility  of  the  SNIEP  having  a  solution  for  certain 
lists  for  which  the  Boyle  and  Handelman  result  guarantees  a  solution  to  the  RNIEP. 

As  a  specific  example,  consider  n  =  6  and  Ai  =  +  e,  A2  =  A3  =  A4  =  1,  and 

As  =  Ae  =  -3  for  small  e  >  0.  For  e  >  0,  {Ai,  A2,  A3,  A4,  A5,  As}  satisfies  the  Boyle 
and  Handelman  conditions.  For  each  m  >  6,  there  is  an  €m  >  0  such  that  at  least 
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m  -  6  O’s  must  be  added  to  obtain  the  nonnegative  realization  A  e  Mm  guaranteed 
by  Boyle  and  Handelman.  Consider  m  -  22.  For  some  £22  >  0,  Johnson,  Laffey 
and  Loewy  show  that  there  is  no  symmetric  nonnegative  matrix  of  any  order  whose 

nonzero  spectrum  is  Ai,  A2,  A3,  A4,  A5,  Ae- 

In  Chapter  6,  we  will  discuss  the  results  that  show  for  lists  of  length  n  ^  4  the 
RNIEP  and  the  SNIEP  are  equivalent.  But  knowing  that  the  RNIEP  and  the  SNIEP 
are  different  in  general,  what  is  the  smallest  list  for  which  they  are  different?  We  will 
show  in  Chapter  6  that  the  RNIEP  and  the  SNIEP  are  different  for  lists  of  length 

five. 
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CHAPTER  6 


Lists  of  Length  Five 

In  this  chapter  we  derive  new  results  concerning  the  RNIEP  and  the  SNIEP.  Recall 
that  the  RNIEP  is  to  determine  the  necessary  and  sufficient  conditions  for  a  list  of  n 
real  numbers,  a  =  {Ai,  Aa, .  •  • ,  A„},  to  be  the  spectrum  of  some  nonnegative  matrix. 
The  SNIEP  is  similar  with  the  only  difference  being  that  a  must  now  be  the  spectrum 
of  some  symmetric  nonnegative  matrix.  In  this  chapter  we  will  consider  the  values 
of  a  to  always  be  ordered  such  that  Ai  ^  Aa  ^  ^  An  where  Ai  ^  |An|-  In  other 

words,  Ai  will  be  the  Perron  root  should  a  nonnegative  matrix  exist.  Of  course  we 
already  have  seen  some  necessary  conditions  for  a  real  a  such  as  the  moments  must 
be  nonnegative,  the  Perron  root  must  be  a  part  of  a,  and  the  JLL  condition  must 
hold.  These  conditions  are  sufficient  in  lists  of  four  or  fewer  real  numbers  and  we  will 
see  this  in  the  section  (6.2).  These  conditions  are  not  sufficient  in  lists  of  five  or  more 
real  numbers,  and  we  will  see  this  in  section  (6.2).  We  begin  by  stating  some  recent 
results  concerning  the  RNIEP  and  the  SNIEP  in  the  case  where  n  ^  4  in  section 
(6.1).  It  is  for  these  smaller  dimensions  that  the  problems  are  equivalent. 
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6.1  The  Soules  Approach  to  the  SNIEP 


Two  recent  papers,  McDonald  and  Neumann  [21]  and  Knudsen  and  McDonald 
[16],  discuss  the  Soules  approach  to  the  SNIEP  for  n  ^  5.  A  list  of  n  real  numbers, 

{1,  A2, . . . ,  A„}  satisfying  the  conditions  1  ^  |Ai|  for  i  =  2, 3, . . . ,  n  and  1  +  A2  H - f- 

An  ^  0  define  a  polytope  in  BP  called  the  trace  nonnegative  polytope  denoted  by 
Tn-  Any  list  of  n  real  numbers  realizable  as  the  spectrum  of  a  symmetric  nonnegative 
matrix  must  be  in  this  polytope  r„,  however  for  n  ^  5,  r„  contains  points  which  are  not 
spectrum  of  any  symmetric  nonnegative  matrices.  For  example,  {1, 1,  — 
and  {1,|,|,— 1  —  1}  are  not  the  spectrum  of  any  nonnegative  matrix  as  the  first  one 
fails  to  be  reducible  as  there  is  a  repeating  Perron  root  and  the  second  one  fails  the 
moment  condition  as  S3  ^  0.  Therefore  we  denote  by  Sn  the  set  of  points  in  that 
are  realized  as  spectrums  of  symmetric,  nonnegative  matrices. 

Soules  [31]  developed  a  method  of  constructing  symmetric  nonnegative  matrices 
whose  eigenvalues  satisfy  Ai  ^  A2  ^  ^  A„.  Starting  with  a  vector  a;  >  0,  Soules 

shows  how  to  .construct  an  orthogonal  matrix  R,  called  the  Soules  matrix,  such  that 
for  A  =  diag  {\i,  X2, . . . ,  Xn)  the  matrix  A  =  RAR'^  has  nonnegative  off-diagonal 
entries.  He  shows  that  all  off-diagonal  entries  are  positive  if  and  only  if  Ai  >  A2. 
Soules  then  finds  sufficient  conditions  on  Aj’s  which  guarantee  the  diagonal  entries  of 
A  =  RAR'^  are  also  nonnegative.  In  this  construction  from  Chapter  5  Theorem  23, 
the  positive  vector  x  becomes  the  Perron  vector  of  A  corresponding  to  Ai. 

Let  the  Soules  set  be  the  set  of  all  lists  {Ai,  A2,  •  •  • ,  A^}  for  which  there  is  a 
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Soules  matrix  R  such  that  RAR^  ^  0  where  A  —  diag  (Ai,  A2,  •  •  • ,  A^).  For  n  3 
and  n  =  4,  McDonald  and  Newmann  [21]  show  that  <Sn  =  Sn  =  In  Knudsen 

and  McDonald  it  is  shown  that  {1,  ™ 

©5  showing  that  65  is  properly  contained  in  S5.  We  saw  earlier  that  S5  is  properly 

contained  in  T5. 

The  following  interesting  lemma  by  McDonald  and  Neumann  [21]  gives  a  condi¬ 
tion  relative  to  the  spectrum  of  a  5  x  5  irreducible  symmetric  nonnegative  matrix. 
Although  it  is  a  simple  condition,  it  is  a  useful  result  for  our  work  in  section  (6.3). 

Lemma  1.  [21]  Let  A  6e  o  5  x  5  irreducible  nonnegative  symmetric  matrix  with 
eigenvalues  {1,  A2,A3,  A4,  A5}  where  1  ^  A2  ^  A3  ^  A4  ^  A5  ^  —1.  Then  —  (A5)  -t- 
trace  (A)  ^  A2.  Moreover,  if  some  eigenvector  of  A  corresponding  to  A2  has  only  one 
positive  or  one  negative  entry,  then  trace  (A)  ^  A2. 

Proof.  If  As  ^  0,  then  — A5  +  trace  (A)  =  1  -t-  A2  +  A3  +  A4  ^  A2  since  A  has  all 
nonnegative  eigenvalues. 

If  As  <  0,  then  let  B  =  A  -  As/  where  A  =  [a^j]  for  i,j  =  1,2, 3, 4, 5.  By  this 
definition,  since  B  has  all  nonnegative  eigenvalues  and  is  Hermitian,  B  is  a  positive 
semidefinite  matrix.  Since  As  is  negative,  B  is  also  a  symmetric  nonnegative  matrix. 
Let  X  be  the  Perron  vector  for  the  Perron  root,  1  —  As,  associated  with  matrix  B. 
Then  by  the  Perron-Frobenius  Theorem  x  >  0  and  can  be  chosen  so  that  |l2;ll2  —  1- 
Let  y  be  the  eigenvector  corresponding  to  A2  —  As  associated  with  matrix  B.  Since  B 
is  symmetric,  x'^y  =  0.  This  forces  y  to  have  at  least  one  negative  entry  and  at  least 
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one  positive  entry.  As  y  has  five  entries,  this  means  y  or  —y  has  at  most  two  positive 
entries.  Let  us  assume  y  has  at  most  two  positive  entries  located  at  the  top  of  the 
vector.  Now  let  y  be  represented  by  the  two  vectors  yi  and  y2  where  yi  contains  the 
positive  part  of  y  and  y2  contains  the  nonpositive  part  of  y.  Partition  B  appropriately 


for  y  so  that 


Consider  the  two  cases: 


1.  Suppose  yi  contains  exactly  two  elements.  Then  we  partition  B  so  that  Bn  is 
2x2  and  B22  is  3  x  3.  Thus  we  have  Bnyi  +  -6122/2  =  (^2  —  ^5)  Hi-  Then 
as  Bi2y2  ^  0,  we  have  Buyi  ^  (A2-A5)yi.  Since  B  is  nonnegative  and 
semidefinite,  so  is  Bn  and  we  can  conclude  that  by  the  definition  oi  B  =  A—X5I 
that 

trace  (A)  —  2A5  ^  an  +  0,22  —  2A5 
=  trace  (^n) 

^  P  (^11) 

^  A2  —  A5. 

Therefore,  —  (A5)  +  trace  (A)  ^  A2. 
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2.  Suppose  2/1  contains  one  element.  Then  by  a  similar  argument  where  Bu  is  a 
1x1  matrix,  we  have 

trace  (.4)  —  As  ^  On  —  As 

=  trace  (Bn) 

^  P  {Bn) 

^  A2  —  As 

giving  trace  {A)  ^  A2. 

□ 

In  section  (6.3)  we  show  that,  except  for  two  partially  known  cases,  every  other 
possible  list  of  five  real  numbers  satisfying  certain  sufficient  conditions  are  in  85- 

6.2  Equivalence  of  RNIEP  and  SNIEP  for  n  ^  4 

In  the  case  of  lists  of  four  real  numbers,  Loewy  and  London  [20]  have  a  construc¬ 
tion  technique  solving  the  RNIEP.  This  is  discussed  in  Theorem  28.  The  following 
symmetric  constructions  are  from  Wuwen’s  paper  [34].  In  these  constructions,  Wuwen 
uses  Theorem  27  due  to  Fiedler. 

1.  Consider  a  list  of  length  one.  Clearly  Ai  ^  0  and  hence  A  =  will  do. 

2.  For  a  list  of  length  two,  a  —  {Ai,  A2}  must  satisfy  Ai  +  A2  ^  0  where  Ai  ^  1A2|- 
An  example  of  a  symmetric  nonnegative  matrix  with  a  as  its  spectrum  is: 

Ai  +  A2  Ai  —  A2 
Ai  —  A2  Ai  -f  A2 
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3.  Now  we  consider  a  list  of  length  three.  In  this  case  where  a  —  {Ai,  A2,  A3}  C  J?, 
+  A2  4-  A3  ^  0,  Ai  >  I  Ail  for  i  =  2,  3  and  Ai  ^  A2  >  A3  the  following 
possibilities  arise: 


(a)  if  Ai  ^  A2  ^  0  ^  A3,  then  an  example  of  a  symmetric  nonnegative  matrix 
realizing  a  as  its  spectrum  is: 


Ai+A3  Ai-Aa  n 

2  2 

Ai  —A3  Ai  H-A3  n 
2  2  ^ 


0  0  A2 

(b)  if  Ai  >  0  >  A2  ^  A3,  then  an  example  of  a  symmetric  nonnegative  matrix 
realizing  a  as  its  spectrum  is: 


Ai  4-A2+A3 

Ai  +A2  “A3 

■y-2AiA2 

2 

2 

2 

Ai+Aa-Aa 

Ai+A2+Aa 

^/-2AlA2 

2 

2 

2 

V— 2Ai  A2 

2 

1 

to 

to 

0 

4.  Next  we  present  two  results  for  the  case  of  n  =  4.  Here  a  —  {Ai,  A2,  A3,  A4}  C 
R,  Ai  ^  A2  ^  A3  ^  A4,  and  Ai  ^  |Ail  for  i  =  2,  3,  4.  The  RNIEP  solution  given 
in  the  next  theorem  is  due  to  Loewy  and  London  [20]. 

Theorem  28.  [20]  Let  a  =  {Ai,  A2,  A3,  A4}  he  a  set  of  four  real  numbers,  and 
assume  that  a  satisfies  Tnaxi^i^4  |Aij  €  a  and  si  =  Ai  +  A2  +  A3  +  A4  ^  0.  Then 
a  is  the  set  of  eigenvalues  of  a  nonnegative  matrix  A. 


Proof  Assume  that  Ai  ^  A2  ^  A3  ^  A4  and  we  have  the  Perron  root  require¬ 
ment,  |Ai|  ^  Ai,  for  i  =  2,3,4. 
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If  Aj,  Aj,  A4  >  0,  -4  is  the  matrix  A  =  diag  (Ai.  Az,  A3,  A4) 

If  A2,  A3,  A4  ^  0.  the  assertion  follows  from  Theorem  16. 

If  A4  ^  0,:  A2,  A3  >  0,  cT  is  realized  by  the  nonnegative  matrix 

!  Aj  4*  A4  Ax  0  0 

-A4  0  0  0 

i  ^4  = 

0  0  A2  0 

0  0  0  A3 

If  A3,  A4  ^  0,  A2  >  0  we  distinguish  between  two  cases: 

(a)  Az  -H  A3  ^  0,  and 

(b)  A2  "H  A3  <  0. 

In  case  (4),  o'  is  realized  by  the  nonnegative  matrix 

j  Aj  +  A3  A3  0  0 

i  '-As  0  0  0 

;  .4  = 

0  0  A;  -f-  A4  Ai 

!  0  0  -A4  0 
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In  case  (b),  a  is  realized  by  the  matrix  A  =  Udiag  (Ai,  A2,  A3,  A4)  U'^  where 
U  is  the  orthogonal  matrix 

1111 
-1-111 
1  -1-11 
-11-11 

The  entries  of  A  are, 

Oil  =  022  —  O33  =  044  =  Ai  +  A2  +  A3  +  A4  ^  0, 

012  —  <^21  —  O34  =  O43  =  — Ai  —  A2  +  A3  +  A4  ^  0, 

Oi3  =  031  =  O24  =  042  =  Ai  —  A2  —  A3  +  A4  ^  0, 

0i4  =  041  =  023  =  032  =  — Ai  +  A2  —  A3  +  A4  ^  0, 

and  hence  A  is  a  nonnegative  matrix.  Cl 

The  solution  to  the  SNIEP  for  a  list  of  four  real  numbers  was  presented  by 
Wuwen  [34].  The  construction  techniques  require  the  same  sufficient  conditions 
as  those  in  Theorem  28. 

(a)  If  Ai  ^  A2  ^  A3  ^  A4  ^  0  then  A  =  diag  (Ai,  A2,  A3,  A4)  realizes  a  as  its 
spectrum. 

(b)  It  Ai  ^  A2  ^  A3  >  0  >  A4,  then  an  example  of  a  symmetric  nonnegative 
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matrix  realizing  a  as  its  spectrum  is: 


"1“  A4  Ai  A4  0  0 

Ai  —  A4  Ai  +  A4  0  0 

0  0  2A2  0 


0 


0 


0  2A3 


(c)  If  Ai  ^  A2  ^  0  >  A3  ^  A4  and  A2  ^  IA3I,  then  an  example  of  a  symmetric 
nonnegative  matrix  realizing  a  as  its  spectrum  is: 


Xi  +  A4  \i  A4  0  0 

Ai  ““  A4  Ai  +  A4  0  0 

0  0  A2  +  A3  A2  “■  A3 

0  0  A2  A3  A2  +  A3 


(d)  If  Ai  ^  A2  ^  0  >  A3  ^  A4  and  A2  <  IA3I,  then  an  example  of  a  symmetric 
nonnegative  matrix  realizing  a  as  its  spectrum  is: 


Ai  +  A2  +  A3  +  A4  Ai  +  A2  +  A3  —  A4 
Ai  +  A2  +  A3  —  A4  Ai  +  A2  +  A3  +  A4 

5  5 


5 

5 

0 

-2A3 


6 

5 

-2A3 

0 


where  S  =  y— (A2  +  A3)  (Ai  +  A3). 
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(e)  If  Ai  ^  0  >  Aa  ^  As  ^  A4,  using 


B  = 


Ai  +  A2  +  A3  +  A4  Ai  +  A2  +  A3  +  A4  yj  2  (Ai  +  A2)  A3 


Ai  +  Aa  +  As  +  A4  Aj  +  Aa  +  As  +  A4  y/  2  (Ai  +  Aa)  A3 
1^—2  (Ai  +  Aa)  As  \/— 2  (Ai  +  Aa)  A3  0 


1 


u  = 


y/2  (Ai  +  Aa  —  As) 


\/Ai  +  Aa 


\/Ai  +  Aa 


V— 2A3 

and  5  =  V—AjAa,  then  an  example  of  a  symmetric  nonnegative  matrix 


realizing  a  as  its  spectrum  is: 

y4  = 


B 

6u^ 


6u 

0 


This  completes  the  case  where  the  RNIEP  and  the  SNIEP  are  equivalent.  Al¬ 
though  Johnson,  Laffey  and  Loewy  [13]  proved  the  two  problems  are  not  equivalent 
for  all  dimensions,  it  was  an  existence  proof  as  described  in  Chapter  5.  In  section 
(6.3)  we  analyze  the  case  where  a  —  {Ai,  Aa,  A3,  A4,  A5}  in  terms  of  the  RNIEP  and  the 
SNIEP.  We  show  that  the  list  of  five  real  numbers  is  the  first  time  the  two  problems 
are  not  equivalent. 


6.3  Lists  of  Five  Real  Numbers 

The  following  addresses  the  SNIEP  in  the  case  where  the  lists  contain  five  real 
numbers.  That  is,  given  a  list  of  five  real  numbers  a  =  {Ai,  Aa,  A3,  A4,  A5}  such  that 
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1.  Ai  ^  A2  ^  As  ^  A4  ^  A5, 

2.  Ai  =  maxi^i^5  lAjj, 

3.  Sk'^  0  for  A:  =  1,  2,  3,  . . .  and 

4.  JLL  condition 

we  determine  sufficient  conditions  for  a  to  be  the  spectrum  of  a  5  x  5  symmetric 
nonnegative  matrix.  In  the  analysis,  we  first  attempt  to  use  the  conditions  from 
Kellogg’s  Theorem  to  determine  if  there  exists  a  nonnegative  matrix  and  thus  a 
symmetric  nonnegative  matrix  due  to  a  result  of  Fiedler.  We  do  so  by  examining 
the  pairing  of  positive  numbers  in  the  list  with  nonpositive  numbers  in  the  list.  For 
example  if  Ai  ^  A2  >  A3  >  0  ^  A4  ^  A5  then  the  pairs  in  question  are  A2  with 
As  and  A3  with  A4.  If  A2  +  A5  <  0  and  A3  +  A4  <  0  then  the  set  K  =  {2,3}  in 
Kellogg’s  conditions.  The  remaining  conditions  are  then  satisfied  if  Ai  +  As  ^  0, 
Ai  +  (A2  +  As)  +  A4  ^  0,  and  Ai  +  (A2  +  A5)  +  (A3  +  A4)  ^  0.  If  there  are  cases  where 
Kellogg’s  conditions  do  not  apply,  for  instance  in  our  example  when  Ai  +  (A2  +  As)  + 
A4  <  0,  then  we  attempt  to  partition  a  and  use  the  conditions  from  Borobia  s  Theorem 
or  Radwan’s  Theorem  to  see  if  there  exists  a  symmetric  nonnegative  matrbc.  The  final 
goal  of  each  case  is  a  construction  technique  showing  a  symmetric  nonnegative  matrix 
realizing  a  as  its  spectrum.  Many  times  we  will  use  Fiedler’s  result  from  Theorem  27. 

The  following  figures  show  the  possible  positive  and  nonpositive  distribution  of 
five  real  numbers.  The  break  down  from  each  are  the  possible  outcomes  based  on 
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Kellogg’s  conditions.  Each  branch  further  divides  out  when  we  apply  Borobia’s  par¬ 
titioning  scheme. 


A,  >Aj>A3>A4>A5 


A,  >Aj  SAj  SA4  >0>A5 
3 


0  >  A5  >  A3  +  A4  0  >  A3  +  A4  +  A5 


A|  +  A3  +  A4  <  0  A3  +  A3  +  A4  >  0 

unknown  3 


Figure  2:  List  of  Five  Real  Numbers 
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A,  >Aj>A3>A,  SA; 


Aj  +  A2  +  +  A5  <  0  Aj  +  A2  +  A4  +  Aj  ^  0 

unknown  3 

Figure  4:  List  of  Five  Real  Numbers 

Our  answers  to  the  SNIEP  for  lists  of  five  real  numbers  lie  in  the  results  of  the 
following  cases. 


6.3.1  Ai  >  0  ^  A2  ^  A3  ^  A4  ^  A5 

By  Suleimanova’s  Theorem  and  Theorem  27,  there  exists  a  5  x  5  nonnegative, 
symmetric  matrix  having  spectrum  {Ai,  A25  A3,  A4,  A5}.  A  construction  is  as  follows: 

|(Ai+A2+A34'A4+A5)  |(Ai4-A2+A3+A4“A5)  |\/-'2(Ai4“A2+A3)A4  a 

Aj  —  |(AiH-A2+A34*A4— As)  ^(Ai+A2+A3H-A4-I-A5)  I •^“2(Ai+A2+A3)A4  a 

I  2{Ai  +A2-f  A3)A4  I  y^—2(\i  +A2+A3)A4  0  b 

a  a  6  0, 
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where  a  =  b  =  Since  y4i  has  eigenvalues 

i/2(Ai+A2+A3-A4)  V2(Ai+A2+A3-A4) 

Ai  +  A2,  A3,  A4,  and  A5  where  Ai  +  A2  is  the  Perron  root,  we  know  that  the  normalized 
Perron  vector  u  of  Ai  +A2  is  nonnegative.  So  using  Theorem  27,  with  matrix  B  =  [0], 


which  has  Perron  root  0  and  Perron  vector  v  =  [1],  we  find  p  —  \/— A1A2  ^  0  which 
makes  the  following  matrix  symmetric  and  nonnegative.  Thus, 

[  Ax  puu^l 


0  J 

is  nonnegative  and  has  the  spectrum  {Ai,  A2,  A3,  A4,  A5}. 


6.3.2  Xi  ^  A2  !!>  0  ^  A3  ^  A4  ^  A5 


1.  If  A2  +  A5  <  0  then  Kellogg’s  conditions  are  satisfied.  That  is  Ai  +  A5  ^  0  and 
A^  (A2  “h  A5)  4“  A3  4“  A4  ^  0  because  the  first  is  the  Perron  root  condition  and 
the  second  is  the  trace  nonnegative  condition.  This  means  {Ai,  A2,  A3,  A4,  A5} 
is  the  spectrum  of  some  nonnegative  symmetric  matrix.  An  example  of  such  a 


matrix  using  Theorem  27  is 


Ai— S+A3+A4  Ai — e+As— A4  ^ — 2(Ai  £^^3  ^  ^ 

Ai— e+A3“A4  Ai“*e+A3+A4  2{Ai  — e)A3  6  S 

yj— 2 ( A 1  —  e ) A3  2(Ai— e)A3  0  p  p 

5  S  p  A2+£+A5  A2~I"£^'~A5 

6  6  p  A2+£— As  A2+£+A5  _ 


where  S  = 


\/e(Xi  —e)\/Xi  — £— A2 


^  >/£V^Ai— £— A2 


and  e  ^  0  satisfying  Ai  —  e  ^  A2 


and  A2  ~f"  £  T  A5  ^  0. 


2.  For  all  remaining  subcases,  A2  +  A5  ^  0  must  be  true.  So  one  of  the  following 
partitions  of  the  set  a  must  apply. 
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(a)  If  Ai  >  A2  >  0  ^  A5  ^  A3  +  A4  and  A2  +  A3  +  A4  ^  0.  Then  the  following 
naatrix  is  symmetric,  nonnegative,  and  has  spectrum  {Ai,  A2,  A3,  A4,  A5}. 

A2  H”  A3  +  A4  A2  +  A3  —  A4  \/— 2A2 A3  0  0 

A2  +  A3  —  A4  A2  +  A3  +  A4  \J —2X2^  0  0 

\f--2X2^  \/““2A2A3  0  0  0 

0  0  0  Ai  +  A5  Ai  •—  A5 

0  0  0  Ai  —  A5  Ai  +  A5 

(b)  If  Ai  ^  A2  >  0  ^  A5  ^  A3  +  A4  where  A2  +  A3  +  A4  <  0  but  Ai  +  A3  +  A4  ^  0, 
then  the  following  matrix  is  symmetric,  nonnegative  and  has  the  desired 
spectrum  {Ai,  A2,  A3,  A4,  A5}. 

A2+A5  A2—A5  000 

A2-'A5  A2+A5  0  0  0 

0  0  A1+A3  +  A4  A1+A3  — A4  2(Ai+A3)A4 

0  0  Ai+As— A4  A14-A3+A4  y/— 2(  Ai  +A3 )  A4 

0  0  y/— 2(Ai-f A3)A4  -y/— 2(Ai+A3)A^  0 

(c)  If  Ai  ^  A2  >  0  ^  As  ^  As  +  A4  and  no  previous  conditions  such  as  (2a) 
or  (2b),  then  we  can  form  certain  conclusions  about  the  problem.  First,  if 
a  solution  exists,  it  must  be  an  irreducible  matrix.  Second,  A3  /  0.  Third 
the  two  restrictions  that  0  >  A5  >  A3  4-  A4  and  A2  +  As  >  0  must  apply. 
For  the  first  conclusion,  we  rule  out  a  repeating  Perron  root.  Suppose 
Ai  =  A2,  then  the  partitioning  {Ai,A5}  and  {A2,A3,A4},  must  be  spectra 
of  2  X  2  and  3x3  nonnegative  matrices  respectively.  But  then  we  have 
Ai  +  A3  +  A4  =  A2  +  A3  +  A4  <  0  which  contradicts  the  necessary  trace 
condition.  In  the  case  of  Ai  7^  A2,  if  the  matrix  is  reducible  some  partition 
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must  occur  among  the  list  of  five  numbers.  However  the  weakest  partition 
of  {Ai,  A3,  A4}  cannot  occur  since  we  are  not  assuming  (2b).  But  any  other 
partition  is  worse  since  A3  +  A4  ^  Aj  +  Xj  for  any  i,j  =  3, 4, 5  with  i  ^  j. 
This  means  Ai  will  never  dominate  the  sum  of  any  two  of  the  nonpositive 
entries;  therefore,  no  other  partition  will  work. 

The  second  conclusion  that  A3  0  is  because  Ai  +  A3  +  A4  <  0.  As  we 
must  satisfy  the  Perron  root  condition  that  Ai  +  A4  ^  0,  adding  a  zero  will 
not  make  Ai  +  A3  +  A4  <  0. 

In  the  third  conclusion,  the  first  restriction  is  derived  from  the  inequalities 
0  >  As  ^  A4  ^  As  and  Ai  +  A3  +  A4  <  0.  As  we  no  longer  have  Perron  root 
domination  with  this  partition  and  A1+A5  ^  0,  we  realize  that  A5  >  A3+A4. 
The  second  restriction  is  apparent  from  the  nonnegative  trace  condition, 
Ai  4“  A2  4~  A3  T  A4  +  As  ^  0.  For  if  A2  T  A5  =  0  then  Ai  A3  T  A4  ^  0 
contradicting  Ai  +  A3  +  A4  <  0.  Hence  A2  +  As  >  0. 

An  additional  observation  is  that  this  case  requires  2A2  >  Ai.  We  find  this 
to  be  true  since  A2  +  A3  >0  and  A2  +  A4  >  0  giving  2A2  >  —  (A3  +  A4)  and 
since  Ai  +  A3  +  A4  <  0  giving  -  (A3  +  A4)  >  Ai.  Putting  it  all  together  we 
have  2A2  >  —  (A3  +  A4)  >  A^. 

Therefore,  the  complete  set  of  restrictions  for  this  case  are: 

i.  Aj  >  A2  >  0  >  A3  ^  A4  ^  As 

ii.  A2  +  As  >  0 
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iii.  A3  +  A4  ■<  A5,  3,nd 


iv.  Ai  +  A3  +  A4  <  0. 


As  of  this  submission,  we  have  examples  of  lists  satisfying  the  prescribed 
restrictions  and  realized  by  nonnegative  matrices  and  also  examples  of  lists 
that  meet  the  restrictions  that  cannot  be  spectrums  of  any  nonnegative 
matrix.  We  prove  that  there  does  not  exist  any  symmetric  nonnegative 
matrix  having  spectrum  with  the  prescribed  restrictions  and  also  prove 
that  if  a  nonnegative  matrix  does  exist,  it  must  be  irreducible.  The  proof 
and  details  are  in  the  conclusion  section  (6.4)  of  this  chapter. 

(d)  If  A4  ^  A2  >  0  ^  A3  +  A4  +  A5  and  A2  +  A3  +  A4  -b  A5  ^  0,  then  the  following 
matrix  is  symmetric,  nonnegative,  and  has  spectrum  {Ai,  A2,  A3,  A4,  A5}. 

A24-A3+A4+A5  A24'A3+A4— As  2(A2+A3)A4  ^  A4  ^ 

2^  A2+A34-A4— As  A2+A3+A4+A5  2(A2+A3)A4  0 

^=2  V-2(^2+A3)A4  V-2(A2+A3)A4  0  0 

v/-A2A3\/2(A2+A3)  -v/-A2  A3  v^2(A2+A3)  2y^A2A3A4  ^  ^ 

^A2  -bA3  — A4  y^A2+A3  — A4  y^A2+A3  — A4 

0  0  0  0  2Ai  ^ 


(e)  If  Ai  ^  A2  >  0  ^  A3  +  A4  +  As,  A2  +  A3  +  A4  +  As  ■<  0,  and  Ai  +  A3  +  A4  + 
As  ^  0,  then  the  following  symmetric  nonnegative  matrix  has  spectrum 
{Ai,  A2,  A3,  A4,  As}. 


Ai+A3+A4H-A5  A1+A3+A4— As  2(Ai+A3)A4  ^  ^  0 

Ai4-A3+A4— As  A1+A3+A4+A5  'y/~2(Ai4-A3)A4  ^  ^ 

V-2(Ai+A3)A4  y'_2(Ai+A3)A4  0  ° 

\/~AiA3\/2(Xi+A3)  y-Ai  A3v^2(Ai+A3)  2yAiA3A4  ^  ^ 

yAi+A3  — A4  yAi+A3  — A4  yAi+A3  — A4 

0  0  0  0  2A2 
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6.3.3 


Ai  ^  A2  ^  As  >  0  ^  A4  ^  As 


(a)  If  A2  “f"  As  <C  0,  A3  -f-  A4  <C  0,  and.  Ai  -I-  A2  4"  A4  As  ^  0,  th.6n  Ksllogg's 
conditions  are  satisfied  and  by  Theorem  27,  there  exists  a  symmetric  non¬ 
negative  matrix  having  spectrum  {Ai,  A2,  A3,  A4,  As}.  An  example  of  such 
is  the  following  matrix. 

Ai-f'A3+A4+A5  Ai+A24'A4 — As  S  SO 

1  A1+A2+A4 — As  A1+A2+A4+A5  6  SO 

A  =  -  s  s  0  -2A4  0 

2  <5  s  -2X4  0  0 

0  0  0  0  2A3  _ 

where  6  =  (A2  +  A4)  (Ai  +  A4). 

(b)  If  A2  +  As  <  0,  A3  -h  A4  ^  0,  then  the  following  matrix  has  spectrum 
{Ai,  A2,  A3,  A4,  As}. 

Aj  As  A4  —  As  0  0  0 

Ai  —  As  Ai  +  As  0  0  0 

0  0  A3  -|-  A4  A3  —  A4  0 

0  0  A3  —  A4  A3  -f  A4  0 

0  0  0  0  2A2 

(c)  If  A2  +  As  ^  0,  A3  +  A4  <  0,  then  the  following  matrix  has  spectrum 
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{Ai,  A2,  A3,  A4,  A5}. 

A2  4*  A5  A2  —  A5  0  0  0 

A2  A5  A2  +  A5  0  0  0 

0  0  Ai  +  A4  Ai  —  A4  0 

0  0  Ai  —  A4  Ai  +  A4  0 

0  0  0  0  2A3 

(d)  If  A2  +  A5  ^  0,  A3  +  A4  ^  0,  then  the  following  matrix  has  spectrum 
{Ai,  A2,  A3,  A4,  A5}. 

A2  +  A5  A2  “  As  0  0  0 

A2  —  As  A2  +  As  0  0  0 

0  0  A3  +  A4  A3  A4  0 

0  0  A3  —  A4  A3  +  A4  0 

0  0  0  0  2Ai 

(e)  If  A2  +  As  <  0,  A3  +  A4  <  0,  and  Ai  +  A2  +  A4  +  As  <  0  where  Aj  ^ 
A2  ^  A3  >  0  ^  A4  ^  As,  then  Kellogg’s  conditions  are  not  satisfied.  For 
Borobia’s  conditions  partitioning  can  only  happen  in  one  way;  namely, 
Ai  ^  A2  ^  A3  >  0  >  A4  +  As- 

Moreover,  we  can  deduce  additional  restrictions.  For  if  Ai  =  A2  it  is  not 
possible  to  partition  the  list  to  get  a  reducible  matrix  realizing  the  lists  as 
its  spectrum.  If  A4  =  0  then  Ai  +  A2  +  A4 + As  <  0  becomes  Ai  +  A2  +  As  <  0. 
Since  Ai  +  As  ^  0  and  A2  >  0  this  cannot  happen.  Therefore  we  have  the 
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following  inequalities  Ai  >  A2  ^  A3  >  0  >  A4  ^  Ag. 

Borobia’s  conditions  are  not  satisfied  in  this  partition  where  we  group  A4 
and  As  together  because  the  Perron  root  does  not  dominate  the  sum  of  A4 
and  Ag.  This  is  apparent  from  the  restriction  Ai  +  A2  +  A4  +  Ag  <  0  which 
implies  Ai  +  A4  +  Ag  <  0.  Therefore  we  still  end  up  with  the  restriction 
of  Ai  +  A4  +  Ag  <0  and  do  not  have  any  nonnegative  solution  from  this 
partition. 

As  of  this  submission,  we  have  examples  of  lists  satisfying  the  restrictions 
and  that  are  realized  by  nonnegative  or  symmetric  nonnegative  matrices. 
Soule’s  construction  of  a  symmetric  nonnegative  matrix  with  the  assump¬ 
tion  that  the  all  one’s  vector  is  the  Perron  vector  does  not  work  nor  does 
Theorem  27  as  this  set  cannot  have  a  reducible  nonnegative  matrix  for  a 
solution.  This  is  apparent  when  we  try  to  partition  the  positive  eigenval¬ 
ues  with  the  negative  ones  and  still  meet  the  necessary  conditions  for  a 
nonnegative  solution.  A  detailed  summary  is  provided  in  the  conclusion  of 
this  chapter. 
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6.3.4  ^  -^2  ^  A3  ^  A4  >  0  ^  As 
As  Ai  +  A5  ^  0,  the  symmetric  nonnegative  matrix  realizing  a  is  given  by 

Aj  4"  A5  Ai  —  A5  0  0  0 

Aj  —  A5  Ai  +  A5  0  0  0 

0  0  2A2  0  0 

0  0  0  2A3  0 

0  0  0  0  2A4 

6.3.5  Aj  ^  A2  ^  A3  ^  A4  ^  A5  >  0 
Then  A  =  diag  (Ai,  A2,  A3,  A4,  A5)  is  a  symmetric  nonnegative  matrix  realizing  a. 

6.4  Conclusion 

There  are  two  cases  from  the  above  where  we  have  examples  but  no  definite 
conclusions.  However,  in  one  of  the  cases  we  prove  there  does  not  exist  a  symmetric 
nonnegative  matrix  having  the  list  of  five  real  numbers  as  a  spectrum  and  in  the  other 
case  we  have  examples  of  symmetric  nonnegative  matrices  with  spectra  satisfying  the 
conditions.  The  following  is  a  summary  of  the  restrictions  on  the  two  cases. 

1.  Ai  >  A2  >  0  >  A3  ^  A4  ^  A5  with  restrictions 

(a)  A2  +  A5  >  0, 

(b)  A3  +  A4  <  A5,  and 

(c)  Ai  +  A3  +  A4  <  0. 
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2.  Ai  ^  A2  ^  As  >  0  ^  A4  ^  As  with  restrictions 

(a)  Ai  >  A2  ^  A3  >  0  >  A4  ^  As, 

(b)  A2  +  As  <  0, 

(c)  A3  +  A4  <  0,  and 

(d)  Ai  +  A2  +  A4  +  As  <  0. 

In  case  (1),  we  prove  that  there  does  not  exist  a  solution  to  the  SNIEP. 

Theorem  29.  If  a  =  {Ai,  A2,  A3,  A4,  As}  he  a  list  of  real  numbers  satisfying  the  fol¬ 
lowing: 

1.  Ai  >  A2  >  0  >  A3  ^  A4  ^  As, 

2.  A2  +  As  >  0, 

3.  As  >  A3  +  A4,  and 
4-  Ai  +  A3  +  A4  <  0, 

then  there  does  not  exist  a  5  x  5  symmetric  nonnegative  matrix  having  spectrum 
{Ai,  A2,  A3,  A4,  As}.  Moreover,  the  only  nonnegative  matrix  that  may  realize  a  as  its 
spectrum  must  be  irreducible. 

Proof  Without  loss  of  generality  we  assume  the  list  of  numbers  is  normalized.  Then 
there  are  two  cases  to  consider.  If  a  nonnegative  matrix  realizes  a  as  its  spectrum,  it 
would  need  to  be  reducible  or  irreducible. 
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1.  The  reducible  case;  Although  a  discussion  was  presented  in  the  prior  section 
concerning  this  case  we  repeat  it  here.  If  a  nonnegative  reducible  matrix  realizes 
(j,  then  we  should  be  able  to  partition  the  list  a  so  that  each  sublist  would  be 
the  spectrum  of  a  smaller  nonnegative  matrix.  But  this  is  impossible.  As  we 
must  put  two  nonpositive  values  with  one  positive  value  for  a  partition  to  work, 
the  smallest  negative  sum  a  positive  value  must  dominate  is  A3  +  A4  with  the 
remaining  negative  value  A5  to  be  dominated  by  the  remaining  positive  value. 
The  first  restriction  comes  from  A3  +  A4  ^  Aj  +  Xj  for  any  combination  of 
i,j  =  3,4,5  with  i  ^  j.  Considering  the  choices,  we  can  either  match  the 
Perron  root  Ai  or  A2  with  A3  +  A4.  Since  Ai  ^  A2  and  A2  +  A5  ^  0  the  best 
partition  is  {Ai,  A3,  A4}  and  {A2,  A5}.  But,  the  first  list  fails  the  trace  condition 
as  Ai  +  A3  +  A4  <  0.  Therefore  since  we  can  do  no  better  than  dominating  the 
sum  of  A3  +  A4  with  a  positive  value,  there  cannot  exist  any  nonnegative  matrix 
having  a  =  {Ai,  A2,  A3,  A4,  A5}  as  its  spectrum. 

2.  Suppose  a  is  realized  by  a  symmetric  nonnegative  matrix  which  is  irreducible. 
By  Lemma  1,  we  have  —  {X5)+trace  (A)  =  Ai-f  A2+A3+A4  ^  A2  but  this  implies 
Ai  +  A3  +  A4  ^0  which  contradicts  the  last  condition  in  the  hypothesis  of  this 
theorem.  Therefore,  there  does  not  exist  an  irreducible  symmetric  nonnegative 
matrix  realizing  a. 

So  for  a  list  of  five  real  numbers  satisfying  the  restrictions  of  Theorem  29  to  be  the 
spectrum  of  a  nonnegative  matrix,  the  solution  set  of  nonnegative  matrices  excludes 
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symmetric  matrices  and  excludes  reducible  matrices. 


□ 


The  above  theorem  does  not  answer  whether  or  not  the  restrictions  are  sufficient  for 
there  to  exist  a  5  x  5  irreducible,  nonnegative,  nonsymmetric  matrix.  We  can  prove 
they  are  not  sufficient  by  the  following  examples.  Consider  the  lists  {5, 4,  —3,  —3,  —3} 
and  {8,7,-4, -5, -6}.  They  satisfy  the  conditions  of  Theorem  29;  however,  the 
refined  JLL  necessary  condition  that  4s4  ^  s\  discussed  in  Chapter  5  fails  as  4496  ^ 
4624  for  list  one  and  69392  ^  144400  for  list  two.  Therefore  these  lists  cannot  be 
spectrums  of  some  nonnegative  5x5  matrix. 

As  discussed  in  Chapter  3,  Laffey  and  Meehan  [18]  discovered  the  necessary  and 
sufficient  conditions  for  a  list  of  5  complex  numbers  summing  to  zero  to  be  the  eigen¬ 
values  of  a  5  X  5  nonnegative  matrix.  In  this  particular  case  we  have  examples  of  lists 
of  five  real  numbers  for  which  the  conditions  of  Theorem  29  are  met  and  for  which 
Theorem  11  can  be  applied.  Consider  the  list  a  =  {97, 71, —44, —54, —70}  which 
satisfies  the  restrictions.  For  this  list  we  get  the  following  nonnegative  matrix  having 
£7  as  its  spectrum.  Mathematica  was  used  to  create  this  matrix  using  Laffey  and 


Meehan’s  construction. 

0 

1 

0 

0 

0 

12101  1  3\/418657 

2  '  2 

0 

1 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

1 

^  (14157354960  +  20548080v/418657) 

0 

228312 

12101  3^418657 

2  2 

0 
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From  this  case  we  derive  the  following  Corollary  4. 


Corollary  4.  In  the  case  of  a  list  of  five  real  numbers,  the  two  problems,  the  RNIEP 
and  the  SNIEP,  are  different. 

Proof.  Consider  the  list  a  —  {97, 71, —44, —54, —70}.  This  list  has  trace  zero  and 
satisfies  Laffey  and  Meehan’s  conditions.  Hence  there  exists  a  matrix  A  ^  0  such 
that  the  spectrum  of  A  is  the  given  list. 

The  list  also  satisfies  conditions  of  Theorem  29.  Therefore  there  does  not  exist 
a  symmetric  nonnegative  matrix  realizing  the  given  list  as  its  spectrum.  Hence  the 
SNIEP  and  the  RNIEP  are  different  for  lists  of  length  five.  □ 

In  case  (2),  we  give  examples  that  there  do  exist  symmetric  nonnegative  matrices 
and  hence  nonnegative  matrices  having  a  =  (Ai,  A2,  A3,  A4,  A5}  as  a  spectrum.  The 
following  examples  are  symmetric. 


0 

903 

43 

151 

641 

903 

94 

634 

138 

24 

43 

634 

96 

767 

1 

151 

138 

767 

44 

691 

641 

24 

1 

691 

0 

has  the  spectrum  a  =  {1657.8406, 439.1386,342.1249,  -1043.4224,  -1161.6816}. 
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40 

436 

838 

2 

167 

436 

0 

0 

177 

730 

838 

0 

18 

840 

319 

2 

177 

840 

97 

851 

167 

730 

319 

851 

4 

has  the  spectrum  a  =  {1827.3721, 358.2611, 210.4551,  -997.7358,  -1239.3515}. 

One  last  example  of  this  case  comes  from  Reams  [29].  The  list  a  =  {6, 1, 1,  —4,  —4} 
satisfies  all  the  conditions  of  Theorem  29  in  this  case  and  does  have  a  trace  zero 
nonnegative  matrix  having  a  as  its  spectrum. 

Theorem  30.  Let  a  —  {Aj,  A2,  A3,  A4,  A5}  be  a  list  of  five  real  numbers  satisfying  the 
following  conditions: 

1.  A2  +  A5  <  0, 

2.  A3  +  A4  <  0, 

3.  Ax  T  A2  T  A4  T  A5  <C  0,  and 
4-  Ai  >  A2  ^  A3  >  0  >  A4  ^  A5. 

A  nonnegative  matrix  realizing  a  cannot  be  reducible. 

Proof.  Let  a  satisfy  the  given  conditions.  For  there  to  exist  a  nonnegative  reducible 
matrix  having  a  as  its  spectrum,  it  must  be  possible  to  partition  a  into  smaller  lists 
which  are  the  spectra  of  smaller  matrices.  As  Ai  +  A5  ^  0  and  we  have  the  given 
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condition  that  Ai  +  A2  +  A4  +  A5  <  0,  the  two  imply  A2  +  A4  <  0.  Therefore  as  A2 
does  not  dominate  either  of  the  negative  values  A4  or  A5,  it  cannot  be  the  Perron  root 
for  any  smaller  matrix.  Since  Ai  +  A4  +  A5  <  0  we  cannot  partition  all  the  negative 
values  with  Ai.  Hence  if  a  matrix  should  exist  with  a  as  its  spectrum,  it  would  have 
to  be  irreducible.  □ 

As  we  have  seen,  this  problem  has  many  different  facets  to  it.  While  the  RNIEP  is 
equivalent  to  the  SNIEP  in  many  cases,  we  now  know  in  case  (1)  the  two  problems  are 
indeed  different.  It  is  unclear  whether  or  not  the  RNIEP  and  the  SNIEP  are  different 
in  the  latter  case  (2).  As  both  of  the  cases  require  irreducible  matrices  for  possible 
solutions,  more  construction  techniques  or  existence  proofs  need  to  be  developed. 
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CHAPTER  7 


Lists  of  Length  Six 

As  we  have  seen  in  Chapter  6,  the  5x5  case  is  the  first  case  in  which  the  RNIEP 
and  SNIEP  are  different.  The  process  of  analyzing  this  difference  came  from  breaking 
down  the  possible  positive  and  nonpositive  numbers  in  the  list  of  five  real  numbers 
and  trying  to  apply  Kellogg’s  and  Borobia’s  conditions.  We  continue  this  process 
on  the  list  of  six  real  numbers,  o  =  {Ai,  A2,  A3,  A4,  A5,  Ag},  that  satisfy  the  following 
conditions.  No  one  has  looked  at  this  scenario  in  the  literature  as  of  this  submission. 

1.  Ai  ^  A2  ^  A3  ^  A4  ^  A5  ^  Ag, 

2.  Ai  =  |Ai|, 

3.  Sfc  ^  0  for  A:  =  1,  2,  3,  . . .,  and 

4.  JLL  condition. 

Our  goal  is  to  identify  sufficient  conditions  that  prove  the  existence  of  symmetric 
nonnegative  matrices  having  the  spectrum  a  =  {Ai,  A2,  A3,  A4,  A5,  Ag}. 
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7.1 


Lists  of  Six  Real  Numbers 


7.1.1  Aj  >  0  ^  A2  ^  A3  ^  A4  ^  A5  ^  Ae 

In  this  case  since  Ai  dominates  the  sum  of  the  negative  eigenvalues,  we  know 
by  Suleimanova’s  condition  that  there  exists  a  nonnegative  matrix  having  a  as  its 
spectrum.  Moreover,  using  Theorem  27  we  can  construct  a  symmetric  nonnegative 
matrix  with  cr  as  its  spectrum. 

7.1.2  Ai  ^  A2  >  0  ^  A3  ^  A4  ^  A5  ^  Ae 

We  begin  our  analysis  of  this  by  applying  Kellogg’s  conditions  which  identify  those 
elements  that  belong  to  the  set  K.  In  this  distribution  of  numbers  we  have  only  the 
pairing  of  A2  +  Ag  to  analyze.  The  possible  results  from  this  lead  to  three  main  cases 
with  subcases  hinging  on  the  remaining  Kellogg’s  conditions. 

1.  If  A2  +  Ag  <  0,  then  we  have  K  =  {2}.  The  remaining  Kellogg’s  conditions  to 
check  are  that  A4  +  Ag  ^0  and  Ai  +  (A2  +  Ag)  -t-  A3  +  A4  ■!"  Ag  ^  0.  Both  of  these 
conditions  must  be  true  as  the  first  is  the  Perron  root  condition  and  the  second 
is  the  nonnegative  trace  condition.  So  A2  +  Ag  <  0  implies  Kellogg’s  conditions 
are  satisfied  so  by  Theorem  22,  we  know  there  exists  a  symmetric  nonnegative 
matrix  having  a  as  the  spectrum. 

2.  If  A2  +  Ag  >  0,  then  we  have  K  =  %.  The  remaining  condition  to  check  is 
whether  Ai  +  A3  +  A4  +  Ag  ^  0.  First  we  will  assume  this  condition  holds. 
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So  for  Ai  +  A3  +  A4  +  A5  ^  0  we  know  by  a  corollary  of  Suleimanova’s  Theo¬ 
rem,  there  exists  a  symmetric  nonnegative  4x4  matrix,  Ai,  having  spectrum 
a  (Ai)  =  {Ai,  A3,  A4,  As}.  We  also  know  there  exists  a  symmetric  nonnegative  2 
X  2  matrix,  A2,  having  spectrum  a  {A2)  =  {A2,  As}  since  A2  +  As  ^  0.  Putting 
the  two  together,  we  find  the  symmetric  nonnegative  matrix,  B, 


0  A2 


such  that  (7  (B)  —  {Ai,  A2,  A3,  A4,  A5,  As}. 


3.  If  A2  -f  As  ^  0  and 

Ai  -t-  A3  -1-  A4  H-  As  <  0,  (3) 

then  Kellogg’s  conditions  are  not  met.  We  now  try  to  partition  and  apply 
Borobia’s  conditions.  Note,  to  apply  Borobia’s  conditions  to  any  partition,  we 
must  assume  the  Perron  root  condition  holds  for  that  partition.  This  means  Ai 
is  greater  than  the  absolute  value  of  the  most  negative  sum  in  the  partition. 
When  the  Perron  root  condition  fails  for  the  partition,  we  will  list  it  as  the  last 
subcase. 


(a)  If  Ai  ^  A2  >  0  ^  As  ^  A3  -t-  A4  +  As.  Then  in  order  to  apply  Borobia’s 
conditions,  we  must  satisfy  the  Perron  root  condition,  that  is  Ai  -|-  A3  -I- 
A4  +  As  ^  0.  But  this  is  not  possible  given  (3).  Therefore  we  conclude  this 
particular  partition  gives  us  no  useful  information. 
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As  an  individual  partition  we  may  gain  nothing,  but  applying  this  idea  to 
the  general  partition  of  Ai  ^  A2  >  0  ^  Aj  ^  Aj+Afe  +  A^  where  i,j,  k,  m  are 
distinct  elements  from  {3,4, 5,6},  we  are  able  to  remove  these  partitions 
from  contention.  This  is  because  A3  +  A4  +  A5  ^  Aj  +  Aj  +  Afe  for  any 
combination  of  the  nonpositive  numbers.  Therefore  since  the  Perron  root 
condition  fails  when  considering  the  most  negative  sum  A3  +  A4  +  A5  it  will 
fail  for  any  other  sum  of  three  distinct  nonpositive  numbers  also. 

(b)  If  Ai  ^  A2  >  0  ^  A3  +  A4  +  As  ^  As.  Then  we  try  to  apply  Kellogg’s 
conditions  to  the  list  of  four  elements  a  1  =  (Ai,  A2,  A3  +  A4  +  A5,  Ag).  Since 
A2  +  As  is  nonnegative  on  the  original  a,  we  know  K  =  0.  However  we  get 
the  next  implication 

A3  +  A4  +  As  ^  Ag 

Ai  +  A3  +  A4  +  As  ^  Ai  +  Ag  ^  0 

contradicting  (3).  Thus  this  partition  cannot  work. 

(c)  If  Ai  ^  A2  >  0  >  A3  +  A4  ^  As  +  Ag.  Here  we  try  to  apply  Kellogg’s 
conditions  on  the  set  a\  =  (Ai,  A2,  A3  +  A4,  Ag  +  Ag}.  If  Ai  +  Ag  +  Ag  ^  0  we 
satisfy  the  required  Perron  root  condition.  So  assuming  this  is  true,  there 
are  two  subcases  dependent  on  the  sign  of  A2  +  Ag  +  Ag. 

i.  If  A2  +  As  +  Ag  <  0,  then  K  =  {2}.  Checking  the  remaining  condition 
of  Ai  +  As  +  Ag  ^  0  which  we’ve  assumed  to  be  true,  we  have  the 
existence  of  a  symmetric  nonnegative  matrix  having  spectrum 
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a  =  {Ai,  A2,  A3,  A4,  As,  As}  by  Radwan  [27]. 
ii.  If  A2  +  As  +  As  ^  0,  then  iiT  =  0.  So  by  Suleimanova’s  condition  and 
Fiedler,  we  know  there  exist  symmetric  nonnegative  3x3  matrices  Ai 
with  (7  (yli)  =  {Ai,A3,A4}  and  j42  with  o'[A2)  =  {A2)As,A6}.  Thus 


the  symmetric  matrix  B  defined  by 


0  A2 


has  spectrum  a. 

iii.  When  the  Perron  root  condition  is  not  satisfied,  that  is  Ai  +  As  +  As  < 
0,  we  cannot  apply  Radwan’s  result.  This  means  the  existence  of  a 
symmetric  nonnegative  matrix  or  a  nonnegative  matrix  is  unknown  in 
this  case. 


(d)  Suppose  A 1  ^  A2  >  0  ^  A3+A5  ^  A4+A6.  Letcri  —  {Ai,  A2,  As+As,  A4+A6}. 

i.  If  A2  +  A4  +  As  <  0,  then  applying  Kellogg’s  conditions  to  cti  we  have 
K  =  {2}.  Then  by  the  assumption  that  Ai  +  A4  +  As  ^  0  all  of 
Kellogg’s  conditions  are  met  and  we  have  the  existence  of  a  symmetric 
nonnegative  matrix  from  Radwan’s  Theorem. 

ii.  If  A2  +  A4  +  Ag  ^0,  then  applying  Kellogg’s  conditions  to  cri  we  have 
K  =  0.  Assuming  Ai  +  A4  +  Ag  ^  0,  we  may  construct  a  symmetric 
nonnegative  matrix  B  having  spectrum  a  from  two  smaller  3x3  sym¬ 
metric  nonnegative  matrices  Ai  and  A2  where  cr{Ai)  =  {Ai,A3,A5} 
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and  (7  (^2)  —  {^2,  A4,  Ag}.  Here 


has  spectrum  a. 


iii.  When  the  Perron  root  condition  is  not  satisfied,  that  is  Ai  +  A4  + 
Ag  <  0,  we  cannot  apply  Radwan’s  result.  This  means  the  existence 
of  a  symmetric  nonnegative  matrix  or  a  nonnegative  matrix  is  still 
unknown. 


(e)  Suppose  Ai  ^  A2  >  0  ^  As+Ag  ^  A4+A5.  Let  ai  =  {Ai,  A2,  Ag+Ag,  A4+A5}. 

i.  If  A2  +  A4  +  Ag  <  0,  then  applying  Kellogg’s  conditions  to  cji  we  have 
K  =  {2}.  Then  by  the  assumption  that  Ai  +  A4  +  Ag  ^  0  we  have  the 
existence  of  a  symmetric  nonnegative  matrix  from  Radwan’s  Theorem. 

ii.  If  A2  +  A4  +  Ag  ^0,  then  applying  Kellogg’s  conditions  to  ai  we  have 
K  =  0.  Assuming  Ai  +  A4  +  Ag  ^  0,  we  may  construct  a  symmetric 
nonnegative  matrix  B  having  spectrum  a  from  two  smaller  3x3  sym¬ 
metric  nonnegative  matrices  Ai  and  A2  where  a{Ai)  =  {Ai.AgjAg} 
and  a  {A2)  =  {A2,  A4,  Ag}.  Here 

Ai  0 

B  = 

0  A2 

has  spectrum  a. 


92 


iii.  When  the  Perron  root  condition  is  not  satisfied,  that  is  Ai  +  A4  +  A5  < 

0,  we  cannot  apply  Radwan’s  result.  This  means  the  existence  of  a 
symmetric  nonnegative  matrix  or  a  nonnegative  matrix  is  unknown. 

(f)  Suppose  Ai  ^  As  >  0  >  A4+A5  ^  As+Ag.  Let  ai  =  {Ai,  As,  A4+A5,  As+Ae}. 

i.  If  As  +  A3  +  Ae  <  0,  then  applying  Kellogg’s  conditions  to  ai  we  have 
K  =  {2}.  Then  by  the  assumption  that  Ai  +  A3  +  Ae  ^  0  we  have  the 
existence  of  a  symmetric  nonnegative  matrix  from  Radwan’s  Theorem. 

ii.  If  As  +  A3  +  Ae  ^  0,  then  applying  Kellogg’s  conditions  to  ai  we  have 
K  =  ^.  Since  As  +  A3  +  Ae  ^  0  it  follows  that  Ai  +  A3  +  Ae  >  0,  so  we 
may  construct  a  symmetric  nonnegative  matrix  B  having  spectrum  a 
from  two  smaller  3x3  symmetric  nonnegative  matrices  Ai  and  As 
where  cr  (Ai)  =  {Ai,  A4,  A5}  and  a  (As)  =  {As,  A3,  Ae}-  Here 

Ai  0 

B  = 

0  A 

has  spectrum  a. 

iii.  When  the  Perron  root  condition  is  not  satisfied,  that  is  Ai  +  A3  +  Ae  < 
0,  we  cannot  apply  Radwan’s  result.  This  means  the  existence  of  a 
symmetric  nonnegative  matrix  or  a  nonnegative  matrix  is  unknown. 

(g)  Note,  that  to  completely  analyze  all  partitions  we  must  consider  the  case 
where  Ai  ^  As  >  0  ^  A4  +  Ae  ^  A3  +  A5.  But  this  is  impossible  given  our 
initial  assumptions  that  A3  ^  A4  and  A5  ^  Ae- 
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(h)  Suppose  Ai  ^  A2  >  0  ^  A3  ^  A4  ^  As  +  As- 

i.  If  A2  +  As  +  Ae  <  0,  then  K  =  {2}  in  Kellogg’s  conditions.  Assuming 
Ai  +  As  +  As  ^  0  and  Ai  +  (A2  +  As  +  Ae)  +  A3  +  A4  ^  0,  there  exists  a 
symmetric  nonnegative  matrix  having  a  as  the  spectrum  by  Radwan’s 


Theorem. 

ii.  If  A2+AS+A6  ^  0,  then  K  =  %.  Since  A2+A5+A6  ^  0  =>  Ai+As+Ae  ^ 
0,  we  get  Ai  +  A3  +  A4  ^  Ai  +  A5  +  Ae  ^  0.  This  guarantees  the 
existence  of  two  symmetric  nonnegative  3x3  matrices,  Ai  and  A2 
where  cr(Ai)  =  {Ai,A3,A4}  and  (7(^2)  =  {A2,A5,A6}.  So  we  may 
construct  the  symmetric  nonnegative  matrix 


0  A2 


which  has  the  desired  spectrum  a. 

iii.  When  the  Perron  root  condition  is  not  satisfied,  that  is  Ai  +  A5  + 
Ae  <  0,  we  cannot  apply  Radwan’s  result.  This  means  the  existence 
of  a  symmetric  nonnegative  matrix  or  a  nonnegative  matrix  is  still 
unknown. 


(i)  Suppose  Ai  ^  A2  >  0  ^  A3  ^  A5  ^  A4  +  Ae. 

i.  If  A2  +  A4  +  Ae  <  0,  then  K  =  {2}.  And  by  Radwan’s  Theorem,  since 
Ai  +  A4  +  Ae  ^  0  by  assumption  and  Ai  +  (A2  +  A4  +  Ae)  +  A3  +  A5  ^  0, 
there  exists  a  symmetric  nonnegative  matrix  having  <7  as  the  spectrum. 
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ii.  If  Aa  +  A4  +  Ae  ^  0,  then  K  Assume  Ai  +  A4  +  Ag  ^  0.  Notice  also 
that  since  A3  +  A5  ^  A4  +  As,  we  have  Ai  +  A3  +  A5  ^  Ai  +  A4  +  As  ^ 


0.  The  conditions  give  the  existence  of  two  symmetric  nonnegative 
3x3  matrices,  Ai  and  A2  where  a(Ai)  =  {Ai,A3,A5}  and  cr(A2)  = 
{Aa,  A4,  Ae}.  So  we  may  construct  the  symmetric  nonnegative  matrix 


0  Aa 


which  has  the  desired  spectrum  a. 

iii.  When  the  Perron  root  condition  is  not  satisfied,  that  is  Ai  +  A4  + 
Ag  <  0,  we  cannot  apply  Radwan’s  result.  This  means  the  existence 
of  a  symmetric  nonnegative  matrix  or  a  nonnegative  matrix  is  still 
unknown. 


(j)  Suppose  Ai  ^  Aa  >  0  ^  A4  ^  A5  ^  A3  +  Ag. 

i.  If  Aa  +  As  +  Ag  <  0,  then  K  =  {2}.  Assume  Ai  +  A3  +  Ag  ^  0.  By 
Radwan’s  Theorem,  since  Ai  +  A3  +  Ag  ^0  and  Ai  +  (Aa  +  A3  +  Ag)  + 
A4  +  Ag  ^  0,  there  exists  a  symmetric  nonnegative  matrix  having  o  as 

the  spectrum. 

ii.  If  Aa  +  A3  +  Ag  ^  0,  then  AT  =  0.  Continuing  with  the  constraints  we 
check  if  Ai  +  A3  +  Ag  ^  0  and  Ai  +  A4  +  Ag  ^  0.  The  first  one  is  true 
because  Ai  +  A3  +  Ag  ^  Aa  +  A3  +  Ag  ^  0  but  the  second  one  has  no 
guarantee.  Now  we  must  consider  the  sign  of  Ai  +  A4  +  Ag. 
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If  A1+A4+A5  >0  then  we  have  two  symmetric  nonnegative  3x3  matri¬ 
ces,  Ai  and  A2  where  a  (Ai)  =  {Ai,  A4,  A5}  and  a  {A2)  =  {A2,  A3,  Ae}. 
So  we  may  construct  the  symmetric  nonnegative  matrix 


B  = 


0 


0  A2 


which  has  the  desired  spectrum  a. 

If  Ai  +  A4  -t-  A5  <0  we  cannot  apply  Radwan’s  result.  This  means  the 
existence  of  a  symmetric  nonnegative  matrix  or  a  nonnegative  matrix 
is  still  unknown. 

iii.  When  the  Perron  root  condition  is  not  satisfied,  that  is  Ai  -f  A3  -I- 
Ag  <  0,  we  cannot  apply  Radwan’s  result.  This  means  the  existence 
of  a  symmetric  nonnegative  matrix  or  a  nonnegative  matrix  is  still 
unknown. 


(k)  Suppose  Ai  ^  A2  >  0  ^  A3  ^  Ae  ^  A4  -1-  A5. 

i.  If  A2  +  A4  -I-  As  <  0,  then  K  =  {2}.  By  Radwan’s  Theorem,  since 

Ai  +  A4  -f  As  ^  0  by  assumption  and  Ai  -H  (A2  +  A4  -f  As)  +  A3  -I-  Ae  ^  0, 
there  exists  a  symmetric  nonnegative  matrix  having  a  as  the  spectrum. 

ii.  If  A2  +  A4  -I-  As  ^  0,  then  RT  =  0.  Continuing  with  the  constraints  we 

check  if  Ai  -f  A4  -f  As  ^  0  and  Ai  -f  A3  4-  Ae  ^  0.  The  first  one  is  true 

because  Ai  -1-  A4  -i-  As  ^  A2  +  A4  -1-  As  ^  0  but  the  second  one  has  no 

guarantee.  Now  we  must  consider  the  sign  of  Ai  -t-  A3  -f-  Ae- 
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If  Ai+As+Ae  ^  0  then  we  have  two  symmetric  nonnegative  3x3  matri¬ 
ces,  Ai  and  A2  where  a{Ai)  =  {Ai,A3,A6}  and  <7(^2)  =  {A2,A4,A5}. 
So  we  may  construct  the  symmetric  nonnegative  matrix 


B  = 


0 

0  A2 


which  has  the  desired  spectrum  a. 

If  Ai  -I-  A3  -I-  Ae  <  0  we  cannot  apply  Radwan’s  result.  This  means  the 
existence  of  a  symmetric  nonnegative  matrix  or  a  nonnegative  matrix 
is  still  unknown. 

iii.  When  the  Perron  root  condition  is  not  satisfied,  that  is  Ai  -1-  A4  -t- 
A5  <  0,  we  cannot  apply  Radwan’s  result.  This  means  the  existence 
of  a  symmetric  nonnegative  matrix  or  a  nonnegative  matrix  is  still 
unknown. 


(1)  Suppose  Ai  >  A2  >  0  ^  A3  ^  A4  -f  As  ^  As- 

i.  If  A2  -f-  Ae  <  0,  then  K  =  {2}.  Our  next  conditions  to  satisfy  are 
Ai  +  Ae  ^  0  and  Ai  +  (A2  +  Ag)  -!-  A3  A4  -H  A5  ^  0.  Both  of  these 
are  true.  The  first  is  because  of  the  Perron  root  condition  on  <7 .  The 
second  is  the  nonnegative  trace  condition.  Here  Radwan’s  Theorem 
applies  and  we  have  the  existence  of  a  symmetric  nonnegative  matrix 
having  o  as  the  spectrum. 
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ii.  If  A2  +  As  >  0,  then  K  =  0.  Checking  the  next  condition  of  Kellogg, 
we  see  Ai  +  A3  +  A4  +  A5  ^  0  because  of  (3). 

Hence,  this  particular  partition  gains  no  further  restrictions  for  the  exis¬ 
tence  of  symmetric  nonnegative  matrices  or  nonnegative  matrices.  We  do 
however  point  out  that  no  other  partition  Ai  ^  A2  >  0  ^  Aj  ^  A^  +  A*,  ^  Ae 
works  for  i,j,  k  being  nonrepeating  elements  from  {3, 4, 5}. 

(m)  Suppose  Ai  ^  A2  >  0  ^  A3  -f-  A4  ^  A5  ^  Ag. 

i.  If  A2  -f  Ae  <  0,  then  K  =  {2}.  Our  next  conditions  to  satisfy  are 
AiH-Ag  ^  0  and  Ai-4-(A2  +  A6)-f-A3-f  A4-1-A5  ^  0.  Both  of  these  are  true. 
The  first  is  because  of  the  Perron  root  condition  on  a.  The  second  is 
the  nonnegative  trace  condition.  Here  Radwan’s  result  applies  and  we 
have  the  existence  of  a  symmetric  nonnegative  matrix  having  a  as  the 
spectrum. 

ii.  If  A2  -t-  Ae  ^  0,  then  X  =  0.  Checking  the  next  condition  of  Kellogg, 
we  see  Ai  -h  A3  -I-  A4  -I-  A5  ^  0  because  of  (3).  Hence,  this  particular 
partition  gains  no  further  restrictions  for  the  existence  of  symmetric 
nonnegative  matrices  or  nonnegative  matrices. 

7.1.3  Ai  ^  A2  ^  A3  >  0  ^  A4  ^  A5  ^  Ae 

In  this  particular  distribution  of  numbers  we  find  many  more  cases  because  we 
have  more  positive  values  to  match  up  with  negative  ones.  For  example,  the  straight 
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application  of  Kellogg’s  conditions  requires  us  to  analyze  A2  +  Ae  and  A3  +  A5  to 
determine  whether  or  not  2,3  E  K .  We  will  consider  all  possible  outcomes  through 
the  following. 

1.  If 

A2  d"  Ag  <  0  (4) 

and 

A3  +  A5  <  0,  (5) 

then  K  =  {2,3}.  We  will  examine  the  different  outcomes  of  the  remaining 
conditions  of  Kellogg,  that  is  Ai  +  Ag  ^  0,  Ai  +  (A2  +  Ag)  +  Ag  ^  0,  and  Ai  + 
+  Ag)  +  {A3  +  As)  +  A4  ^  0.  The  first  and  third  are  true  by  the  necessary 
conditions  on  g  with  regard  to  the  Perron  root  condition  and  the  nonnegative 
trace  condition.  The  second  condition  need  not  be  true. 

(a)  If  Ai  +  (A2  +  Ag)  +  As  ^  0,  then  as  we  satisfy  all  of  Kellogg’s  conditions, 
we  have  by  Fiedler’s  result  a  symmetric  nonnegative  matrix. 

(b)  If 

Ai  +  (A2  +  Ag)  +  A5  <  O5  (6) 

then  Kellogg’s  conditions  are  not  satisfied.  Note  that  Ai+(A2  +  A6)+A5  <  0 
gives 

A3  +  A4  ^  0  C^) 
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due  to  the  nonnegative  trace  requirement.  In  this  case  where  Kellogg  s 
conditions  are  not  met,  we  try  partitioning  to  use  Borobia’s  and  Radwan’s 
results  to  further  resolve  the  problem. 

i.  If  Ai  ^  A2  ^  -^3  >  0  ^  A4  ^  As  +  Ag.  Then  to  make  this  partition 
viable,  we  must  have  Ai  +  A5  +  As  ^  0.  But  this  clearly  contradicts  (6) 
as  A2  >  0.  Therefore  this  partition  fails  to  give  new  information. 

ii.  If  Ai  >  A2  >  As  >  0  ^  As  ^  A4  +  As,  then  the  following  cases  analyze 
the  possible  K  sets  for  Borobia’s  conditions. 

A.  If  A2  +  A4  +  Ae  <  0  and  A3  +  A5  <  0,  then  K  =  {2,3}.  Then 
Borobia’s  and  Radwan’s  conditions  are  satisfied  if  Ai  +  A4  +  Ae  ^  0 
and  Ai  +  (A2  +  A4  +  Ae)  +  As  ^  0.  The  condition  Ai  +  A4  +  Ae  ^  0 
may  or  may  not  be  true.  For  this  case  we  will  assume  it  is,  but  we 
will  address  when  it  does  not  in  the  final  case  of  this  partition.  But 
the  second  condition  Ai  +  (A2  +  A4  4-  Ae)  +  A5  ^  0  cannot  be  true 
as  it  contradicts  (6).  Therefore  this  partition  when  A  =  {2,3} 
cannot  happen. 

B.  If  A2  +  A4  +  A6  <  Oand  A3  +  A5  ^  0,  then  AT  =  {2}.  But  A3  +  A3  ^  0 
contradicts  (5).  Therefore  this  partition  cannot  happen. 

C.  If  A2+A4+Ae  ^  0  and  A3+A5  <  0,  then  AT  =  {3}.  But  A2+A4+Ae  ^ 
0  contradicts  (4).  Therefore  this  partition  cannot  happen. 

D.  If  A2  +  A4  +  Ae  ^  0  and  A3  +  A5  ^  0,  then  AT  =  0.  But  again,  each 
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condition  contradicts  either  (4)  or  (5).  Therefore  this  partition 
cannot  happen. 

E.  If  Ai  +  A4  +  Ae  <  0,  then  the  Perron  root  condition  is  not  satisfied 
in  this  partition.  Then  this  partition  does  not  resolve  whether  or 
not  there  exists  a  symmetric  nonnegative  matrix  or  a  nonnegative 
matrix  having  a  as  a  spectrum. 

iii.  If  Ai  >  A2  ^  A3  >  0  ^  As  ^  A4  +  As,  then  the  following  cases  analyze 
the  possible  K  sets  for  Borobia’s  conditions. 

A.  If  A2  +  A4  +  As  <0  and  A3  +  As  <  0,  then  K  =  {2,3}.  Then 
Borobia’s  and  Radwan’s  conditions  are  satisfied  if  Ai  +  A4  +  As  ^  0 
and  Ai  +  (A2  +  A4  +  As)  +  As  ^  0.  The  condition  Ai  +  A4  +  As  ^  0 
mav  or  may  not  be  true.  For  this  case  we  will  assume  it  is,  but  we 
will  address  when  it  does  not  in  the  final  case  of  this  partition.  But 
the  second  condition  Ai  +  (A2  +  A4  +  As)  +  As  ^  0  cannot  be  true 
as  it  contradicts  (6).  Therefore  this  partition  when  K  =  {2,3} 
cannot  happen. 

B.  If  A2  +  A4  +  A5  <  0  and  As  +  As  ^  0,  then  K  =  {2}.  But  As  +  As  ^  0 
contradicts  (5).  Therefore  this  partition  cannot  happen. 

C.  If  A2+A4+AS  ^  0  and  As+As  <  0,  then  K  =  {3}.  But  A2+A4+AS  ^ 
0  contradicts  0  >  A2  +  As  ^  A24-A4  +  A5  using  (4)  and  the  partition 
(iii).  Therefore  this  partition  cannot  happen. 
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D.  If  A2  +  A4  +  As  ^  0  and  A3  +  Ae  ^  0,  then  X  =  0.  But  again, 
each  condition  contradicts  known  prior  conditions.  Therefore  this 
partition  cannot  happen. 

E.  If  Ai  +  A4  +  As  <  0,  then  the  Perron  root  condition  is  not  satisfied 
with  this  partition.  Then  this  partition  does  not  resolve  whether  or 
not  there  exists  a  symmetric  nonnegative  matrix  or  a  nonnegative 
matrix  having  a  as  a  spectrum. 

iv.  If  Ai  ^  A2  ^  A3  >  0  ^  A4  +  As  ^  Ae,  then  we  show  that  the  only 
viable  case  for  K  is  when  K  =  {2,3}.  This  means  the  inequalities 
A2  +  Ae  <  0  and  A3  +  A4  +  As  <  0  are  true.  We  know  they  must  be 
because  of  (4)  and  (5).  These  conditions  rule  out  any  other  possible 
K.  So  in  this  only  viable  case,  for  the  remaining  conditions  to  be  true 
we  need  Ai  +  Ag  >  0  and  Ai  +  (A2  +  Xe)  +  A4  +  As  ^  0.  Unfortunately 
the  second  one  contradicts  (6).  Therefore  even  this  case  fails. 

V.  If  Ai  ^  A2  ^  A3  >  0  ^  A4  +  As  +  Ae,  then  to  discuss  Borobia’s  and 
Rad  wan’s  conditions  we  must  have  Ai  +  A4  +  As  +  As  ^  0.  But  clearly 
this  contradicts  (6),  so  this  partition  gains  us  nothing. 
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The  previous  partitions  were  based  on  the  original  distribution  of  ^ 
A2  ^  A3  >  0  ^  A4  ^  Ao  ^  Ae  where 

A2  +  Ae  <  0, 

^^3  +  -^5  <  0) 

Ai  +  (A2  +  Ag)  +  Ag  <  0,  and 

A3  +  A4  ^  0. 

2.  Let  us  now  address  the  case  w'here  we  have  K  =  {2}  with  regard  to  the  original 
spread  of  Ai  ^  A2  ^  A3  >  0  ^  A4  ^  Ag  ^  Ag.  If 

A2  +  Ag  <  0  (8) 

and 

A3  +  Ag  ^  0,  (9) 

then  K  =  {2}.  We  will  examine  the  different  outcomes  of  the  remaining  condi¬ 
tions  of  Kellogg,  that  is  Ai  -t-  Ag  ^  0,  and  Ai  +  (A2  +  Ag)  -f  A4  ^  0.  The  first  is 
true  by  the  necessary  condition  on  a  with  regard  to  the  Perron  root  condition. 
The  second  condition  need  not  be  true.  We  will  analyze  this  in  the  following. 

(a)  If  Ai  -b  Ag  ^  0  and  Ai  +  (A2  +  Ag)  +  A4  ^  0,  then  Kellogg’s  Theorem  ap¬ 
plies  and  by  Fiedler’s  result  we  know  there  exists  a  symmetric  nonnegative 
matrix  having  a  as  the  spectrum. 

(b)  If 

Ai  +  (A2  +  Ag)  -f  A4  <  0,  (10) 
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Kellogg’s  conditions  fail.  This  leaves  us  with  Borobia’s  and  Radwan’s 
results  which  allow  for  partitions.  The  following  cases  consider  all  possible 

partitions. 

i.  If  Ai  ^  A2  ^  As  >  0  ^  As  ^  A4  +  As,  then  consider 

=  {Ai,  A2,  As,  As,  A4  +  As}.  We  will  analyze  the  possible  K  sets  from 
applying  the  conditions  to  ai . 

A.  If  A2  +  A4  +  As  <  0  and  As  +  As  <  0,  then  K  =  {2,3}.  Combine  this 
with  the  Perron  root  condition  by  assuming  that  Ai  +  A4  +  As  ^  0, 
then  we  need  only  consider  the  inequality  Ai+(A2  +  A4  +  As)+A6  ^ 
0.  Unfortunately  this  final  inequality  contradicts  (10).  Therefore 
this  case  cannot  happen. 

B.  If  A2  +  A4  +  AS  <  Oand  As  +  Ae  >  0,  then  K  =  {2}.  But  As  +  Ae  ^  0 
contradicts  (8),  so  this  case  cannot  happen. 

C.  If  A2  +  A4  +  As  ^  0  then  using  (8)  and  the  partition  we  get  0  > 
^2  +  As  ^  A2  +  A4  +  As  which  is  a  contradiction.  So  we  have  ruled 
out  K  =  {3}  and  K  =  0. 

D.  If  A1+A4+ As  <  0,  then  we  do  not  satisfy  the  Perron  root  condition 
with  this  partition.  Therefore  we  do  not  know  whether  or  not  there 
exists  a  symmetric  nonnegative  matrix  or  a  nonnegative  matrix 
having  cr  as  a  spectrum. 
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ii.  If  Ai  ^  As  ^  As  >  0  ^  A4  ^  As  +  Ae,  then  consider 

(Ti  -  {Ai,  As,  As,  A4,  As  +  As}.  But  this  set  fails  the  Perron  root  condi¬ 
tion,  that  is  Ai  -l-  As  +  As  ^  0  because  it  contradicts  (10).  Therefore, 
this  partition  does  not  give  us  anything  new. 

iii.  If  Ai  ^  As  ^  As  >  0  ^  As  ^  A4-f  As.  Then  consider  d  =  {Ai,  As,  As,  As,  A4-t- 
Ae}.  We  will  analyze  the  possible  K  sets  from  applying  Kellogg’s  con¬ 
ditions  to  (Ti. 

A.  If  A3  +  A4  +  A6  <  0  and  As  +  As  <  0,  then  K  =  {2,3}.  Combine  this 
with  the  Perron  root  condition  by  assuming  that  Ai  +  A4  -t-  Ae  ^  0, 
then  we  need  only  consider  the  inequality  Ai-l-(As  +  X4  +  A6)+As  ^ 

0.  Unfortunately  this  final  inequality  contradicts  (10).  Therefore 
this  case  cannot  happen. 

B.  If  As  -h  A4  -t-  Ae  <  0  and  A3  +  A5  ^  0,  then  K  =  {2}.  But  to  apply 
Borobia’s  and  Radwan’s  conditions,  we  need  Ai  -1-  A4  +  Ae  ^  0. 
This  inequality  contradicts  (10)  so  this  case  cannot  occur. 

C.  If  AS  +  A4-HA6  ^  0,  then  we  contradict  (8).  It  follows  that  K  =  {3} 
and  =  0  cannot  occur. 

D.  If  Ai  A4-f- Ae  <  0,  then  we  do  not  satisfy  the  Perron  root  condition 
with  this  partition.  Therefore  we  do  not  know  whether  or  not  there 
exists  a  symmetric  nonnegative  matrix  or  a  nonnegative  matrix 
having  cr  as  a  spectrum. 
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iv.  If  Ai  ^  A2  ^  A3  >  0  ^  A4  +  As  ^  Ae,  then  consider 


=  {Ai,  A2,  A3,  A4  +  As,  Ae}.  We  will  analyze  the  possible  K  sets  from 
applying  Kellogg’s  conditions  to  ai. 

A.  If  A2  +  Ag  <  0  and  A3  +  A4  +  As  <  0,  then  K  =  {2, 3}.  Here  the 
necessary  Perron  root  condition  is  met  since  that  Ai  +  Ag  ^  0  for 
a.  Hence  we  need  only  consider  the  inequality  Ai  +  (A2  +  Ag)  + 
A4  +  As  ^  0.  Unfortunately  this  final  inequality  contradicts  (10). 
Therefore  this  case  cannot  happen. 

B.  If  A2  +  Ag  <  0  and  A3  +  A4  +  As  ^  0,  then  K  =  {2}.  Then 
we  may  subdivide  ai  into  <72  =  {Ai,Ag},  <73  =  {A3,A4,As},  and 
0-4  =  {A2}.  Each  of  these  lists  satisfy  Suleimanova’s  condition  to 
be  the  spectrums  of  symmetric  nonnegative  matrices,  A2,  A3,  and 
A4  respectively.  Therefore  a  symmetric  nonnegative  matrix  having 
a  as  its  spectrum  is 


C.  If  A2  +  Ag  ^  0  it  contradicts  (8).  So  it  follows  that  K  =  {3}  and 
K  =  0  cannot  occur. 

D.  If  Ai  ^  A2  ^  A3  >  0  ^  A4  +  As  +  Ag.  Then  the  only  possible  group¬ 
ing  under  question  is  A2  +  A4  +  As  +  Ag.  As  we  must  satisfy  (10), 
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this  expression  must  be  negative,  that  is  A2  +  A4  +  A5  +  Ae  <  0. 
Therefore  the  only  element  in  the  K  set  for  this  partition  is  2.  As 
+  A4  +  A5  +  Ae  <  0  because  of  (10),  we  find  the  Perron  root 
condition  fails  for  this  partition.  Therefore  this  partition  gives  no 
new  information  to  the  problem. 

The  previous  partitions  were  based  on  the  original  distribution  of  six 
real  numbers  where 

Ai  ^  A2  >  A3  >  0  ^  A4  ^  Ae  >  ^6, 

A2  +  ^6  0) 

A3  +  A5  ^  0, 

Ai  +  (A2  +  Ae)  +  A4  <  0,  and 
A3  +  As  ^  0. 

3.  Let  us  now  address  the  case  where  we  have  K  =  {3}  with  regard  to  the  original 
spread  of  Ai  ^  A2  ^  A3  >  0  ^  A4  ^  Ae  ^  Ae-  If 

A2  +  Ae  ^  0  ^ 


and 


A3  +  Ae  <  0, 


(12) 


then  K  =  {2>}.  We  will  examine  the  different  outcomes  of  the  remaining  con- 
ditions  of  Kellogg,  that  is  Ai  +  A,  S  0,  and  A,  +  (As  +  As)  +  A,  >  0.  The  first 
is  true  by  the  Perron  root  condition  on  o.  The  second  condition  need  not  be 

true.  We  will  analyze  this  in  the  following. 
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(a)  If  Ai  +  As  S  0  and  A,  +  (A3  +  As)  +  A4  »  0,  then  Kellogg’s  Theorem  ap¬ 
plies  and  by  Fiedler’s  result  we  know  there  exists  a  symmetric  nonnegative 
matrix  having  a  as  the  spectrum. 


(b)  If 


Ai  +.  (A3  +  As)  +  A4  <  0, 


(13) 


then  Kellogg’s  conditions  fail.  This  leaves  us  with  Borobia’s  and  Radwan’s 


results  which  require  partitions.  The  following  cases  consider  all  possible 
partitions. 


i.  If  Ai  ^  A2  ^  Ag  >  0  ^  A4  ^  As  +  Ae.  This  partition  does  not  apply  as 
the  Perron  root  condition  fails,  that  is  Ai  +  (As  +  As)  ^  0  by  (13). 


ii.  If  Ai  ^  Aa  ^  A3  >  0  ^  As  ^  A4  +  Ae.  This  partition  does  not  apply  as 
the  Perron  root  condition  fails,  that  is  Ai  +  (A4  +  Ae)  ^  0  by  (13). 

iii.  If  Ai  ^  Aa  ^  A3  >  0  >  Ae  ^  A4  +  As-  This  partition  does  not  apply  as 
the  Perron  root  condition  fails,  that  is  Ai  +  (A4  4-  As)  0  by  (13). 

iv.  If  Ai  ^  Aa  ^  A3  >  0  ^  A4  +  As  ^  Ae-  This  partition  cannot  occur 
because  Ai  +  A4  +  As  ^  Ai  +  Ae  ^  0  is  a  contradiction  to  (13). 

V.  If  Ai  ^  Aa  ^  A3  >  0  ^  A4  +  As  +  Ae-  This  partition  cannot  occur  given 
the  restriction  from  (13).  That  is  0  >  A1+A2+A4+AS  >  A1+A4+A5+A6 
fails  the  Perron  root  condition. 
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The  previous  partitions  were  based  on  the  original  distribution  of  six  number 
where 

Ai  >  A2  ^  As  >  0  ^  A4  ^  As  ^  ^^6’ 

A2  +  As  ^  0) 

A3  +  As  <  0,  and 
Ai  +  (A3  +  As)  +  A4  <  0. 

4.  Let  us  now  address  the  case  where  we  have  X  =  0  with  regard  to  the  original 
spread  of  A^  ^  A2  ^  A3  >  0  >  A4  ^  As  >  Ae-  For  K  =  0  then  we  must  have  A2  + 

As  >  0,  and  As+As  ^  0.  But  since  A4+A4  ^  0  we  satisfy  the  sufficient  conditions 
necessary  for  there  to  exist  three  2  x  2  symmetric  nonnegative  matrices,  A4,  A2, 
and  A3  having  spectrums  <7  (Ai)  =  {Ai,A4},  cr  (A2)  -  {A2,A6},  and  a  (A3) 

{A3,  As}.  So  we  may  construct  the  symmetric  nonnegative  matrix 

Ai  0  0 

B  =  0  A2  0 

0  0  A3 

realizing  o  as  its  spectrum. 

7.1.4  Ai  ^  A2  ^  A3  ^  A4  >  0  >  As  ^  As 

In  this  particular  distribution  the  straight  application  of  Kellogg’s  conditions  re¬ 
quires  us  to  analyze  As  -I-  Ae  and  A,  -h  A,  to  determine  whether  or  not  2, 3  £  K.  We 
will  consider  all  possible  outcomes  through  the  following. 
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A2  +  Ae  <  0 


(14) 


1.  If 


and 

As  +  A5  <  0,  (^^) 

then  K  =  {2, 3}.  We  will  examine  the  different  outcomes  of  the  remaining 

conditions  of  Kellogg,  that  is  Ai  +  Ae  ^  0,  Ai  +  (A2  +  Ae)  +  A5  ^  0,  and  Ai  + 

+  Ag)  +  (As  +  As)  ^  0.  The  first  is  true  by  the  necessary  condition  on  cr  with 
regard  to  the  Perron  root  condition.  The  second  condition  need  not  be  true.  The 
third  condition  depends  on  the  second.  Notice  that  if  Ai  +  (A2  +  Ag)  +  As  ^  0, 
then  Ai  +  (A2  +  Ag)  +  (As  +  A5)  <  0  cannot  happen  for  As  ^  0.  So  we  may 
conclude  the  only  viable  condition  is  the  second  one,  that  is  Ai  +  (A2  +  Ag)+A5  ^ 
0.  If  this  condition  is  satisfied  there  is  a  symmetric  nonnegative  matrix  whose 

spectrum  is  a. 

Suppose 

Ai  +  (A2  +  Ag)  +  As  <  0,  (I®) 

we  find  that  a  does  not  satisfy  the  Kellogg  conditions.  Let  us  consider  the 
following  partition  to  see  if  this  distribution  may  be  solved  using  Borobia’s  and 
Radwan’s  conditions.  If  Ai  >  A2  ^  A3  ^  A4  >  0  ^  As  +  Ag.  This  partition  does 
not  apply  as  the  Perron  root  condition  fails,  that  is  Ai  +  As  +  Ag  <  0  because 

of  (16). 


no 


A2  +  Ae  <  0 


(17) 


A3  +  A5  ^  0, 


then  K  ~  {2}.  Again,  using  the  Perron  root  condition  from  a  and  (16)  we  can 
find  symmetric  nonnegative  matrices,  Ai  and  A2  such  that  a  (Ai)  =  {Ai,  Ae}  and 
a{A2)  =  {A3,  As}.  Using  these,  we  can  construct  the  symmetric  nonnegative 


matrix 


Ai  0  0  0 

0  A2  0  0 

0  0  A2  0 

0  0  0  A4 


having  spectrum  a. 


A2  +  Ae  ^  0 


A3  +  A5  <  0, 


then  K  -  {3}.  Again,  using  the  Perron  root  condition  from  cr  and  (19)  we  can 
find  symmetric  nonnegative  matrices,  Ai  and  A2  such  that  a  (Ai)  =  (Ai,  As}  and 
a{A2)  =  {A2,A6}.  Using  these,  we  can  construct  the  symmetric  nonnegative 
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matrix 


having  spectrum  o. 


Ai  0  0  0 


B  = 


0  A2  0 
0  0  A3 


0 

0 


0  0  0  A4 


4.  If 


A2  +  Ae  ^  0 


(21) 


and 


A3  +  A5  ^  0, 


(22) 


then  K  =  %  Again,  using  (21)  and  (22)  we  can  find  symmetric  nonnegative  ma¬ 
trices,  Ai  and  A2  such  that  a  (Ai)  =  {A2,  Ae}  and  a  (A2)  =  {A3,  As}-  Therefore, 


we  can  construct  the  symmetric  nonnegative  matrix 

Ai  0  0  0 

0  A2  0  0 

B  = 

0  0  Ai  0 

0  0  0  A4 


having  spectrum  a. 


7.1.5  Ai  ^  A2  ^  A3  ^  A4  ^  As  >  0  ^  Ag 

In  this  particular  distribution  we  find  a  simple  solution  from  the  necessary  Perron 
root  condition,  that  is  A,+As  >  0.  With  A,,  A,  as  the  eigenvalues  of  a  2  x  2  symmetric 
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nonnegative  matrix,  we  can  easily  find  a  symmetric  nonnegative  matrix  having  a  as 


its  spectrum.  The  following  is  one  such  matrix. 

-^1  +-^6 

2 

Ai— Ae 

2 

0 

0 

0 

0 

Ai— Ae 

2 

Ai  -f-Ae 

2 

0 

0 

0 

0 

0 

0 

A2 

0 

0 

0 

0 

0 

0 

A3 

0 

0 

0 

0 

0 

0 

A4 

0 

0 

0 

0 

0 

0 

A5 

7.2  Additional  Results 

In  chapter  6  we  have  a  useful  result  from  McDonald  and  Neumann  [21]  in  regard 
to  an  irreducible  5x5  nonnegative  symmetric  matrix  having  eigenvalues  1  ^  A2  ^ 
As  ^  A4  ^  A5  ^  —1  where  1  is  the  Perron  root.  A  similar  result  is  obtained  in  higher 
dimensions.  We  give  a  proof  in  the  same  lines  of  McDonald  and  Newmann  for  a  list 
of  six  real  numbers  and  then  generalize  it  to  a  list  of  n  real  numbers. 
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Theorem  31.  Let  A  be  a  6  x  6  irreducible  nonnegative  symmetric  matrix  with  eigen¬ 
values  {1,A2,A3,A4,A5,A6}  where  1  ^  A2  ^  A3  >  A4  ^  As  ^  As  ^  -1-  Then 
-2Xq  +  trace  (A)  ^  A2. 

Moreover,  if  some  eigenvector  of  A  corresponding  to  X2  has 

1.  exactly  two  positive  or  exactly  two  negative  entries,  then  -Xe  +  trace{A)  >  A2 
or 

2.  exactly  one  positive  or  exactly  one  negative  entry,  then  trace  (A)  ^  A2. 

Proof.  Let  A  =  K]  where  i,j  =  1,2, 3, 4, 5, 6,  be  an  irreducible  nonnegative  sym¬ 
metric  matrix  with  eigenvalues  {1,  A2,  A3,  A4,  A5,  Ae}  where  1  ^  A2  ^  A3  >  A4  ^ 

6 

As  ^  Ae  ^  -1-  First  suppose  Ae  ^  0.  Then  -As  -t- trace  (A)  -  -As  +  ^^A^ 

1  +  A2  -t-  A3  4-  A4  +  As  ^  A2,  since  A*  ^  0  for  t  =  1, 2, . . . ,  6. 

Now  suppose  Ae  <  0.  Then  let  B  =  A  -  X,L  B  is  positive  semidefinite  since  B 

is  symmetric  and  has  nonnegative  eigenvalues,  that  is  a  (B)  =  {1  -  Ae,  A2  Ae,  A3 
Ae,  A4  -  Ae,  As  -  Ae,  0}.  And  since  A  is  nonnegative,  we  find  that  A  -  Ael  =  5  is  also 
nonnegative.  Therefore  B  is  positive  semidefinite,  nonnegative,  and  symmetric. 

Let  X  be  the  Perron  vector  corresponding  to  the  Perron  root  1  -  Ae  for  B.  By 
the  Perron-Frobenius  Theorem  x  >  0.  So  we  may  choose  x  such  that  Hxlb  =  L  Let 
y  be  an  eigenvector  corresponding  to  A2  -  Ae-  Since  B  is  symmetric,  x  y  -  0  so  we 
may  conclude  that  y  has  at  least  one  negative  entry  and  at  least  one  positive  entry. 
As  y  has  six  entries,  then  either  y  or  -y  has  at  most  three  positive  entries.  Without 
loss  of  generality  let  y  contain  at  most  three  positive  entries  positioned  at  the  top 
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of  the  vector.  Let  Vi  be  the  positive  part  of  y,  and  y2  be  the  nonpositive  part  of  y. 
Therefore 

= 


yi 


2/2 

Partition  B  to  match  y.  So  when  yi  is  3  x  1  then 

yi  Buyi  +  Bi2y2 

y2  B2iyi  +  ^222/2 


Bn  B\2 
B21  B22 


where  B\i  is  3  x  3  and  B22  is  3  x  3. 

Then  Buy,  +  BnV2  =  (Aj  -  A,)  y,  and  since  BnV2  <  0  it  follows  that  B^^y,  > 
_  x,)yi.  Now  as  B  is  nonnegative  and  semidefinite,  so  is  Bn,  therefore  as 


we  have 


On  —  Ag 

012 

Ol3 

Bn  = 

021 

022  Ae 

O23 

031 

032 

O33  ~  Ai 

trace  (Bn) 

=  ail  +  CI22  +  ^33  “ 

-3A6 

<  Oil  +  022  +  O33  +  O44  +  <^55  “  SAe 
=  trace  (A)  —  3A6. 

Now  the  Perron  root  p  (Bn)  satisfies  p  (Bn)  ^  A2  -  Ag.  Therefore  -3A6 +troce  (A)  ^ 
trace  {Bn)  >  P  (Bn)  >  A2  -  Ae  giving  our  desired  conclusion  -2A6  +  trace  (A)  ^  A2. 

A  similar  process  applies  when  yi  is  2  x  1  or  1  x  1.  In  the  first  case,  this  implies 
Bn  is  a  2  X  2  matrix  and  B22  is  a  4  x  4  matrix,  so  Bnyi  >  (A2  -  \)  Vi  so  we  have 
-2A6  +  trace  (A)  ^  trace  (Bn)  ^  P  (Bn)  ^  A2  -  Ae  giving  -  As  +  trace  (A)  ^  A2. 
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When  yi  is  1  X  1,  then  —As  +  trace  {A)  ^  trace  {Bn)  ^  p{Bii)  ^  A2  Ag. 
Therefore,  trace  (A)  ^  A2.  ^ 

An  easy  generalization  of  the  previous  result  is 

Theorem  32.  Let  Abe  annxn  irreducible  nonnegative  symmetric  matrix  with  eigen¬ 
values  {1,  A2, . . . ,  A„}  where  1  >  A2  ^  ^  A„  ^  -1,  then  (4  -  n)  A„  +  trace  (A)  > 

A2. 

Another  simple  result  uses  the  ideas  from  Perfect  and  Reams  [29].  Let  Ai,  A2,  •  ■  - ,  Ag 
be  the  eigenvalues  of  a  6  x  6  matrix,  A.  We  may  represent  by  Pi  the  sum  of  all 

products  of  the  eigenvalues  as  follows: 

6 

1.  Pi  =  Z)  Ai, 

1=1 

2.  P2  =  A1A2  +  A1A3  +  A1A4  +  A1A5  +  AiAg  +  A2A3  +  A2A4  +  A2A5  +  A2A6  +  A3A4  + 
A3A5  A3A6  “1“  A4A5  +  A4A5  “h  AsAgj 

3.  ps  =  A1A2A3  +  •  •  *  +  A4A5A6, 

4.  p4  =  A1A2A3A4  +  •  •  •  +  A3A4A5A6, 

5.  P5  =  A1A2A3A4A5  +  A1A2A3A4A6  +  A1A2A3A5A6  +  A1A2A4A5A6  +  A1A3A4A5A6  + 

A2A3A4A5A6,  and 

6.  pe  =  A1A2A3A4A5A6. 

Notice  that  pi  is  the  trace  of  A,  pg  is  the  determinant  of  A,  and  the  characteristic 
polynomial  equation  of  A  is  A®  —  piA®  +  P2A‘^  —  P3A®  +  P4A^  —  psA  +  pg  —  0. 
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Recall  Newton’s  identities  for  symmetric  functions  [11],  state 

Si  -  Pi  =  0, 

S2  —  Pl^l  +  2p2  —  Oj 

33  -P1S2  +  P2S1  -  3p3  =  0, 

S4  —  P1S3  +  P2S2  —  P3S1  +  P4  =  0, 

S5  —  P1S4  +P2S3  ~  P3S2  +P4S1  —  Ps  =  0) 

S6  -  P1S5  +  P2S4  -  P3S3  +  P4S2  -  P5S1  +  P6  =  0. 

Now  solving  for  the  coefficients  in  terms  of  the  moments  when  the  trace  of  A  is  zero^ 

that  is  Pi  =  0,  we  find: 

Pi  =  Si  =  0 

P2  =  -^52 

P3  =  ^^3 

2pl  -  S4 


4 


S2  2s4 

""  8 

— P3S2  +  P2013  +  S5 

P5  =  5 

^  +  35 

^  g 

P3S3  -  Se  —  P2S4  —  P4S2 
P6  =  - - g 

e-  5?— 2S4 

^S3  —  Sg  +  Y'^4  ~  8 

"  -  g 

8S3  —  Ss®  —  24s6  +  I8S2S4 

"  144  ” 
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Using  the  above  with  the  knowledge  that  pe  is  the  determinant  of  A  we  have  the 
following  special  case  of  the  JLL  condition. 


Theorem  33.  Let  A  be  a  singular,  trace  zero,  6x6  matrix  with  spectrum  a  - 
{Ai,  A2,  A3,  A4,  A5,  Ag}.  Then  3sq  ^  s^. 

Proof.  As  A  is  singular  we  know  pg  =  0-  Therefore, 

— 245g  A  I8S2S4  +  8S3  3s2  0 

3s2  (654  S2)  —  24sg  853 

but  S2  ^  0  and  the  JLL  inequality  gives  6S4  -  S2^  ^  0,  so 

24s6  -Ssl^O 

3sg  ^  S3. 


□ 


7.3  Conclusion 

This  chapter  analyzes  the  possible  outcomes  for  lists  of  six  real  numbers  in  order 
to  solve  the  SNIEP  and  the  RNIEP,  and  resolves  65  out  of  a  possible  79  cases  for  the 

SNIEP. 

We  now  turn  to  the  final  chapter  of  new  results  regarding  the  NIEP  where  we  have 
added  a  twist  to  the  problem.  Here  we  address  whether  or  not  a  list  of  n  numbers 
can  be  the  spectrum  of  a  nonnegative  matrix  that  is  subordinate  to  a  bipartite  graph. 
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CHAPTER  8 


NIEP  for  Matrices  Subordinate  to  a  Given  Graph 

8.1  Graph  Version  of  the  NIEP 

Although  most  of  this  dissertation  has  addressed  the  NIEP  in  a  general  setting 
mathematically  it  is  natural  to  consider  a  graph  version  of  the  NIEP.  Recall  a  graph, 
G,  is  a  set  of  vertices,  {yi,V2,  ■ .  • ,  Vn}  and  a  set  of  edges,  {ei,  62, ... ,  Cm}  that  connect 
vertices.  We  will  consider  simple  graphs  with  no  loops  or  multiple  edges.  A  directed 
graph  is  one  in  which  there  is  a  direction  associated  with  each  edge.  The  graph 
theoretic  version  of  the  NIEP  requires  that  if  two  vertices  do  not  have  an  edge  between 
them  in  a  given  graph  G,  then  the  matrix  realizing  tr  =  {Ai, . . . ,  A„}  must  have  a 
zero  entry  in  that  position.  So  for  a  directed  graph  where  the  edge  has  direction  from 
Vk  to  Vs,  the  realizing  matrix  A  =  (ay)  must  have  a  zero  entry  in  ask-  If  fhe  graph 
is  undirected,  then  a  realizing  matrix  A  =  (ay)  must  have  zeroes  located  in  the  aks 
and  ask  entries  when  Vk  and  Vs  are  not  joined  by  an  edge.  We  say  an  n  x  n  matrix 
A  —  [aij)  is  subordinate  to  a  graph  G  if  aij  ^  0  for  i  ^  j  implies  that  there  is  an  edge 
in  G  connecting  vertex  Vi  with  Vj.  Note  that  the  diagonal  entries  are  unrestricted. 

The  Graph  Nonnegative  Inverse  Eigenvalue  Problem  (G-NIEP)  relative  to  a  graph 
G  on  n  vertices  asks  to  identify  the  necessary  and  sufficient  conditions  for  lists  of  n 
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complex  numbers  to  occur  as  the  eigenvalues  of  a  nonnegative  matrix  subordinate  to 
the  given  graph  G.  Similarly  we  define  the  G-RNIEP  and  the  G-SNIEP.  Notice  that 
in  the  G-SNIEP  the  graph  G  must  be  undirected  and  the  list  of  n  numbers  must  be 
real  as  the  realizing  matrices  are  symmetric.  Leal-Duarte  and  Johnson  point  out  in 
their  paper  [19]  that  requiring  the  realizing  matrix  to  be  subordinate  to  G  rather  than 
having  the  graph  of  the  realizing  matrix  be  exactly  G,  guarantees  that  the  solution 

set  to  each  problem  is  always  a  closed  set. 

When  we  select  the  type  of  graph  we  open  a  gamut  of  new  problems.  In  our  case 
we  will  restrict  ourselves  to  the  G-NIEP  when  the  graph  in  question  is  a  bipartite 
graph.  A  bipartite  graph  is  an  undirected  graph  G  in  which  the  vertices  may  be 
partitioned  into  two  parts  in  such  a  way  that  all  edges  have  a  vertex  in  each  part. 
A  complete  bipartite  graph  is  a  bipartite  graph  having  all  allowed  edges.  Associated 
with  each  (bipartite)  graph  is  a  matching  number  denoted  m  =  m  (G).  The  matching 
number  of  any  graph  is  the  maximum  number  of  vertex  disjoint  edges  of  the  graph. 
In  the  case  of  a  bipartite  graph,  m  (G)  is  less  than  or  equal  to  the  smaller  number 
of  the  vertices  in  its  parts.  In  the  case  where  G  is  a  complete  bipartite  graph,  m  (G) 
will  equal  the  smaller  cardinality  of  the  two  parts. 

Some  answers  are  already  known  regarding  these  problems.  Leal-Duarte  and 
Johnson  [19]  resolved  the  G-SNIEP  for  matrices  subordinate  to  a  bipartite  graph. 
We  will  discuss  their  result  in  the  next  section. 
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8.2  G-SNIEP  for  Matrices  Subordinate  to  a  Bipartite  Graph 

In  the  paper  by  Leal-Duarte  and  Johnson  [19],  we  find  that  simple  linear  inequal¬ 
ities  are  all  that  are  necessary  and  sufficient  to  describe  a  solution  to  the  G-SNIEP 
for  matrices  subordinate  to  a  bipartite  graph.  Their  solution  for  a  —  {Ai,  A2, .  ■  • ,  A„} 
where  Ai  ^  A2  ^  ^  A„  and  Aj  +  A„  ^  0,  is: 

Ai  q-  An  ^  0 

A2  +  ^  0 

Am  ”1”  Ayi— m-t-1  ^  0 
Am+1  ^  0 

An— m  ^  0 

where  m  =  m  (G).  If  n  is  even,  where  m  (G)  =  |n,  the  lower  single  entry  inequalities 
do  not  occur.  The  following  lemmas  walk  us  through  the  proof  given  by  Leal-Duarte 
and  Johnson  that  (23)  gives  a  complete  solution  to  the  G-SNIEP  for  matrices  subor¬ 
dinate  to  a  bipartite  graph.  The  first  lemma  shows  the  necessity  of  these  inequalities, 
while  the  second  lemma  shows  the  sufficiency  for  the  existence  of  a  symmetric  non¬ 
negative  matrix  subordinate  to  G. 

Lemma  2.  [19]  Let  G  be  a  bipartite  graph  on  n  veHices  with  matching  number  m  and 
suppose  that  A  is  a  symmetric  nonnegative  matrix  subordinate  to  G.  The  eigenvalues 
of  A  satisfy  the  inequalities  of  (23). 
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Proof.  Let  A  =  D  +  B  where  is  a  nonnegative  diagonal  matrix  and  5  is  a  non¬ 
negative  symmetric  matrix  with  diagonal  entries  all  zero.  Then  rank  (B)  ^  2m  {G). 
This  is  because  no  submatrix  whose  size  is  larger  than  2m  {G)  can  have  a  nonzero 
term  in  its  determinant.  Also  as  B  is  similar  to  —  B  via  the  signature  matrix  with  1  s 
in  the  diagonal  entries  corresponding  to  one  of  the  parts  of  G  and  -1  s  in  the  other 
diagonal  entries,  we  have 

/^n  ~  /^1 

/in-1  =  -1^2 

/in-m+1  f^m 
/im+l  ® 

l^n—m  ^  0 

in  which  /Ui  ^  >  /in  are  the  eigenvalues  of  B.  But,  since  D  is  positive  semidefinite, 

we  have,  using  the  Weyl  inequalities  [11],  that 

Ai  ^ /ii,  for  i  =  1, . . . ,  n 

in  which  Ai  ^  ^  A„  are  the  eigenvalues  of  A.  It  follows  that  the  Aj’s  satisfy  the 

inequalities  (23).  ^ 

Lemma  3.  [19]  Let  G  be  a  bipartite  graph  on  n  vertices  with  matching  number  m 
and  suppose  that  Ai  ^  ^  A„  satisfy  the  inequalities  of  (23).  Then  there  exists  a 

nonnegative  symmetric  n  x  n  matrix  subordinate  to  G  with  eigenvalues  Ai, . . . ,  An. 
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Proof.  As  we  have  seen  in  Chapter  5,  for  any  two  real  numbers  a  and  ^  such  that 
a  +  P  0  there  is  a  2  x  2  symmetric  nonnegative  matrix  having  eigenvalues  a  and 
p.  Now,  consider  a  collection  of  edges  of  G  that  realize  m{G).  Suppose,  without 
loss  of  generality,  that  they  are  {1,2},  {3,4},  . . .,  {2m  —  1,2m}.  Construct  a  2  x 
2  symmetric  nonnegative  matrix  with  eigenvalues  Ai,A„  and  call  it  Ai,  and  repeat 
this  process  constructing  a  2  x  2  symmetric  nonnegative  matrix  with  eigenvalues 
Afc,  An_fc+i  and  calling  it  A/t  for  fc  =  2, . . . ,  m.  Further  let  the  1x1  matrix  Aj  =  [Aj], 
for  j  =  m  +  1, . . . ,  n  —  m.  Now  Ai  0  •  •  •  0  Am  ®  Am+i  0  •  •  •  0  A^-m  is  a  symmetric 
nonnegative  n  x  n  matrix  with  eigenvalues  Ai, . . . ,  A^  that  is  subordinate  to  G.  □ 

Now  combining  the  two  lemmas  together,  Leal-Duarte  and  Johnson  give  the  simple 
yet  elegant  result. 

Theorem  34.  [19]  Let  G  be  a  bipartite  graph  on  n  vertices  with  matching  number  m. 
There  is  a  symmetric  nonnegative  n  x  n  matrix  subordinate  to  G  and  with  eigenvalues 
Xi,...,Xn  if  and  only  if  Xi,...,Xn  satisfy  the  inequalities  (23). 

We  state  a  few  more  results  derived  by  Leal-Duarte  and  Johnson  from  Theorem  34 
and  the  notion  of  independence  number  of  a  graph.  The  independence  number, 
denoted  i  (G),  of  an  undirected  graph  G  is  the  maximum  number  of  vertices  of  G 
among  which  there  are  no  edges.  In  the  case  of  a  bipartite  graph  G,  i  (G)  =  n—m  (G). 
The  first  result  gives  an  upper  bound  on  the  number  of  negative  and  nonnegative 
eigenvalues  a  symmetric  nonnegative  matrix  subordinate  to  a  bipartite  graph  can 
have.  The  second  result  puts  a  lower  bound  on  the  number  of  nonnegative  eigenvalues 


123 


a  symmetric  nonnegative  matrix  subordinate  to  an  undirected  graph  can  have. 

Corollary  5.  [19]  Annxn  symmetric  nonnegative  matrix  subordinate  to  the  bipartite 
graph  G,  has  at  most  m  {G)  negative  eigenvalues  and  at  least  n-m{G)  nonnegahve 
eigenvalues. 

Theorem  35.  [19]  An  nx  n  symmetric  nonnegative  matrix  subordinate  to  the  undi¬ 
rected  graph  G  has  at  least  i  (G)  nonnegative  eigenvalues. 

Since  this  section  wraps  up  the  solution  to  the  G-SNIEP  for  matrices  subordinate 
to  a  bipartite  graph  and  does  so  with  relative  simplicity,  we  consider  next  the  G- 
NIEP  problem  for  a  bipartite  graph  G.  We’ll  see  in  the  next  section  that  although 
the  G-NIEP  is  a  natural  extension,  a  solution  may  not  have  the  same  simplicity  and 
elegance  that  the  G-SNIEP  solution  does. 

8.3  G-NIEP  for  Matrices  Subordinate  to  a  Bipartite  Graph 

Let  us  consider  a  list  of  complex  numbers  a  =  {Ai,  Aa, . .  • ,  A„}  where  Ai  ^  [Aj]  for 
alH  =  2, . . .  ,n.  We’ll  try  to  find  necessary  and  sufficient  conditions  on  such  a  list  to 
solve  the  G-NIEP  for  trace  zero  nonnegative  matrices  subordinate  to  a  bipartite  graph 
G.  We  already  know  certain  necessary  conditions  on  a  namely:  a  is  closed  under 
complex  conjugation,  the  Perron  root  must  be  an  element  of  a,  the  moments  must 
be  nonnegative,  and  the  JLL  condition  must  hold.  In  addition  to  these,  the  following 
theorem  gives  two  more  necessary  conditions  for  a  spectrum  of  a  non-symmetric, 
nonnegative,  trace  zero  matrix  subordinate  to  a  bipartite  graph. 
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Theorem  36.  Let  Abe  a  non- symmetric,  nonnegative,  trace  zero  n  x  n  matrix  subor¬ 
dinate  to  a  bipartite  graph.  Then  a  (A)  =  {Ai,  A2, . . . ,  A„}  must  satisfy  the  following. 

1.  u  (A)  is  closed  under  multiplication  by  -1  and 

2.  f  n  is  odd,  then  0  G  cr  (A) . 

Proof.  Let  cr  (A)  =  {Ai,  A2,  •  •  ■ ,  A„}  be  the  spectrum  of  a  non-symmetric,  nonnegative, 
trace  zero  n  x  n  matrix,  A,  subordinate  to  a  bipartite  graph.  Then  a  is  a  subset  of 
the  complex  numbers  and  must  contain  the  Perron  root  because  A  is  nonnegative. 

n 

We  also  have  that  Aj  =  0  by  the  trace  zero  assumption. 

As  a  is  the  set  of  roots  of  the  characteristic  equation  for  A  which  has  real  coeffi¬ 
cients,  we  know  it  must  be  closed  under  complex  conjugation. 

We  may  consider  A  to  be  permutation  similar  to  the  nonnegative  matrix  having 

the  form  _ 

Oil  ^12 

A21  O22 

where  On  is  the  A:  x  A:  zero  matrix  and  O22  is  the  (n  -  A:)  x  (n  -  k)  zero  matrix.  This  is 
because  A  is  subordinate  to  a  bipartite  graph  where  we  are  assuming  the  two  vertex 
disjoint  subsets  have  cardinalities  k  and  n  —  k.  So  without  loss  of  generality  let  A  be 
of  this  form.  Then  by  Theorem  6,  we  know  that  the  eigenvalues  of  this  matrbc  are 
the  square  roots  of  the  eigenvalues  of  the  product  A12A21  which  is  a  A:  x  A:  matrix. 
This  means  we  have  at  most  2k  nonzero  eigenvalues  where  the  remaining  n  —  2k 
eigenvalues  must  be  zero  for  A:  <  f .  Now  since  the  product  A12A21  is  a  nonnegative 
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matrix,  we  have  the  spectrum  of  ^12^21  to  be  closed  under  complex  conjugation.  So 
for  re*®  6  cr  {A12A21),  it  follows  that  re~*®  €  a  (^12^21)-  Therefore  the  square  roots  of 
these  two  numbers  will  be  in  a  (A).  That  is,  iy^e  2  6  cr  (A)  and  -iz^/re  2  g  cr  (A). 
As  this  is  at  most  four  eigenvalues  and  contains  the  additive  inverse  of  each  one,  we 
find  that  a  (A)  is  closed  under  multiplication  by  -1. 

Now  due  to  the  restriction  that  cr  (A)  is  closed  under  multiplication  by  -1,  we  find 
that  the  natural  pairing  of  G  a  (A)  with  -A^  G  cr  (A)  can  be  done.  So  if  n  is  odd, 
trace  (A)  =  0  requires  that  the  remaining  unpaired  eigenvalues  must  be  zero.  Thus 

0  G  <7  (A).  ^ 

The  next  result  allows  us  to  construct  a  matrix  having  a  =  {±v^>  •  •  •) 
when  we  know  of  a  nonnegative  matrix  having  Ai, . . . ,  A„,  as  its  eigenvalues. 


Theorem  37.  Let  Abe  a  nonnegative  n  yin  matrix  having  eigenvalues  Ai,  A2, . . . ,  An- 
Then  the  matrix  B  defined  by 

[0  A 


B  = 


I  0 


5 


where  I  is  the  n  x  n  identity  matrix,  is  a  nonnegative,  trace  zero,  2n  x  2n  matrix 


subordinate  to  a  bipartite  graph  with  matching  number  n.  Moreover,  the  eigenvalues 
of  B  are  ±^/X^,  ±>/X2,  . . .,  ±\/A^- 


Proof.  Using  Theorem  6,  let  A21  =  I  and  A12  =  A  where  each  one  is  an  n  x  n  matrix 
.  Then  the  eigenvalues  of  A21A12  =-IA  =  A  are  indeed  the  eigenvalues  of  A.  Then  we 
find  by  Theorem  6,  that  the  eigenvalues  of  B  are  indeed  ±^A^,  ±\/^,  •  •  •>  ° 
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Theorem  36  and  Corollary  37  are  used  for  solving  the  G-NIEP  subordinate  to 
a  bipartite  graph  in  the  case  where  a  =  {Ai,  A2, . . . ,  An}  contains  seven  or  fewer 
elements.  We  can  construct  solutions  to  these  cases  using  the  results  from  prior 
chapters. 

Let  us  consider  the  following  cases  for  the  G-NIEP  for  trace  zero  matrices  subor¬ 
dinate  to  a  bipartite  graph.  Let  Ai  be  the  Perron  root  in  all  cases. 


1.  In  the  case  where  cr  =  {Ai},  the  only  solution  is  when  Ai  =  0.  Then  such  a 
nonnegative,  trace  zero  matrix  is  the  zero  matrix,  and  it  is  subordinate  to  a 
bipartite  graph  with  one  vertex  and  no  edges. 


2.  Consider  the  case  where  a  =  {Ai,A2}.  First  of  all,  to  satisfy  Theorem  36, 
A2  =  -Ai.  This  means  A2  is  real  since  Ai  is.  The  solution  to  such  a  set 


cr=  {Ai,-Ai}  is: 


A2 


0  A? 

1  0 


This  matrix  has  eigenvalues  i:Ai  and  is  subordinate  to  a  bipartite  graph  on  two 
vertices. 


3.  Consider  the  case  where  cr  =  {Ai,  A2,  A3}.  As  we  know,  since  Ai  G  <7  it  follows 
that  -Ai  e  a  too.  Without  loss  of  generality,  let  A2  =  -Ai.  Since  the  trace  is 
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zero,  we  have  A3  =  0.  Thus  using  the  result  from  (2),  we  have 


A?,  — 


0  A?  0 
1  0  0 
0  0  0 


This  is  a  nonnegative  matrix  subordinate  to  a  bipartite  graph  having  three 
vertices  with  matching  number  1. 


4.  Consider  the  case  where  a  =  {Ai,  A2,  A3,  A4}.  We  may  assume  A2  =  — Ai  and 
A4  =  -A3.  Note  that  A3  could  be  a  pure  imaginary  number  or  a  real  number. 
Since  A3  could  be  negative,  a  necessary  requirement  for  a  2  x  2  matrix  with 
eigenvalues  A?,  A|  to  exist  is  that  A?  +  A^  ^  0.  But  this  requirement  follows 
from  the  necessary  condition  on  a  that  Ai  ^  IA3I. 


Therefore  using 

Aj  —  Ag  Aj  +  Ag 

we  find  a  nonnegative  matrix  subordinate  to  a  bipartite  graph  on  four  vertices 
with  matching  number  2  using  Theorem  37,  namely 


A4  — 


0  B 
I  0 


5.  Consider  the  case  where  a  =  {Ai,  A2,  A3,  A4,  A5}.  Using  Theorem  36  we  know 
that  one  of  the  values  in  a  must  be  0,  so  assume  A5  =  0.  Then  the  remaining 
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numbers  must  satisfy  the  conditions  from  (4).  Let  A.4  be  the  nonnegative  ma¬ 
trix  having  eigenvalues  Ai,  A2,  A3,  A4.  Then  the  following  nonnegative  matrix  is 
subordinate  to  a  bipartite  graph  on  five  vertices  with  matching  number  2  and 
has  eigenvalues  Ai,  A2,  A3,  A4,  A5. 


^5  — 


0 


0 

0 


6.  Consider  the  case  where  a  =  {Ai,  A2,A3,A4,  AsjAe}.  This  is  the  first  time  a 
complex  number  with  nonzero  real  part  and  nonzero  imaginary  part,  may  be 
an  element  of  cr.  This  is  because  cr  must  be  closed  under  complex  conjugation 
and  additive  inverse,  and  have  the  Perron  root  as  an  element.  For  notation 
purposes,  let  Aj  and  Aj+i  for  i  odd,  be  the  additive  inverses  of  each  other. 

Therefore  we  can  use  the  3x3  symmetric  nonnegative  solution  from  Theorem  9 
to  construct  the  nonnegative  matrix  A  having  the  eigenvalues  Ai,A3,A5.  Then 
by  Theorem  37,  we  have  the  nonnegative  matrix  subordinate  to  a  bipartite 
graph  on  six  vertices  with  matching  number  3.  For  example. 


0  0  0  2  3  4 
0  0  0  6  7  8 
0  0  0  4  6  7 
1  0  0  0  0  0 
0  1  0  0  0  0 
0  0  1  0  0  0 
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is  a  matrix  subordinate  to  the  complete  bipartite  graph 
a  (B)  =  {±4.03997,  ±  (0.408884  -  0.572597i) ,  ±  (0.408884  +  0.572597i)}. 

7.  Consider  the  case  where  a  =  (Ai,  A2,  A3,  A4,  A5,  Ae,  A7}.  By  Theorem  36  as  seven 
is  odd,  we  know  that  one  of  the  eigenvalues  of  a  must  be  zero.  So  we  may  use 
the  construction  technique  of  (6)  to  create  a  nonnegative  matrix,  B,  having 
spectrum  a  —  {0}.  Then  our  solution  for  a  nonnegative  matrix  subordinate  to 
a  bipartite  graph  on  seven  vertices  with  matching  number  3  is 

B  0 
0  0 


8.4  Conclusion 

While  this  chapter  gives  a  partial  solution  to  the  G-NIEP  for  bipartite  graphs, 
it  becomes  apparent  that  continuing  on  in  this  fashion  requires  the  solution  to  the 
NIEP  in  the  general  setting. 
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